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1.0 INTRODUCTION AND SUWARY OF RESULTS 


1.1 INTRODUCTION 

This report discusses the results of a study completed by personnel of 
the Electro-Optical and Data Systems Group, Hughes Aircraft Company In El 
Segundo, CA for NASA/GSFC in Greenbelt, Md under contract NAS 5-27139. The 
intent of the study was to determine the mutual influence of a laser communi¬ 
cation system and Its host spacecraft and to what degree the mutual influence 
limited the acquisition, tracking and pointing processes so important to a 
laser communication link between thn two spacecraft. The study used as a 
baseline a laser link between a low earth orbiting (LEO) satellite and a 
geosynchronous earth orbiting (GEO satellite. The LEO satellite was modeled 
as a LANDSAT satellite and the GEO as a TORSS satellite. Vibrational and 
disburbance data from the LANDSAT satellite was used to model the baseframe 
noise environments for the study* 

The laser link between satellites was a generic channel transferring 500 
M>ps data from the LEO to GEO using the GaAIAs laser as the laser light 
source. The LEO optical subsystem was fashioned as a 6 inch diameter optical 
telescope for transmit and the GEO used a 20 inch diameter optical telescope 
for receive functions. The probability of bit error was specified as 1 x 
10 . The bit probability of burst error due to mispointing was not 
specified but was shown by the study to be of basic concern. Burst errors 
arise when the realtime mispoint angle is greater than the nominal mispoint 
angle. The occurrance of such mispoint is described by a orobability function 
and has been evaluated herein. 

The major aspects of pointing and tracking with a satellite-borne optical 
system have been evaluated. We have evaluated the major influencing areas. 
Including: (1) orbital aspects such as spacecraft relative motions, point 
ahead, and sun and moon optical noise, (2) burst errors introduced by the 
electronic and optical noise levels, (3) servo system design and configura¬ 
tions, and the noise sources such as, sensor noise, base motion disturbances. 
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gimbal friction torque noise, (4) an evaluation of the tracking and beacon 
link and the type of sensors used, (5) the function of the acquisition 
procedure and an evaluation of the sensors employed, and (6) an estimate of 
the size, wight and power needed for the satellite system. The system was 
evaluated as a two axis gimballed configuration with all optical, electro- 
optical, laser and receiver components borne on-gimbal. Support electronics 
such as servo control electronics, communication electronics, power supply and 
thermal control were placed off-gimbal. The laser source was a GaAlAs laser 
diode operating at about 8300 A. An avalanche photodiode communication sensor 
was assumed. 

The major study results are briefy summarized in Table 1.1-1. The study 
indicated that 1 urad rms pointing may be achieved on a satellite with base- 
frame motion disturbances as measured on the LANDSAT satellite. The design of 
the servo control subsystem is critical and those aspects have been identi¬ 
fied. The system weight and power values have been estimated for a non- 
redundant configuration and for redundant electro-optical and electronic 
components. The redundant system weight for a LEO satellite is about 175.2 lb 
and that for a GEO satellite system is 211.5 lb. The associated electrical 
power requirements are 137.8 watt and 143.1 watt, respectively. Since the 
system incorporates moving electro-optical and optomechanical elements, the 
torque imparted to the satellite is of concern. We estimate that the torque 
will be 0.03 ft lb. Spatial acquisition bewteen the LEO and GEO with an 
angular uncertainty of 0.5 degree may be accomplished within 20 minutes if a 
solid state array sensor is employed. 

In general, we find that the pointing/tracking and acquisition require¬ 
ments can be met with an appropriately design system, and with components that 
are currently commerically available. 

1.2 ORBITAL ASPECTS 

The orbital aspects of the two satellite configuration was evaluated with 
respect to the influence upon the operation of the laser communication system. 
The laser system was modeled as a two axis gimballed configuration and it was 
determined that the gimbal would have requirements as summarized in Table 
1 . 2 - 1 . 
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TABLE 1.1-1. MAJOR STUDY RESULTS 


Parameter 

Study Result 

• 

Pointing Accurary 

• 

1 urad rms possible with 
properly designed servo 
control subsystem, with LANDSAT 
type base motion disturbances 

• 

Weight (with redundancy) 

• 

LEO * 175.2 lb 



• 

GEO * 211.5 lb 



• 

Redundant components assumed 

t 

Power 

• 

LEO s 137.8 



• 

GEO = 143.1 

• 

Torque imparted to 
spacecraft by moving 
optomechanical elements 

• 

0.03 ft-lb 

• 

Acquisition time 

• 

2 minutes, with area array 
sensor 
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TABLE 1.2-1. SUMMARY OF GIMBAL REQUIREMENTS (MAXIMUM VALUES) 

DETERMINED BY THE LEO AND GEO ORBITS. 


Pointing Error 

i 

Angular Velocity 

Angular Acceleration 

Spacecraft 

Az 

El 

Az 

EL 

Az 

El 

LEO 

♦ 175* 

♦ 140* 

200 urad 

sec 

200 urad 

sec 

a 

! 

t 

1 urad 
sec* 

0.5 urad 
sec^ 

GEO 

♦ 20* 

♦ 20* 

1 u^ad 

sec 

200 urad 

sec 

1 urad 
sec* 

1 urad 
sec^ 


The various parameters associated with pointing and tracking from one terminal 
to the other were graphically illustrated. These parameters included; the 
pointing angles in azimuth (Az) and elevation (El), the pointing angle rates 
(Az-dot and EL-dot), the angular acceleration (Az-double dot and El-double 
dot), and point ahead angles. One example of this type of orbital calculation 
is shown in Figure 1.2-1. 

The influence of the sun and moon as inhibiting agents to maintenance of 
the satellite to satellite link due to being high background optical niose 
sources was also evaluated. The angular position of the sun and moon relative 
to the line of sight between the two satellites was determined with the 
orbital program. The amount of operational time lost per month when the sun 

w 

and moon were within a small angle of the line of sight was evaluated. The 
angles were 4*, 1* and 100 urad. The system outage was noted graphically and 
in tabular form In hours per month over the year. An example Is shown in 
Figure 1.2-2. 
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The clear line of sight between satellites was determined to be about 71 
percent of the time as about 29 percent of the time the LEO satellite Is 
obscurred by the earth from the GEO satellite. The amount of time loss per 
month due to glmbal lock was also determined. This Is principally due to 
glmbal lock In the LANDSAT satellite. It was found, for example, that the 
maximum amount of time lost due to the sun vs moon In 4* field of view along 
the line of sight is less than a maximum of about 166 minutes per month with 
an average of only about 60 minutes per month. The maximum time lost per 
month due to glmbal dock was about 3300 minutes and the average was about 2600 
minutes. 

1.3 SYSTEM CONFIGURATION 

The optical and mechanical configuration of the optical communication 
terminal for the LEO and GEO satellites was developed. Two concepts were 
evaluated. One used a flat scan mirror to direct the laser light to the 
distant receiver. Common transmit and receive optics were assumed. The 
second used separate transmit and receive optics to avoid the complexity of 
dichrolc elements. All optics and laser related components are borne on a two 
axis glmbal. This latter approach Is the baseline concept. The selected 
configuration is noted in Figure 1.3-1. 

The basic system block diagram (transmitter) Is noted In Figure 1.3-2. 
The figure notes at the laser transmitted, receiver, associated electronics 
and all optics are borne on-gimbal. The large scale, low frequency pointing 
and tracking errors are corrected by the two axis gimbal. The small, high 
frequency angular disturbances are correct by the small image motion compensa¬ 
tion (IMC) mirrors. The normal processing and support electronics are placed 
off glmbal. A blfuncating mirror Is used to provide the separation between 
the large acquisition field of view and the smaller tracking field of view. A 
6 Inch telescope Is used In the transmitter and a 20 inch telescope Is used In 
the receiver. 
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Figure 1.3-1. Baseline Configuration (LEO and GEO). 
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Figure 1.3-2. Basic System Block Diagram (transmitter). 






The design parameters of the baseline approach are noted in Table 1.3-1 
and the concept design functions are noted in Table 1.3-2. The size, weight 
and power summary for the LEO terminal is a 6 inch diameter telescope, 175.2 
pounds and 137.8 watts. The associated factors for the GEO terminal is a 20 
inch diameter telescope, 211.5 pounds and 143.1 watts. 


TABLE 1.3-1. TERMINAL CONCEPT DESIGN PARAMETERS 



Property 

GEO 

1 LEO 

• 

Concept 

• Separate Transmit 
and Receive Optics 

• Separate Transmit 
and Receive Optics 

• 

Transmit Optics 

• 2.3 Inch Diameter 

• 6 Inch Diameter 

• 

Receive Optics 

• 20 Inch Diameter 

• 6 Inch Diameter 

• 

Primary Material 

• Beryllium 

• Beryl 11 urn 

• 

Acq Sensor 

• SI-Quad 

• Si-Array 

1 

! • 

Track Sensor 

• Si-Quad (APD) 

t Si-Quad (APD) 

• 

Comm Sensor 

• SI-APD 

t SI-APD 

t 

Laser Source 

• GaAlAs 

• GaAlAS 

• 

Transmit Wavelength 

• Beacon (LDR) = 0.85 pm 

• Comm (HDR) ■ 0.83 pm 

• 

Receive Wavelength 

• Comm (HDR) ■ 0.83 pm 

• Beacon (LDR) ■ 0.85 pm 

• 

Gimbal 

• Two Axis-Optics on 

• Two Axis - Optics on 


4 

Gimbal 

Gimbal 

• 

Az Range 

• ♦ 20* 

• + 175* 

• 

El Range 

• 1 20 * 

• ♦ 140* 

• 

Weight (with 
redundancy) 

211.5 1b 

175.2 lb 

t 

Power 

• 143.1 watt 

• 137.8 watt 









TABLE 1.3-2. CONCEPT DESIGN FUNCTION DEFINITIONS 


Function 

Design 

Parameter 

Acquisition 

1 

1 

1 

i 

t Concept 

• Fleld-of-Regard 

• Angular Scan Rate 

• Laser Wavelength 

• Laser Power 

• Beamwldth 

• Acquisition Sensor 

• Acquisition Time 

• Initiated by GEO 

t 0.5 deg Cone for GEO and LEO 
e 1.6 deg/sec 

• 0.85 um 

t 50 mW - cw 

• 130 urad (full Airy Lobe) 

• Nominal 400 x 500 area array 
on LEO 

Quad Cell on GEO 

• £ 2 minutes 

i 

Tracking 

i 

i 

i 

• 

• Concept 

• Field-of-Vlew 

• T C Wavelength 

• T C Data Rate 
t GEO Beacon Power 

• GEO Beacon Beamwldth 

• GEO and LEO Sensor 

• GEO Tracks on Data Uplink 

• LEO Tracks on T C Downlink 

• 140 prad full Sensor, 

13.7 urad Closed Loop 

e 0.81 um 

• 10 Kbps 

• 20 mW - cw 

• 35 urad (full Airy Lobe) 

. • APD - Quad Cel1 

Communications 

Point Ahead 

i 

• Concept 

• Power 

• Beamwldth 

• Data Rate 

• Wavelength 

t Concept 

• Multi source Summing Device 

• 200 mW - Average 

• 10 urad (1/E 2 )* 

• 500 MBPS * 

• 0.8 - 0.9 um 

• Closed Loop Update via T $ C 
Link for LEO Uplink 

• Minor Update for GEO Downlink 

* Not consistent with burst error requirements and 1 urad pointing capability 
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1.4 BURST ERRORS 


The performance of a free space, satellite-borne optical communication 
system Is describe! by two major parameters: the probability of bit error and 
the probability of burst error. The probability of bit error is normally 
specified as a system requirement. It determines the nominal signal Intensity 
required to recover the signal from the communication channel with its con¬ 
tributory noises. The nominal probability of bit error Is proportional to the 
signal and noise counts received during the bit time interval. The probability 
of burst error on the other hand, Is frequently Ignored as It is not directly 
dependent upon the optical signal and noise levels. The probability of burst 
error establishes the percent of time that the servo system will maintain 
pointing within some error limit. Should the pointing error c become greater 
than some value c*, the signal intensity at the distant receiver would be 
decreased, and the system performance described by the probability of bit 
error would be degraded. Many data bits could be lost since the pointing 
servo response time is much longer than the bit time. 

The single axis probability distribution for mispointing (P *) Is 
described as a gaussian function. It is shown that the gausslan assumption is 
reasonable for the complementary filter servo approach as sensor noise 
(gaussian In form) dominates the servo control system performance. The 
probability of burst error Is defined as 



where the required pointing error must be maintained at e < c*. 

The resulting equation for the burst error Is related to the laser 
beamwidth by the standard range equations. The nominal mispoint loss which 
should be included In the link budget Is shown to be about -4 dB. The 
resulting equation may be manipulated to show that the diameter of the optical 
antenna is limited by the burst error. This Is given by 
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X 

D * —— 

2.4 0 sj - In P* 

where x Is the laser wavelength and o Is the standard deviation of the single 
axis gausslan probability distribution. 

The diameter of the optical antenna used to provide the optlmlum beam- 
width In the pointing err>r environment described by the standard deviation o 
Is related to the probability of burst error as shown In Figure 1.4-1. 
Examples of the system requirements under burst error limitations are noted In 
Table 1.4-1. 
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Table 1.4-1. EXAMPLE SYSTEM REQUIREMENTS UNDER 

BURST ERROR LIMITATIONS. 


1 

System 

Availability 

(percent) 

PBE 

t/o 

a/ o 

90 

io-i 

2.13 

8.86 


10-2 

3.03 

12.60 


99.9 

1 

10-3 

3.70 

15.39 

99.99 

io - 4 

4.28 

17.80 

99.999 

10-5 

4.79 

19.93 

99.9999 

10-5 

i 

5.25 

21.84 

1 


Optimum Antenna 
Diameter for 
o * 1 urad and 
x = 0.83 nm 
(inch) 


10.2 

9.3 


Relative 
Power 
Required 
(dB) 


3.1 

4.8 

6.1 

7.0 

7.8 


& 









1.5 SERVO ANALYSIS 


A servo analysis of the baseline system configuration was completed. 
Both azimuth and elevation axes were represented by a single axis analysis. 
All error sources were assumed identical for both axes except for base 
motion. Spacecraft base motion coupling was in accordance with the LANDSAT 
satellite disturbance data provided by NASA/GSFC. The worst case envelope of 
all the spacecraft disturbances were used for the study. 

Three servo configurations were considered for the study and included; 
(a) gyro-stabilized, (b) mass stabilized and (c) complementary filter. The 
gyro stabilized configuration has an outer low bandwidth track loop with two 
higher bandwidth inner loops. The mass stabilized design replaces the two 
gyros with a reference platform that points along the desired LOS. The 
complementary filter configuration closes a single wide bandwidth loop about 
the track sensor with no inner loops. Two bandwidths were evaluated; 1 Hz for 
the gyro stabilized servo and 400 Hz for the complementary filter servo. 

These were the contributing noise sources from the servo analysis; base 
motion, gimbal friction and track sensor noise. Other noise sources were 
shown to be neglibible or may be removed by appropriate calibration. Sensor 
noise was represented as white gaussian noise and its contribution to the 
residual pointing error evaluated for each servo configuration. Similarly, 
the contribution by base motion disturbances was evaluated. The measured 
satellite LANDSAT base motion dlsburbance provided by NASA/GSFC were used for 
the evaluation. The measured disturbances from 0 to 125 Hz were converted to 
a power spectral density (PSD) curve by squaring the amplitude at each 
frequency and a 1 Hz bandwidth was assumed. 

The friction disturbance was modeled as a Coulomb friction torque acting 
at the gimbal bearings. Friction and motor hysteresis were not addressed as 
they will be dependent upon the particular devices selected for the satellite 
system. A PSD for the friction disturbance was developed from a probability 
density function (PDF) which was derived from modeling the Coulomb friction as 
as a random process. THe PSD was applied to the three configurations and the 
effect documented. 


* 
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A sumnary of the contributions due to the three main error sources is 
noted in Table 1.5-1. The contributions due to four other error sources were 
also evaluated but noted to be less significant. These included; point ahead 
error, dynamic following error, linear variable displacement transformer 
(LVDT) non-linearity error and gyro noise errors. The resulting data were 
included with the other results for the single axis evaluation and evaluated 
for a two axis (Az, El) configuration. 


TABLE 1.5-1. TWO AXIS ERROR SUMMARY 


Servo 

Aperture 

Mean LOS 1 

la LOS 

Total LOS 

Configuration 

Diameter 

Error (urad) 

Error (urad) 

Error {urad) 

Gyro Stab 

6 inch 

0.078 

0.63 

0.71 


20 inch 

0.078 

0.49 

0.57 


6 inch 

0.078 

11.65 

11.73 

| Mass Stab 

20 inch 

0.078 

11.65 

11.73 

Compl. 

Fi1 ter 

6 inch 

0.066 

0.93 

1 

1.0 

« 

20 inch 

0.066 

0.52 

0.59 


Two conclusions were reached as a result of the servo analysis. First, 
it was determined that a 1-nrad pointing and tracking accuracy could be 
maintained with either the gryo-stabilized or complementary filter approaches. 
Although the gyro-stabilized approach had better tracking performance and 
lower sample rates from the optical sensor, a space-qualified gyro and greater 
number of control loops are required. The second conclusion is to recommend 
the complementary filter approach since it is the simplest design that meets 
the pointing requirements. 
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1.6 TRACKING AND BEACON LINK 


The tracking and beacon link between the LEO satellite and the GEO 
satellite was evaluated. The pointing accuracy required for the LEO satellite 
is more severe than that required for the GEO link since the LEO terminal has 
a narrower beamwidth and is the high data rate transmitter. The study focused 
on the LEO link, recognizing that if it may be accomplished, the GEO return 
link could easily be maintained. 

This analysis focused on the problem of tracking a data beacon from the 
GEO satellite to less than one microradian rms accuracy in order to derive 
pointing commands for a high data rate LEO transmitter. The noise equivalent 
angle of a quadrant tracking detector was derived explicitly stating the 
assumed signal and noise parameters per element. The tracking detector 
quadrant gap has been identified as an area needing further technology 
development, but which does not preclude operational systems with currently 
available devices. The angular dynamic range over which closed loop tracking 
can be maintained was shown to be proportional to the rms SNR in an element of 
the quadrant detector. 

Photodetection analysis and link equations were presented to allow 
calculation of tracking (and acquisition) link budgets. It was shown that a 
PIN quadrant detector was not sensitive enough to close the tracking link 
budget to 1.0 prad rms accuracy unless a gyro-stabilized tracking servo is 
used or laser of 200 mW or greater are available. The quadrant APO detector 
is about 19 dB more sensitive than the PIN, which allows 1.0 prad tracking 
accuracy with the complementary filter servo approach, 35 urad GEO beam- 
width, a commercial 20 mW laser diode, and 6 dB link margin. 

For tracking a GEO 10 Kbps beacon data link three different configura¬ 
tions were shown to be viable. A 50-50 optical beamsplit between an APO 
quadrant tracking detector and single element APD data detector satisfied both 
functions (with adequate margin) with a commercial 20 mW laser. Or the 
quadrant APO tracking detector may also provide the data detection by summing 
the four outputs. This may be necessary if the tracking sample rate has 
significant spectral overlap with the beacon data, as would be the case with a 
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4.0 KHz complementary filter sample rate and 10 Kbps beacon data. Finally, It 
was shown that If an array were selected as the acquisition detector, It could 
also perform the tracking function. A CCD (or CIO) array with standard video 
frame rates is compatible with the gyro-stabilized tracking servo. A CID 
array with a circuit to create a 40 x 40 pixel tracking gate could provide 

updates at the 4.0 KHz complementary filter servo sample rate. 

* 

1.7 ACQUISITION 

Acquisition for optical coirmunications is the process by which two 

terminals initially locate one another and progress to mutual fine tracking 

* 

and two way communication. Each satellite terminal has some knowledge of the 
distant location of the other in terms of angular data relative to an inertial 
coordinate system. Its position is then an uncertainty solid angle which must 
be interrogated to solicit response from the opposite terminal. 

Several desirable features of an acquisition system were early 
established. These included: (1) minimum number of acquisition peculiar 
components and complexity, (2) reliable acquisition within several minutes of 
acquisition start, (3) single stage transfer from wide field of view (FOV) 
acquisition to narrow field of view tracking and (4) avoidance of mechanically 
moving elements in the optical path to accomplish the FOV reduction. 

Spatial acquisition between the LEO and GEO is critical as it must be 
accomplished first before communications may begin. Additionally, it must be 
accomplished periodically as the LEO to GEO line of sight is broken as the LEO 
passes behind the earth in its orbit. In general, the positional uncertainty 
of one satellite from the other is dependent upon several factors. These 
include: (1) the altitude uncertainty of the satellite itself, (2) orbital 
ephemeris uncertainties of the position and velocity of the other satellite 
from the first, (3) relative timing errors and point ahead errors. The 
relative significance of these errors was evaluated with the use of actual 
TDRSS satellite data. It was found that satellite attitude errors dominate 
and could be as. large as about ^ 0.23 deg, 3 axis, rms. We therefore used a 
value of 0.5 degree diameter for the angular uncertainty cone- from one 
satellite to the other at the beginning of acquisition. 
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The key parameters associated with acquisition performance, once the 

» 

angular uncertainty is known, is the acquisition time (the time from Inltlaton 
to fine tracking), the laser power required to accomplish acquisition, and the 
probability that successful acqulslton takes place. The acquisition sensors 
and modulation/detection alternatives were evaluated with emphasis on these 
parameters. 

Acquisition sensors were Investigated, Including CCD, CID and photodiode 
arrays, as well as PIN and APD quadrant detectors. Commercially available 
arrays seem to be compatible with acquisition requirements, while quadrant 
detectors may require development in the areas of smaller quadrant gap, and 
smaller device size to reduce capacitance and dark current. 

Acquisition sources were evaluated. Including the single mode cw lasers 
and the high peak power Laser Diode Labs incoherent array. Although the 1.0 
kW peak power of the LDL array was appealing, it was found that Its large 
source size and divergence made It difficult to collect and collimate Its 
emission. Single mode cw lasers (or coherent multi-stripe arrays, if they 
become available) are recommended as the choice for acquisition source. 

Issues for quadrant detector acquisition were examined, including field 
of view reduction and servo angular step reduction. The probability of 
successful acquisition was derived based on detection probability when 
illuminated by the transmitter, and false alarm probabilities while not 
illuminated. Three different modulation/detection approaches for acquisition 
were evaluated for quadrant detector acquisition. These are the cw Inte¬ 
grating (CWI) approach, sinusoidal Intensity modulation (SIM), and time-gated 
acquisition (TGA). The CWI approach was found to yield the best shot noise 
performance for quadrant detector acquisition, but Is seriously limited by the 
DC component of optical background in a practical system. SIM avoids the DC 
background problem, but Is about 4.5 dB worse In shot noise performance than 
CWI. Acquisition times for SIM exceed ten minutes with a quadrant APD 
detector and a 100 mW source, with no power margin, the TGA scheme was found 
to have no advantage over SIM. The quadrant detector acquisition sequence was 
Identified. 
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The field of view reduction for array acquisition was derived, as well as 
the probability of a successful acquisition. Commercially available CCD 
arrays were found to allow acquisition times under two minutes, with a 50 mW 
source and 6.0 dB power margin. The acquisition sequence using an array was 
described. 

Acquisition times were found to be directly proportional to the 
uncertainty solid angle, and Inversely proportional to laser power, for both 
array and quadrant detector acquisition. 

Issues for scanning the acquisition transmitter beam were identified. 
Including optimum beamwldth, beam overlap, scanning with a pulsed laser, and 
spacecraft uncompensated drift. 

A recommended acquisition/tracking configuration was proposed, which used 
an acquisition array on the LEO, and a quadrant APD for both tracking and 
beacon data demodulation on the LEO. Further investigation Is recommended to 
evaluate performing the tracking funtion with the acquisition array and a 
video tracker. 

1.8 SYSTEM WEIGHT AND POWERS 

Weight and power estimates were made for the baseline concept - the two 
axis gimballed optical assembly. The estimate rationale was based on the use 
of beryllium as well as light weighted cervlt for the material for the main 
optical elements. Beryllium optics is baseline. The weight and power 
estimates were made for both the LEO and GEO satellites. Redundancy of the 
electronic, electro-optical and some optical components was considered as 
necessary for a reliable system. A summary of the results are noted In Table 
1 . 8 - 1 . 

TABLE 1.8-1. WEIGHT AND POWER ESTIMATE 



Weight (lb) 

Power 

Satellite 

with Reduncancy 

(watt) 

LEO 

175.2 

137.8 

GEO 

211.5 

143.1 








2.0 REQUIREMENTS 


2.1 STUDY TASKS 

* 

The study tasks were three fold. Task 1 - determine the Influence of the 
spacecraft on the communication system. The Spacecraft can Introduce link 
degradation due to the Influence of attitude and ephemerls errors on the LEO 
to GEO and GEO to LEO pointing accuracy. Additionally, the spacecraft can 
Influence the link acquisition time and the complexity of the acquisition and 
tracking subsystems by Introducing pointing or tracking errors. Task 2 - 
determine the Influence of the communication system on the spacecraft. 
Torques will be generated by moving optical elements and the torques must be 
recognized to allow proper design of the spacecraft control system. Task 3 - 
develop a concept design for the LEO and GEO termihals. 

The task requirements are summarized in Table 2.1-1. 

2.2 LINK REQUIREMENTS 

* 

The communication system link requirements are summarized In Table 
2.2-1. The system is configured as a 500 megabit per second (Mbps) channel 
originating at the LEO (Landsat) satellite and relayed to the GEO (TDRSS) 
satellite. The laser source Is a GaAlAs laser diode operating at a nominal 
wavelength of 0.8 urn. The probability of bit error was specified as 1 x 
10 . Burst error was not specified but its Influence Is discussed in 
section 4.0. A 6 inch diameter telescope is specified for the LEO and a 20 
inch diameter telescope is specified for the GEO. It Is shown in section 4.0 
that these diameters are consistent with a pointing error of 1 urad. 














TABLE 2.1-1. STUDY TASKS 



Spacecraft Impact on the 






System Impact on the 



System Concept Design 


Specifics 


(a) Determine the Influence on the performance 
or the communication link due to spacecraft 
related errors: 

• Attltude errors 

• Attitude Instabilities 

t Ephermls errors 

• Attitude commands uplinked from ground 

control 

(b) Determine the Influence on the system per¬ 
formance due to spacecraft related errors: 

• Acquisition time 

t Tracking accuracy 

• Pointing accuracy 


(c) 


(a) 



(b) 


(a) 



Determine the system configuration anc 
complexity due to spacecraft related errors: 

• Number of acquisition stages 

• Number of track stages 


Determine the Influence on the spacecraft 
due to the operation of the communication 
system during acquisition, tracking, 
pointing and conrnunlcatlon: 

• Moving optical elements 

e Moving mechanical elements 

Determine the feedback to the communication 
system from the spacecraft originally Influenced 
by the operation of the system. 


Develop concept designs for the LEO and GEO 
terminals 


Configuration 
Block diagram 
Weight 
Size 
Power 













TABLE 2.2-1. SYSTEM LINK REQUIREMENTS 


ITEM 

REQUIREMENT 

LEO 

• LANDSAT 

GEO 

• TDRSS 

DATA RATE 

• 500 MPBS 

PBE 

• 1 x lO" 6 

DT 

• 6 INCH 

OR 

t 20 INCH 

WAVELENGTH 

• 0.80 Min 















3.0 ORBITAL ASPECTS 


3.1 INThCAUCTION 

The optical communication system spacecraft architecture under evaluation 
Is a Landsat to TORSS link. The primary high data rate communication terminal 
Is assumed to be borne on the low earth orbiting (LEO) LANDSAT spacecraft and 
the high data rate receiver borne In the geosynchronous earth orbiting (GEO) 
TDRSS spacecraft. The orbital parameters for the two spacecraft are noted In 
Table 3.1-1 and the orbital configuration noted In Figure 3.1-1. 


TABLE 3.1-1. ORBITAL ASSUMPTIONS 


SPACECRAFT 

j 

PRIMARY TERMINAL 

ORBIT 

Landsat 

Transmitter 

• low earth orbit 


• 

e circular orbit 

e 705.3 Km altitude 

# 98.2* inclination 

• sun-synchronous 

TDRSS 

Receiver 

• geosynchronous orbit 

• 43,900 Km altitude 

• 41* W position 

! 


The spacecraft gimbal coordinate system Is defined in Figure 3.1-2. The 
coordinate system Is applicable to both the low earth orbit (LEO) and 
geosynchronous earth orbit (GEO) satellites. The X axis Is nadir pointing, Y 
Is opposite the velocity vector In the orbital plane, and Z is perpendicular 
to the orbit plane following a right hand rule. 
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CIRCULAR ORBIT 



Figure 3.1-1. Laser Crosslink Configuration. 









AXIS 

DEFINITION 

A 

X 

NADIR POINTING 


A 

Y 

OPPOSITF VELOCITY VECTOR 


A 

Z 

A A 

PERPENDICULAR TO X AND Y, RH-RULE 


AZ 

A A 

IN ORBITAL PLANE, POSITIVE FROM X TO Y 


EL 

A A 

OUT OF ORBITAL PLANE NEGATIVE FROM X TO Z 


Figure 3.1-2. Spacecraft Gimbal Coordinate System. 
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3.2 CROSS LINK PARAMETERS 


The applicable link parameters have been defined for the LEO (Landsat) to 
GEO (TDRSS) and GEO (TDRSS) to LEO (Landsat) configuration using an orbital 
program previously developed. The program determines the pointing angles for 
azimuth (Az) and elevation (El) for each satellite terminal, the pointing 
angle rates for Az and El, pointing angle accelerations for Az and El, and 
point ahead angles in Az and El. An example of the expected Az and El 
pointing angles (with the GEO designated as 1 and the LEO designated as 2) for 
the monthof September are shewn in Figure 3.2-1. In general, the Az gimbal 
coverage required for the GEO satellite is about ^20* Az and about +20* El. 
That required for the LEO satellite is about ^175 Az and +140 El. 

The results of the orbital program are summarized in Table 3.2-1. Plots 
for the angular velocity for September are shown in Figure 3.2-2, plots for 
the angular acceleration for September are shown in Figure 3.2-4. The 
resulting point ahead is shown in Figure 3.2-4. 


TABLE 3.2-1. SUMMARY OF GIMBAL REQUIREMENTS (MAXIMUM VALUES) 

DETERMINED BY THE LEO AND GEO ORBITS. 


Pointing Error 

Angular Velocity 

i 

Angular Acceleration 

Spacecraft 

i 

Az 

El 

Az 

EL 

Az 

El 

LEO 

i 

+ 175* 

+ 140* 

200 urad 

sec 

200 urad 

sec 

1 urad 

0.5 urad 
sec^ 

! 

j 

GEO 

+ 20* 

1 ♦ 

in> 

o 

© 

1 urad 

sec 

200 urad 

sec 

1 urad 

1 urad 
sec^ 
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3.3 INFLUENCE BY THE SUN AND MOON 


The position of the sun and moon relative to the line of sight between 
the LEO and GEO spacecraft are of concern. If either is within a few degrees, 
link outage could occur due to increased background noise on the data channel 
or due to link disruption caused by automatic sensor surge current protection 
mechanisms. The position of the sun and moon relative to the line of sight 
may be calculated with the orbital program. Their orientation for September, 
for example, is shown in Figure 3.3-1 and that for December is shown in Figure 
3.3-2. 

We can use this orbital information to determine the disruption periods 
to the communication link due to the sun or moon in the field of view. We 
show in Figure 3.3-3, for example, the total time that the sun is within 1* of 
the line of sight between the LEO and GEO. The link outage for the sun in a 
4* fleld-of-view of the line of sight is shown in Figure 3.3-4. 

Link disruption due to the sun in the field-of-view is not the only cause 
for link outage. The link will be disrupted when the LEO is eclipsed by the 
earth and also when the gimbal is within gimbal-lock range. We assume the two 
axis (Az, El) gimbal configuration as it is less complex than a 3 axis gimbal 
which would have zero time loss due to gimbal lock. We have tabulated these 
efficts in Tables 3.3-1 and 3.3-2 for each month of the year. In Table 3.3-1, 
we show the percent of the time that the link is available from geometry (not 
earth eclipsed), time loss due to gimbal-lock and time loss for the sun within 
4* of the line of sight. The associated times in time counts (1 count = 
0.9375 min) are noted in Table 3.3-2. In general, link outage due to gimbal 
lock is less than about 33000 seconds per month (< 2.8 percent per month). 
Outage due to the sun within 4* of the line-of-sight is less than 180 sec per 
month (< 0.55 percent of the month). 

Orbital constraints will not preclude operation of the link but will 
cause some outage. Once the actual orbits are known, the outage may be 
predicted. 
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4.0 SYSTEM CONFIGURATION 


4.1 CONCEPT 

The configuration concept for the optical communication terminal on the 
LEO (Landsat) and GEO (TROSS) satellites is needed to complete the task 
requirement. The configuration and design, for example, will Influence the 
way In which the system responds to spacecraft related disturbances. In a 
like manner, the design will Influence the response of the spacecraft to 
moving elements In the terminal. Similarity, size, weight and power estimates 
will be directly Influenced by the design concept. 

An optical communication system has two major subsystems; one electro- 
optical In nature and the other mechanical and electrical In nature. A basic 
block diagram is shown In Figure 4.1-1. The implementation and configuration 
of the system can take several forms. Since the main concern for rutur- 
systems Is accurate acquisition, pointing and tracking, the concepts for the 
system configuration were primarily based on the consideration. 

Two basic concepts were evaluated. One has all optics and electro¬ 
optics, such as the laser source, receiver and associated electronics borne 
on-glmbal. The laser light Is directed to the distant target by moving the 
two axis gimbal which points the electro-optical portion of the terminal. The 
second has only a gimbal-borne flat mirror and the electro-optical components 
are mounted on a nearby baseplate. The laser light is reflected and directed 
by the flat mirror to the distant receiver. Two axis gimbal ling Is required 
as a general feature due to the orbits associated with the LEO (Landsat) and 
GEO (TDRSS) spacecraft. 

The gimbal range for the LEO spacecraft was found by orbital analysis In 
section 2.0 to be ^175* Az and ^140* El. Such a range would preclude a 2 axis 
glmballed flat mirror concept for the LEO spacecraft. The gimbal range for 

the GEO spacecraft was only +20* In Az and El and a glmballed flat mirror 
could be used. 


































































The gimballed flat mirror must be larger than the primary telescope 
diameter (20 inches) by at least a factor of JT) Similarly, the mirror 
support structure and gimbal mechanism must all be larger and more massive. 
As a result, the weight of either approach is the same to first order. We 
selected the two axis gimballed optics approach as our baseline design. In 
this way, the design and some subsystems and components would be Identical in 
the LEO and GEO terminals. Recurring and non-recurring costs would be 
minimized. 

In order to complete the baseline concept design, we evaluated the use of 
common versus separate optics for transmit and receive. Common optics 
requires dichroic elements to separate the transmit and receive wavelengths to 
avoid self jamming. The reguired isolation is of order 90 dB. Our experience 
with dichroic elements Intended for space applications and long life has not 
indicated that a 90 dB isolation can be achieved and maintained. 
Consequently, the baseline configuration is a gimballed electro-optical system 
employing separate transmit and receive optics. An added benefit is that only 
one wavelength is required for both data link and return beacon. 

The baseline configuration is noted in Figure 4.1-2. The block diagram 
and optical schematic of the transmitter are shown in Figure 4.1-3 and 4.1-4, 
respectively, and that of the receiver are shown in Figures 4.1-5 and Figure 
4.1-6, respectively. The use of separate transmit and receive optics allows a 
rather simple optics layout as shown in the optical schematics. The point 
ahead optics in the transmitter and the tracking optics in the receiver are 
similar moveable electro-optical elements. A two axis image motion compensa¬ 
tion (IMC) mirror assembly is used. . The unit has been developed at the Hughes 
Electro-Optical facility. An exploded view of the two axis IMC is shown in 
Figure 4.1-7. 

The basic block diagram for the components of the LEO spacecraft noted in 
Figure 4.1-3 and that for the GEO spacecraft noted in Figure 4.1-5 are similar 
in nature, indicating the gimballed .electro-optical components approach. All 
support electronics are shown off-gimbal. The various electronic units 
required are identified. The servo control electronics are used to determine 
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Figure 4.1-4. Optical Schematic-Data Transmitter. 
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the glmbal pointing and remove baseframe angular jitter from the llne-of- 
slght. The acquisition electronics are used to deduce the angular position of 
the distant terminal within the large fleld-of-vlew of the acquisition 
sensor. Telemetry and command electronics, thermal control and power 
conditioning are Included as required for any space terminal. 

4.2 TERMINAL DESIGN FEATURES 

The evaluation of terminal acquisition are described In Section 8.0. The 
evaluations of terminal tracking are described in Section 7.0. As a result of 
these evaluations, specific features of the terminal design have been Identi¬ 
fied. Other features have been identified after evalution of the requirements 
and orbital constraints. The resulting terminal design parameters are sum¬ 
marized in Table 4.2-1. Th concept design function definitions are summarized 
In Table 4.2-2. 

It is Instructive to note two features of the parameters noted In Table 

4.2- 1. As a result of the pointing error and burst error evaluations in 
Section 5.0, a maximum transmit optical diameter of 6 Inches Is associated 
with a 1-iirad (la) pointing error. A lesser pointing error cannot be easily 
maintained wiht the error baseframe vibrational noise Inherent on the LANDSAT 
spacecraft. Smaller optics may be used for transmission of the acquisition 
and beacon beam from the GEO spacecraft. We use two 6 inch optical antennas 
in the LEO for commonality of design; one for transmit, one fore receive. 

4.3 ACQUISITION CONFIGURATION 

A block diagram of the acquisition/tracking function is noted in Figure 

4.3- 1. Note that an area array has been selected as the acquisition sensor. 
A bifurcating mirror provides the fleld-of-vlew reduction between acquisition 
and tracking. Silicon sensors are assumed since GaAlAs laser wavelengths 
(about 0.8 urn) are used. The component selection for the acquisition and 
tracking functions are identified in Sections 7.0 and 8.0. 
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TABLE 4.2-1. TERMINAL CONCEPT DESIGN PARAMETERS 


Property 

1 GEO 

1 LEO 

Concept 

• Separate Transmit 
and Receive Optics 

• Separate Transmit 
and Receive Optics 

Transmit Optics 

• 

• 2.3 Inch Diameter 

t 6 Inch Diameter 

Receive Optics 

• 20 Inch Diameter 

• 6 Inch Diameter 

Primary Material 

• Beryllium 

• Beryllium 

Acq Sensor 

• SI-Quad 

• 

• Si-Array 

Track Sensor 

• Si-Quad (APD) 

• Si-Quad (APD) 

Comm Sensor 

• SI-APD 

• SI-APD 

Laser Source 

S 

• GaAIAs 

• GaAIAS 

Transmit Wavelength 

• Beacon (LDR) = 0.85 urn 

• Comm (HDR) = 0.83 pm 

Receive Wavelength 

• Conn (HDR) « 0.83 pm 

• Beacon (LDR) = 0.85 

Gimbal 

• Two Axis-Optics on 
Gimbal 

• Two Axis - Optics on 
Gimbal 

Az Range 

• ♦ 20* 

• ♦ 175* 

El Range 

• ♦ 20* 

• ♦ 140* 

Weight (with 
redundancy) 

211.5 lb 

175.2 lb 


Power 


• 143.1 watt 


• 137.8 watt 






























TABLE 4.2-2. CONCEPT DESIGN FUNCTION DEFINITIONS 


Function 

Design 

Parameter 

Acquisition 

• Concept 

• Field-of-Regard 

• Angular Scan Rate 

• Laser Wavelength 
e Laser Power 

• Bean./idth 

• Acquisition Sensor 

• Acquisition Time 

a Initiated by GEO 

• 0.5 deg Cone for GEO and LEO 

• 1.6 deg/sec 

• 0.85 pm 

e 50 mW - cw 

• 130 prad (full Airy Lobe) 

• Nominal 400 x 500 area array 
on LEO 

Quad Cell on GEO 

• £ 2 minutes j 

Tracking 

• Concept 

• Field-of-View 

• T C Wavelength 

• T C Data Rate 

• GEO Beacon Power 

• GEO Beacon Beamwidth 

• GEO and LEO Sensor 

1 

• GEO Tracks on Data Uplink 

• LEO Tracks on T C Downlink 
e 140 prad full Sensor, 

13.7 prad Closed Loop i 

• 0.81 pm > 

• 10 Kbps 

• 20 mW - cw j 

t 35 prad (full Airy Lobe) . 

0 APD - Quad Cell 

_ ! 

Conmunications 

Point Ahead 

• Concept 

• Power 

• Beamwidth 

• Data Rate 

• Wavelength 

• Concept 

i 

0 Multisource Suirming Device ' 

0 200 mW - Average 

0 10 prad (1/E 2 )* 

0 500 MBPS 

00 .8-0.9 pm 

: 

0 Closed Loop Update via T $ C 
Link for LEO Uplink j 

0 Minor Update for GEO Downlink 

* Not consistent with burst error requirements and 1 prad pointing capability 
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Acquisition/Tracking Basic Block Diagram 



































5.0 BURST ERROR 


5.1 INTRODUCTION 

Free space optical communications between satellites can be considered to 
Involve three major system control functions: a) the laser aided acquisition 
of one satellite with the other leading to the fine pointing of the communica¬ 
tion beam to the receiver, b) the maintenance of the cooperative tracking and 
tracking between opposite terminals and c) the continued operation of the data 
transfer channel. These system functions are Interrelated, Involving many of 
the same electro-optical components. 

The data transfer channel and Its electro-optical components have been 
extensively developed and tested In the laboratory over the last decade [1-16]. 
The pointing and tracking features have also received attention, but with not 
as much emphasis [17-23]. The lack of detailed attention, In comparison to 
that given to the communication subsystem. Is apparently the result of the 
opinion that the pointing and tracking function Is limited more by engineering 
than by technology. The pointing and tracking subsystem uses gimbals, servo 
control electronics, optical sensors and beam motion mirror devices. Most of 
these components have been used In other systems, are not necessarily state- 
of-the-art, and their performance Is predictable by appropriate analysis. 
However, although the subsystems have received some ground and aircraft study, 
little is known of the satellite environment. Information regarding satellite 
attitude control and platform stability at the appropriate accuracies and 
angular motion rates (prad and mrad/sec, respectively) have only recently 
become available. 

As such, the pointing and tracking simulations In the terrestrial envi¬ 
ronment can be considered as only order of magnitude correct for any space 
application. We develop the relations which connect the main parameters of 
the pointing and tracking subsystem with those of the communication subsystem. 
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We also Identify the pointing and tracking related constraints on the applica¬ 
tion of an optical communication system to space. The constraints lead to 
design limits for the optical coriinunlcatlon terminal and associated restric¬ 
tions In Its successful Implementation for space missions. 

5.2 SYSTEM CONFIGURATION 

Angular pointing and tracking between two satellite-borne optical com¬ 
munication terminals may be accomplished only by cooperative means as a result 
of the large distance and rather low laser power levels available. Each 
satellite terminal emits a laser signal which Is a tracking beacon for the 
other terminal. Only after an accurate 1 Ine-of-slght (LOS) Is established and 
accurate pointing maintained, may successful data transfer occur. Mlspolntlng 
between terminals leads to system performance degradations and continuous 
accurate pointing Is required for optical communication systems as they employ 
rather narrow beam'widths. 

The accuracy of beam pointing Is dependent upon many parameters, some of 
which are not specifically related to the optical communication system. Such 
parameters are Inherent with the host satellite and Include: a) satellite 
attitude Instabilities, b) satellite ephemerls algorithm errors, and c) 
system baseframe vibrational noise. Accurate beam pointing Is also dependent 
upon system related factors such as: a) the tracking sensor noise, b) signal 
timing errors, c) microprocessor computational errors ?iid d) mechanical 
misalignments. These factors are normally combined In rms fashion to yield 
the system pointing error value. As a result, the accuracy of beam pointing 
will be recognized as a primary system design requirement rather than a 
parameter of choice. 

Although the pointing and tracking control system Is distinct from the 
communication system, the two are Interrelated In defining the performance 
possible with a system operating In the disturbance environment of the host 
satellite. This Interrelationship Is best seen by considering the basic 
system block diagram shown earlier In Figure 4.1-1, and the basic optical 
layouts which were shown schematically In Figures 4.1-4 and 4.1-6 [20]. It is 


5-2 











shown In the figures that the system Incorporates two functional parts: the 
opto-mechanlcal part responsible for transmission and reception of the laser 
light, and the electronic control part responsible for the dynamic control 
(that Is, pointing and tracking) of the optomechanical part. 

In Figure 4.1-4, the light path for tracking of the beacon light from a 
distant terminal Is Illustrated by the dark line and the reverse light path 
for data transmission Is noted by the dark line In Figure 4.1-6. The basic 
optical, mechanical and opto-mechanlcal elements required for the two Inter¬ 
related functions are noted In the figure. 

Two opto-mechanlcal elements are especially Important for accurate 
pointing. The first Is the Image motion mirrors which are used to orient and 

compensate for LOS angular jlter Induced by mechanical and baseframe distur¬ 

bances on the host satellite. The second Is the point ahead mirrors which are 
used to provloe the angular pointing angles required due to relative motion 
between the two satellites. 

As shown in Figure 4.1-4 and Figure 4.1-6, portions of the optical path 
are common to the pointing and tracking functions while other portions are 
dedicated to only one function. Consequently, the errors associated with 
m^stracklng are different from those associated with mlspointlng. The 
pointing errors are greater than the tracking errors since the former are 
dependent upon the latter and must also Incorporate the errors associated with 
additional opto-mechanlcal devices, such as the point ahead elements. 

Point ahead error Is basically a low frequency, almost static offset 

related to the angular bias from the LOS between the two satellites due to the 

range related time delay for light transfer from one terminal to the other. 
The point ahead error Is associated with ephemerls errors, algorithm errors 
and alignment errors. Additionally, the point ahead device must also Incorp¬ 
orate higher frequency compensation unless another such element Is used to 
correct for satellite Induced pointing errors associated with attitude 
Instabilities and baseframe vibrations. 

The manner In which the various noises contribute to the tracking error 
Is Indicated by the simplified servo control block diagram shown In Figure 
5.2-1. The mein noise contributions to the tracking function may be 
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sumnarized as: a) mechanical and thermal base motion disturbances to the 

coarse and fine track loops, b) electronic noise associated with process 
logic and basic electronic circuitry and c) sensor noise which includes dark 
current, excess noise due to gain, thermal noise and optical background 
noise. The sensor noise may be interpreted as a noise equivalent angle (NEA) 
thereby establishing a lower limit to the tracking and pointing error. 

5.3 POINTING ERRORS 

The performance of a free space, satellite-borne optical communication 
system is described by two major parameters: the probability of bit error and 
the probability of burst error. The probability of bit error is normally 
specified as a system requirement. It determines the nominal signal intensity 
required to recover the signal from the communication channel with its contri¬ 
butory noises. The nominal probability of bit error is proportional to the 

signal and noise counts received during the bit time interval [7, 3, 24]. The 
probability of burst error on the other hand, is frequently ignored as it is 
not directly dependent upon the optical signal and noise levels. The probab¬ 
ility of burst error establishes the percent of time that the servo system 
will maintain pointing within some error limit. Should the pointing error e 
become greater than some value e*, the signal intensity at the distant 
receiver would be decreased, and the system performance described by the 

probability of bit error would be degraded. 

There are several parameters which are commonly used to describe the 

nature of mispointing. These include: a) the radial angular mispoint error 
denoted by e and b) the rms standard deviation o of the probability distri¬ 
bution for mispointing, P(e). We will first determine the proper relationship 
between e and the standard deviation c. We will examine the relationship 
between the probability distribution P(e) and the desired communication 
channel performance as measured by the probability of bit error. We will then 
determine the relationship between the radial pointing error and the beamwidth a. 
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Figure 5.2-1. Simplified servo control block diatram indicating the 

major functions and the major noise sources leading to 
mistracking and mispointing. 




































The parameter Is most importance is a/a since the pointing error e is 

normally defined as some value y at the no level (N = 1, 2 or 3). The 

★ 

pointing error probability distribution P in terms of the ration a/a. 
the parameter a/a can also be used to establish the practical limit on the 
diameter D of the transmit optical antenna since a is proportional to 1/D. 

We assume a two axis pointing and tracking control system incorporating a 
quadrant tracking detector and a gimballed optical subsystem. We describe the 
single axis probability distribution for mispoint as a gausslan function with 
single axis pointing errors of e and e [22]. A gausslan distribution is 

X / 

assumed since the major noise source to the pointing error is usually sensor 
noise. The gaussian function has rms standard deviations o ¥ and o . We 

* 9 ^ 91/9 

assume radial svnmetry so that o = o = a and e = (e x c ) 

r a> x y x y 

and require a normalized P (t)\ Q J P(e) de * 1. The radial probability dis¬ 
tribution, in normalized form, may be recognized as the Rayleigh distribution. 



2 

,-e /2a 



The nominal pointing error is given by 
tion, c/a = /»/2 . 


0 



CO 

p(c) e dc and has the solu- 


The probability of burst error is of pr<T;e Importance here and is given by. 
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The solution of Equation 2 yields. 




(-In P *) 
E 



which essentially establishes the probablity of burst error possible given the 
single axis gausslan pointing error distribution Inherent with a system 
operating within the limits Imposed by the host satellite. In a like manner, 
one may specify P^* and given the pointing error probability distribution 
with its characteristics a, the required pointing error must be maintained at 
6 £ e*. 

The ratio e*/a can be related to the beamwidth a, and consequently, the 
system may be designed to operate at the desired performance level for all 
e £ n*. We remember that the bandwidth so defined will be consistent with the 
burst error probability. The Intensity of the received beam power at e = c* 
must be sufficiently large enough to maintain the nominal probability of bit 
error at the desired level. 

We determine the relationship between the pointing error E = e* and the 
basic system optical design parameter (the beamwidth a), by considering the 
optical power collected at the distant satellite receiver, 

I(e,a) LAP 

- -*4^ (4) 

irR a 


where L represents the various optical losses in the transmitter and recei- 
ver (l_ 0 £ 1.0), = D^/4 Is the area of the receiver optical antenna, 

D is Its diameter, R Is the distance between satellites, P„ is the laser 
power, and I(e,a) is the far field beam profile. We use the symbol - 
here for the angular mispoint at its Interpretlon as c or e* is not Important. 
The far field beam profile is the Airy function given as 


r 2J, (2.44* e/a); 

«•••> * 'A-br. ——] 


(5) 


2.44-ir e/a 
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where Jj Is the Bessel function of the first kind. The full Airy beamwidth 
to the first dark ring In the diffraction pattern Is given by 

a = 2.44X/D (6) 

It Is frequently thought that by minimizing the pointing loss relative to 
the beamwidth selected to achieve the desired probability of bit error per¬ 
formance, the communication channel would be optimized. However, Equation 4 
can be solved for the optimum value of the beamwidth relative to the pointing 
error which is inherent In the Implementation of the system on the host 
satellite. It is optimized by setting. 


dP(e,q) 
da 



The solution to Equation 7 relates the channel beamwidth to the mispoint angle 
used in the link budget as, 

t*/a = 0.24 (8) 

The point at which this falls on an Airy far field beam profile is shown 
in Figure 5.3-1. The nominal mispoint loss which should be used in link 
budget summarizes is about -4 dB. 

The primary optical system design parameter was noted earlier as the ratio 
a/a which defines the optical beamwidth needed to maintain system operation at 
the desired probability of bit error. Combining Equations 3 and 8, we find, 

0.17 

-f- = ,- - = (9) 

0 /* 

/<-!" P E I 
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where it will be remembered that c Is the rms standard deviation of the 
pointing error distribution as determined by the dynamics of the host 
satellite and the operation of the optical system. 

The degree of accurately pointing from one satellite-borne terminal to 
the other, as noted above and Illustrated in Figure 3, Is dominated by noise 
error sources associated with the host satellite and with the operation of the 
optical system. As we have placed practical limits on accurate pointing, we 
can establish practical limits on accurate tracking. 

The tracking sensor In an optical conmunication system Is normally a 
quadrant type sensor, and as such, has a limit to Its tracking accuracy a 
given by [25], 

ir 6 

a * -. (10) 

8 yS/N 

5.4 SYSTEM EXAMPLE 

The system pointing error is the rms value of various error sources 

contributed by the host satellite and the optical system. It is not possible 

here to provide an allocation to each noise source as such is application 

dependent. We will specify the pointing error in general terms as a value e = 

★ 

e* relative to the desired burst error P- over a range of representative 
values. 

We Illustrate the significance of the probability of the burst error In 
terms of a/a by plotting the relationship in Figure 5.4-1. We note that 
Equation 9 specifies the optimum power collected at the receiver. We have 
also summarized a number of the relevant example parameters in Table 5.4-1. 
We use a wavelength of x = 0.83 urn consistent with the Al(GaAs) laser diode 
optical source. 

It may at first seem surprising that the ratio of the beamwldth to the 
standard deviation of the radial pointing error distribution is larger than 
about 5 to 1 as was assumed or Implied by earlier evaluations [7,8]. The 
cause of this increase Is that the communication system Is directly influenced 
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by the probability of burst error and was not so considered In the earlier 

* -3 

evaluations. A burst error of say, P r = 10 , would require a beam- 

» g 

width a = 15.4 prad when o * 1 urad, whereas a burst error of P r = 10 

L * 

would require a = 22.1 urad when o > 1 prad. A burst error of, say, P r 

-3 L 

» 10 means that a system Initially configured to have a probability of bit 

error of, for example, P- . 1 x 10 , would have a temporal availability 

of only 0.999 as a result of the burst error. It would have degraded 

performance for 0.1 percent of the time and the probability of bit error would 

be > 1 x 10"® since Pr Is strongly dependent upon signal intensity 

[23,26]. Hence, operation of the system at the optimum beamwldth once o and 
★ 

P- are known Is necessary to maximize performance. 


TABLE 5.4-1. EXAMPLE SYSTEM REQUIREMENTS UNDER 

BURST ERROR LIMITATIONS 


System 

AvallabllIty 
(percent) 

PBE 

e/a 

a/a 

Optimum Antenna 1 
Diameter for 
o s 1 prad and 
x > 0.83 pm 
(Inch) 

Relative 

Power 

Required 

(dB) 

90 

io-i 

2.13 

8.86 

22.8 

0.0 

99 | 

10-2 

3.03 

12.60 

16.1 

3.1 

99.9 

10-3 

3.70 

15.39 

13.2 

4.8 

99.99 

10-4 

4.28 

17.80 

12.4 

6.1 

99.999 

10-5 

4.79 

19.93 

10.2 

7.0 

99.9999 

10-6 

5.25 

I 

21.84 

9.3 

7.8 
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We also note that the correct pointing error loss which should be used In 
link budget and performance tabulations is not about ~0.2 dB ([7,8] but rather 
about -4 dB. This results from the Intensity decrease In the Airy profile as 
given by Equation 5 when the optimum ratio o/a * 0.24 per Equation 8 Is 
Included. This -4 dB pointing error loss can significantly Influence the long 
term performance of a communication system If Ignored and Insufficient power 
margin Is allocated. 

As an Illustration of the influence of mlspolntlng on the optical system 
design, we have combined Equations 6 and 9 to relate o to the telescope dia¬ 
meter D, so that 

oD ~ - =L== — (ID 

2.4 / (- fnP £ *) 


We Illustrate the result In Figure 5.4-2 for x s 0.83 um. We have 
Included curves for several representative values of the burst error 

* - . . ‘ * .5 

P^. We note, for example, that for = 10 and o = 1 urad, the 
optimum transmit telescope diameter to maximize the power collected at the 
receiver Is about 10 cm. The curious but Inescapable conclusion is that once 
the la value of the pointing error distribution has been determined for the 
optical communication system due to the environment on the host satellite and 
Its performance In that environment, and the desired burst error selected, the 
optimum value of the transmit antenna is specified by Equation 11, and larger 
or smaller diameters only lead to degraded performance. 

We can also evaluate the implications to the tracking sensor and Its 
servo control loop by combining Equations 6 and 10 

2x 

aD - - (12) 

/S/N 


5-13 






















We Illustrate the result in Figure 5.4-3 assuming a diffraction limited spot 
. size for x ■ 0.83 pm on a quadrant type sensor. Curves for three selected 
values of the antenna diameter are shown. We note that It Is required that 
a/c < 1.0. As Indicated by the curves In Figure 5.4-3 there Is a lower limit 
to the S/N In the tracking sensor which must be maintained. The so determined 
value of the S/N In the tracking sensor will Impact the design of the tracking 
channel per the normal range relation. 

5.5 SYSTEM IMPLICATIONS 

% • 

The mlspolntlng of an optical conmunlcatlon beam from one host satellite 
to another has been evaluated. We have optimized the optical power collected 
at a distant receiver to show that an optimum ratio of the pointing error to 
beamwldth (c/a) may be deduced. We find that the resulting pointing error 
loss Is, In general, larger than that previously used In many system designs 
and a pointing error loss of about -4 dB must be allocated to the link budget. 

We have modeled the single axis pointing error as having a gausslan dis¬ 
tribution and deduced that the radial pointing error has a Rayleigh distribu¬ 
tion. This radial pointing error distribution was used to define the probabi¬ 
lity of burst error which was shown to be dependent upon the ratio of the 
pointing error and the standard deviation of the pointing error distribution. 

This relationship was combined with the earlier derived optimum ratio of 
c/a to define the design parameter o/a as a function of the pointing error. 
Once a burst error probability Is selected for the system operation and the lo 
value of the radial pointing error distribution of the host satellite known, 
the optimum beamwldth to maximize power collected at the receiver may be 
determined. 

The design ratio a/a defines the appropriate transmit telescope diameter 
to optimize system performance. This result Is In contrast to many who have 
erroneously assumed that a larger transmit antenna diameter and hence smaller 
beamwidth is beneficial to the system. We have also shown that deviations 
from the optimum o/a ratio will lead to an Increase In the optical power 
needed to maintain the probability of burst error at the desired value. 
Additionally, the pointing requirements may be used to establish the design 
limits for the tracking channel. 
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RATIO TRACKING ERROR TO POINTING ERRORS. A/t 


Figure 5.4 


3. A limit to the signal to noise ratio In the tracking 
servo control loop is defined by the ratio of the 
tracking error to the pointing error for several values 
of the receive antenna diameter. 











As an example, we show a link budget summary for a 500 Mpbs data channel 
at 0.8 pm wavelength In Figure 5.5-1. The example assumes a low earth 
orbiting (LEO) satellite transmitting to a receiver on a geosynchronous earth 
orbiting (GEO) satellite. The LEO antenna Is 15 cm (5.9 inches) In diameter. 
The calculations assume a pulse position modulation (PPM) format with m « 4, 
and a maximum a priori (MAP) receiver configuration. A link margin of only 
about 3 dB Is available. 

It has been determined that an optical conmunlcatlon system cannot be 
properly design without knowledge of the pointing error distribution associ¬ 
ated with the host satellite and the operating system. We believe the 
dominant design parameter to be the proballIty of burst error and, hence, the 
c/a ratio. We believe that Its Influence Is far reaching and can control the 
eventual application of otlcal communication systems for space mlslons. 

















TABLE 5.5-1. COMMUNICATION LINK BUDGET SUMMARY 

(500 Mbps, PPM, m * 4, MAP DETECTION) 


• 

ITEM 

FACTOR 

dB 

Peak Laser Power 

0.8 watt 

-0.97 

LEO Optics/Combining 

0.45 

-3.47 

LEO Antenna Gain 

15 cm diameter 

115.3 

Pointing Loss 

0.4 

-3.98 

Range 

44,500 km 

-296.36 

GEO Antenna Gain 

50 cm diameter 

125.33 

GEO Optics. 

0.32 

-4.95 

Received Peak Power 

i 

-69.16 

Probability Bit Error 

1 x 10- 6 

- 

Probability Burst Error 

1 x 10- 3 

- 

Required Peak Power 

- 

-72.36 

Margin 

’ 

3.2 dB 
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6.0 SERVO ANALYSIS 

6.1 Introduction 

# 

This section addresses all the servo requirements that were either well 
defined at the start of the study or evolved as the study progressed. It also 
defines a system baseline design and mentions those significant assumptions that 
qualify the results. 

6.1.1 Servo Design Requirements 

The primary purpose of this study was to establish a servo design confi¬ 
guration that would maintain a prescribed tracking accuracy between two earth 
orbiting satellites. Each satellite's servo system must track the other sate¬ 
llite's laser communication beam with an accuracy (2-axis) of 1 urad. 

Of seconcary concern was to design a servo that minimizes the reaction 
torque imparted to the spacecraft. This was accomplished by minimizing the 2- 
axis giiffbal dynamics. The RMS reaction torques for the various servo designs 
were compared for base motion and friction disturbances. 

Finally, a third requirement that evolved during the study emphasized 
design simplicity and reliability. The requirement for all components and devices 
to be space qualified, led toward a simplistic design which minimizes costs and 
boosts reliability while still meeting system performance goals. 

6.1.2 System Baseline Design and Assumptions 

The baseline system configuration is the same for both spacecraft trackers. 

It consists of two Cassagrain telescopes - one receiver and one transmitter - 
each mounted to the inner gimbal of a two gimbal set (outer azimuth and inner 
elevation). The optical train of each telescope contains a small (1 1/2 inch 
diameter) 2-axis IMC (image motion compensation) mirror which is capable of operating 
at high servo bandwidths ( hundreds of hertz). The receiver's optical train also 
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contains a track sensor that generates a 1ine-of-sight (LOS) error. This error 
is nulled by simultaneously torquing the azimuth and elevation gimbals, and the 
2-axis IMC mirror. 

Since both telescopes are mounted on the same gimbal set, the transmitter's 
LOS pointing error is very similar to the receiver's LOS tracking error. There¬ 
fore, an error analysis was performed on a receiver servo and then modified to 
reflect the transmitter's error sources (point ahead). Figure 6.1-1 shows a 
block diagram relating the receiver and transmitter functions. 



Figure 6.1-1: Transmitter/Receiver Functional Block Diagram 
The schematic diagram for the baseline on-gimbal receiver configuration is 
presented in Figure 6.1-2. This diagram was used to generate the angular rela- 
tionships of the optical LOS (output space) with respect to gimbal and IMC mirror 
positions (input space). 

*This filter functionally represents the frequency splitting of the track error 
that occurs in each of the candidate serve configurations in Section 6.2 
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Figure 6.1-2: Receiver Telescope Schematic Diagram 

The track sensor and telescope are rigidly mounted to the elevation fix- 

* 

ture. Any gimbal motion directly transmits into track sensor and telescope 
motion. The IMC mirror, however, is shown to be soft-mounted to the fixture to 
minimize gimbal motion coupling to the mirror. A mathematical justification for 
soft-mounting the IMC mirror is discussed in Appendix 6A.1. 

Both elevation and azimuth axes are represented by a single-axis analysis. 
All error sources were assumed identical for both axes except for base motion. 
Spacecraft base motion coupling into the optical LOS is discussed in Section 
6.3.2. To simplify the analysis, azimuth and elevation angles were assumed to be 
small sucn that their cosine functions equal one and their sine functions equal 
the angle. 

In performing the tracking analysis, various assumptions were made to 
bound the scope of this study and to account for the limited amount of satellite 
data available to us. NASA/GSFC provided HAC with base motion disturbance data 
measured on the LANOSAT spacecraft. To maintain a conservative approach, a worst 
case envelope was selected from this data to represent angular base motion distur¬ 
bances about each of the three spacecraft axes - roll, pitch, and yaw. 






















Linear disturbances* were not considered in this analysis because of their 
negligible effect on tracking accuracies. Linear accelerations can be transformed 
to angular accelerations if the telescope centers-of-gravity (C.G.'s) are not 
coincident with the gimbal axes. This phenomenon! called "mass-unbalance" is 
minimized by precisely balancing the on-gimbal weights about each of the gimbal 
axes. 

6.2 Servo System Configurations 

Three servo configurations were considered for this study - gyro-stabilized, 
mass-stabilized, and complementary filter. All three configurations control the 
system's inertial LOS angle such that the received laser beam is maintained near 
the tracker's center null position. 

The gyro-stabilized configuration is a conventional design that has been 
implemented on many HAC built tracking systems (for ground, air, and space appli¬ 
cations). It has an outer low bandwidth track loop with two higher bandwidth 
inner loops for attentuating base motion disturbances. It is a proven candidate 
design that can also be used as a baseline to compare to the two other configurations. 

The mass-stabilized design was considered as a simplification to the gyro- 
stabilized design. It replaces the two rate-integrating gyros (RIG's) with a 
reference platform that points along the desired LOS. By replacing the gyros, 
this design improves reliability and reduces system weight and power requirements 
but, as pointed out in Section 6.3, it compromises tracking accuracy. 

Finally, the complementary filter configuration was considered for its 
utmost simplicity of design. It differs from the other two by the fact it only 
closes a single wide bandwidth loop about the track sensor with no requirements 
for inner loops. Although its simplicity reduces the device count, its wide 
bandwidth allows more sensor noise to contaminate the track LOS. 


*The angular accelerations induced by the relative satellite orbital dynamics are 
addressed in Section 6.3.4 with the dynamic equations derived in Appendix 6A.4. 
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The following sections discuss in detail each servo configuration and the 
philosophies behind their designs. A two-axis servo block diagram is presented 
for each design in Figures 6.2-3, 6.2-4, and 6.2-5. It illustrates the various 
servo loop dynamics and their interrelationships and is extremely helpful when 
discussing base motion coupling. Appendix 6A.2 derives the various servo control 
equations. 

6.2.1 Servo Configurations and Their Design Philosophies 

The gyro-stabilized configuration uses three Type II servo loops to control 
the azimuth and elevation LOS axes (see Figure 6.2.1a Single-Axis Functional 
Block Diagrams). The outermost track loop is low bandwidth and obtains its 
error signal from a centroid tracker. This LOS error is compensated by a lead- 
lag network and one electronic integrator as derived in Appendix 6A.2. The low 
bandwidth track loop minimizes sensor noise transmission to the LOS and is adequate 
for tracking targets having small angular accelerations (see Sections 6.3.1 and 
6.3.4 for discussions of sensor noise and track following error). 

The track loop compensates the track error to generate a rate command for 
the inner coarse gimbal (CG) rate loop. This rate command is manifested as a 
torque applied to the gyro's output axis. The angular displacement between the 
output axis and the gyro case establishes the inertial LOS error that is used by 
the.rate loop to torque the azimuth and elevation gimbals. The gimbal motion 
about the gyro’s input axis produces a torque about the output axis which nulls 
the position error. When this occurs, the gimbal (telescope) LOS is coincident 
with the desired track LOS. 

The CG rate loop bandwidth has upper and lower bounds due to stability 
considerations. Because it is an inner loop, its bandwidth must exceed the outer 
track loop bandwidth by about a factor of 10 to achieve a stable track.* Otherwise 

*The CG loop bandwidth can probably be less than 10 times the track loop bandwidth 
because the IMC loop decouples it from the track loop. 























































the phase lag at the track loop crossover frequency will not satisfy minimum 

* 

phase margin requirements ( 30 deg) for stability. 

On the upper bound, the CG loop bandwidth must be less than five times the 
lowest structural resonant frequency of the gimbal. It must also be less than 
five times the RIG bandwidth. Both constraints satisfy CG loop stability require¬ 
ments. Because of these upper constraints, a third inner loop is required to 
attenuate high frequency base motion disturbances. 

The IMC loop is a very high bandwidth position loop ( few hundred hertz) 
that is needed to attenuate high frequency disturbances. If a torque disturbance 
is applied to the gimbal base, the CG loop will remove that portion of the error 
that falls within its bandwidth. The remaining error appears in the gyro error 
signal that commands the IMC loop. This remaining error is the angular difference 
between the desired LOS and the gimbal's LOS. Because it is a correction to the 
gimbal LOS, the IMC loop error is formed by comparing the RIG error output to the 
relative angle between the IMC mirror and the coarse gimbal as measured by the 
LVOT (linear variable differentiating transducer). The IMC loop attenuates that 
portion of the error that does not exceed its bandwidth. 

During the final presentation, there was some doubt expressed by NASA in 
building a 400 hertz IMC loop. An example of a wide bandwidth steering mirror 
built and tested by HAC is the AABISM used in the Air Force's Airborne Pointer 
Tracker (APT) alignment system. 

The AABISM has a 21 cm diameter and an inertia of 0.2 in-lb-sec^ and is 
water-ccoled to align a high energy laser beam. The mirror is flexure-mounted 
and driven by voice coil actuators with a reactionless mass used to minimize 
• torque coupling to the base. 

The IMC mirror proposed by HAC for this study is 1 1/2 in. in diameter and 
has an inertia of 6.4 X 10*3 in-oz-sec^. The mirror is not water-cooled but is 
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flexure-mounted and driven by voice coil actuators. The necessity of using a 
reactionless mass to isolate the base (telescope) from torque disturbances is 
beyond the scope of this study. If required, the only limitation would be the 
volume constraints on the mirror. 

Figure 6-2-2 shows a photograph of the AABISM steering mechanism along 
with the laboratory-measured open loop frequency response of the steering mirror 
loop. The crossover frequency (0 OB) in the plot occurs at 500 hertz with 
structural resonances measured at a few thousand hertz. 



Figure 6.2-2: AABISM Steering Mechanism and Open Loop Response 


The CG and IMC loops effectively split the LOS error into low and high 
frequency components, respectively. This feature is quite desirable when mini¬ 
mizing the reaction torque imparted to the spacecraft by the gimbal motion. 
Because the gimbal is considerably more massive than the IMC mirror, and the IMC 
mirror can use a reactionless drive, reducing the CG loop bandwidth effectively 
reduces the spacecraft reaction torque. However, the smallest attainable CG loop 
band*idth is limited by the constraints mentioned above in addition to the travel 
limits of the IMC mirror (the less error removed by the coarse gimbal increases 
the amount of error to be removed by the IMC mirror). 
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Figure 6.2-lb shows the mass-stabilized functional block diagram. It is 
quite similar to the gyro-stabilized configuration since It closes an outer low 
bandwidth track loop and two inner high bandwidth position loops. The reference 
platform replaces the two RIG's as the LOS reference. The track loop now provides 
torque commands to voice coil actuators which reposition the soft-mounted reference 
platform to point along the desired LOS. The platform provides Inertial dynamics 
to the track loop such that it eliminates the need for an electronic integrator 
that is required in the gyro-stabilized approach. The platform's relative 
angular position to the gimbal is measured by two LVDT's. This slave error 
drives the coarse gimbal position until it is aligned with the reference platform 
(and the slave error is nulled). The slave error output is used to command the 
IMC loop in the same manner as the gyro output was used in the gyro-stabilized 
configuration. 

A shortcoming to the mass-stabilized approach is that any disturbance 
entering the inner CG loop is not totally isolated from the track loop. Gimbal 
base motion couples into the track loop through the single-post flexure at fre¬ 
quencies that are within the CG loop bandwidth. (It can he shown that the mass- 
stabilized LOS response to base motion approaches the gyro-stabilized response as 
the flexure spring rate approaches zero). 

The complementary filter configuration is a simple servo design that closes 
a single wide bandwidth track loop with a minimal number cf active devices (see 
Section 6.4.3 for component list). This design directly splits the track error 
signal into its low and high frequency components via a complementary filter (see 
Appendix 6A.2 for complementary filter derivation). The low frequency error 
component is removed-by torquing the azimuth and elevation gimbals, and the high 
frequency component is removed by actuating the IMC mirror. 
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This approach has a number of advantages over the other configurations: 

1) The only limit on splitting the error into high and low frequencies is 
the travel limit on the IMC mirror. Lowering the low pass filter 
(Figure 6.2-1C) corner frequency Increases the error to be removed by 
the IMC mirror. 

2) There are no inner loops or bandlimlted devices other than the track 
sensor that must be considered for track loop stability. 

3) Changing the error signal's frequency split only requires resetting 
the low pass filter corner frequency to the desired value. The multi¬ 
loop approach requires the CG loop gain and compensation to be redesigned 
and the overall track loop stability reverified. 

A disadvantage to the complementary filter design is that it does not 
filter as much sensor noise from the loop as does a low bandwidth design. For 
this reason, a key tradeoff in this analysis centers around the sensor's noise- 
equivalent-angle, o N ££ (derived in Section 7.2), and its effect on the optical 
LOS (Section 6.3.1). 

6.2.2 Servo Block Diagrams 

This section defines the variables that appear in the servo block diagrams 
in Figures 6.2-3 through 6.2*5. The loop bandwidth values were selected in 
Section 6.3 and were used to generate the error analysis results. The gain and 
compensation equations for each servo loop are derived in Appendix 6A.2 as a 
function of loop bandwidth. 

The proposed IMC mirror for this study is being built by HAC for the HNVS 
program. Its flexure damping ratio and spring rate characteristics were taken 
from actual test data to be 

C m * damping ratio » .01 

co m s spring rate = 2n (1.25 hertz) 
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For the azimuth channels in each block diagram, the coarse gimbal gain is 

adjusted by the term 1 This term transforms the CGerror from LOS to gimbal 

COSCp 

coordinates so that the gimbal torque commands are in the proper coordinate 
system. The 2-axis gimbal piogram in Figure -6.3-4 shows this relationship as 

u> d = gj 0 cos e o (6.2.2.1) 

In the same respect, the base motion disturbances that enter the loops through 
friction and kinematic coupling were derived from the respective 2-axis and 4- 
axis gimbal piograms. 

6.3 Error Analysis 

During the study, various error sources were investigated as to their 
effect on system performance. As a result, three sources dominated - base motion, 
friction, and sensor noise - while the others were considered negligible or 
assumed to be removable through calibration. The following sections examine base 
motion, friction, sensor noise, and other less significant sources to determine 
their contributions to the LOS error for our three candidate, servo designs. 

The error analyses also determined the required servo loop bandwidths that 
would satisfy the system error budget. The sensor noise and following error 
analyses determined an acceptable track loop bandwidth for the multi-loop servo 
designs and the base motion analysis determined the desired IMC loop bandwidth. 
6.3.1 Sensor Noise 

The amount of sensor noise the track loop responds to is a function of the 
track loop bandwidth. The track loop filters those frequency components of noise 
that exceed its bandwidth. This justifies why low bandwidth track loops are more 
desirable in systems where sensor noise is dominant. To determine appropriate 
track loop bandwidths for this study, an analytical model was used which deter¬ 
mines the LOS RMS error output of the track loop due to sensor noise-equivalent- 
ang!e, C N£A . 
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The expression for the LOS RMS error due to sensor noise is derived in 
Reference 3. It assumes uncorrelated sensor error signals at the sampling outputs 
by bandlimiting the sensor noise to half the sampling frequency. The servo 
block diagram representing this analysis appears below. 
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Figure 6.3-1: Track Loop Model for Sensor Noise Analysis 
Sensor noise is represented by a white noise process with amplitude ° N £a 
that is passed through a low pass filter with corner frequency fs/2 (half the 
sampling rate). The filtered noise enters the track loop (represented by H(S)) 
which samples the error signal at f s hertz. The error is held constant for 1/f s 
seconds with a zero-order-hold function. The loop gain and compensation are 
designed to achieve a desired bandwidth of f c hertz with adequate phase margin 
for loop stability (see Appendix 6A.2 for derivation of servo control equations) 
To solve for the RMS LOS error, the wt1 i te hoise PSD is multiplied 

by the magnitude of H(S) squared and integraved over all frequencies. 

fe 


0 


RMS 


[ 1 Js N (f) | H(2irjf) H(-2irjf) | d 


1/2 


(6.3.1.1 


where c / f \ .2 - s < e < f s 

n (f) 0 NEA’ ~1 ~ f " ~l 

Substituting the c'osed loop transfer function for H(S) into 6.3.1.1 yields, 

1/2 
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(6.3.1 .2) 
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where, f c = track loop bandwidth (hertz) 
f s = sampling rate (hertz) 

C ■ track loop damping ratio 

°NEA = sensor noise-equivalent-angle (urad) 

Equation 6.3.1.2 shows to be a function of the ratio of track loop 
bandwidth to sampling rate, f c /f s . For this report, it was decided to fix this 
ratio at 1/10* and determine the effect on as the sensor sampling rate varies 
w’th track loop bandwidth. 


NEA 

is derived 

in Section 7.2 to be 



= it 1.22 X 

1 

(6.3.1.3) 

NEA 

8 0 

5NR 



where X = laser beam wavelength (urn) 

D = aperature diameter (in) 

SNR = signal-to-noise ratio 

If we assume that we maintain a fixed average laser power for all choices of 
sampling rate, then 

1 

SNR a 1/2 (6.3.1.4) 

^s 

But since f r = 1 f s , then 

10 

SNR a —J_ 

1/2 (6.3.1.5) 

f C 

Because of the proportionality relationship in b.3.1.5, a value for SNR 
must be assumed for a given value of f c . It was assumed that the sensor SNR 
equals 5 for a track loop bandwidth of 400 hertz and a sample rate of 4 kilo- 
hertz. Witn these initial conditions and 


*In Section b.2.1, it was stated that for track loop stability, the track loop 
sampling rate must be about 10 times the bandwidth. 
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X - 0.9 urn 


and £ = .725 

• * 

Equation 6.3.1.2 reduces to 



( 6 . 3 . 1 . 6 ) 


This equation was applied to the three configurations for the 6 inch and 
20 inch diameter aperatures yielding the following results. 



Figure 6.3-2: Sensor Noise Results 


6.3.2 Spacecraft Base Motion Disturbances 

As previously stated, LANDSAT base motion disturbance data was provided to 
HAC by NASA/GSFC. On board disturbances were measured by rate integrating gyros 
and angular displacement sensors. Together, they measured spacecraft roll, 
pitch, and yaw angular displacements over a bandwidth of 0 to 125 hertz. This 
data was converted by HAC to power spectral densities (PSD's) by squaring the 
angular amplitudes. 

A worst case envelope of all the data was selected to represent the angular 
disturbances about each spacecraft axis - roll, pitch, and yaw - which appears in 
Figure 6.3-3. 
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Figure 6.3-3: Spacecraft Base Motion Disturbance PSD 

Figure 6.3-3 shows the measured data abruptly terminating at 125 hertz 
whereas the true data may extend beyond this frequency. However, the plot 
suggests that most of the power resides at the lower frequencies since the PMS 
error at 1 hertz equals 86% of the total RMS error at 125 hertz. Because of the 
uncertainty of disturbances beyond 125 hertz, the IMC loop bandwidth was chosen 
to exceed the highest base motion frequencies. 

The base motion disturbance enters the tracking system at the base of the 
2-axis gimbal. These disturbances take the form of roll, pitch, and yaw rates 
and are related to the gimbal's LOS rates by means of a piogram which appears in 
Figure 6.3-4. 

The spacecraft base motion couples into the gimbal LOS through friction 
and kinematic coupling. (Section 6.3.3 discusses friction modelling in detail. 

We will only discuss nere how the base motion influences friction torque.) 

Coulomb friction torque is a function of the relative azimuth and elevation 

gimbal rates between the gimbal base and its bearings. These rates are repre- 

• • 

sented by n Q and Cq, respectively and when integrated, are the azimuth and 
elevation gimbal angles, n 0 and Cq in the piogram. For example, the azimuth 
gimbal base is rotated at the rate, due to spacecraft yaw angle motion. This 
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motion induces a f-iction torque disturbance into the inner CG loop (see Figure 
6.3-4.) Jhe servo loop responds to the disturbance by torquing the azimuth and 
elevation gimbals yielding inertial gimbal rates of u 01 ^ andu) 0 , respectively. 


SPACECRAFT LOS 
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Figure 6.3-4: Two Axis Gimbal Rate Piogram 
The gimbal angle rates are computed as folio./s: 

• r 

\ - w o n - 

and e 0 = “oe - “b 

where u b = ^cos^ - “isinbo 

Base motion that couples :to the LOS directly, without entering through 
bearings or flexure pivots is called kinematic coupling. It is very easy to 
determine the kinematic coupling terms from the rate piogram - it is the only 
path that is not broken for a gimbal bearing. For our 2-axis gimbal, this path 
orginates from affecting the LOS and azimuth rates as follows: 

“d = w ASin G 0 
“ r = U A C0S e o 

or “(j = f°iCos n 0 + U jSin n 0 ) sin € 0 
and “ r s ( w iCOS n 0 + “ jSin n 0 )cos $ 


6-18 













For n Q = u° and e 0 * 90°, the azimuth rate,w d , equals the spacecraft roll 
rate, o>j, and the roll LOS rate equals zero. Since we are tracking a symmetric 

beam, LOS rate is not servo controlled. 

• • 

The gyro-stabilized and complementary filter configurations are two gimbal 
systems. Kinematic coupling is shown to enter the CG loop of the azimuth channel 
in Figures 6.2-3 and 6.2-5. 

The mass-stabilized configuration is a four gimbal system where the two 
additional angles, and e-j, are the reference platform deflections from the 
coarse gimbal. Its rate piogram is more complicated than the two gimbal system, 
but the same rules and characteristics apply. Figure 6.3-5 reveals four broken 

I 

lines - two for the azimuth and elevation gimbal bearings and two for the platform 
deflections. 

SPACECRAFT GIMBAL LOS PLATFORM LOS 



Figure 6.3-5: Four Axis Gimbal Rate Piogram 
It is interesting to note that the kinematic coupling for the four gimbal 
system influences the elevation and azimuth channels. The elevation term is 
quite small because it is multiplied by the sine of 'H -j which is maintained 
near zero (slave error) by the CG loop. The azimuth term affects the CG loop 
in the same manner as the two gimbal configurations. 
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The flexure pivots are represented by the break points In Figure 6.3-5 

• • 

that generate the angle rates n-j and £■(. The flexure pivots have negligible 
friction torque loads but do transmit attenuated base motion disturbances to the 
platform. These disturbances are sensed by the LVOT and enter the track loop 
through the LVOT error signal (see Figure 6.2*4). This coupling transmits distur- 
bances to the platform that exceed the track loop's bandwidth and degrades LOS 
pointing accuracy. 

A comparison of the servo configuration performances appears in Figure 
6.3-6. The left-hand graphs are the Bode plots of the error rejection transfer 
functions* and the right-hand plots are the LOS error residuals due to base 
motion disturbance inputs. 

The mass-stabi1ized error rejection curve shows less base motion rejection 
then the gyrostabi1ized and complementary filter curve for frequencies less 
than 10 hertz. The IMC loops in each configuration equally attenuate any errors 
with frequencies above 10 hertz. This accounts for the identical error rejection 
curves above the 10 hertz frequency. For the mass-stabilized approach, any error 
less than 10 hertz couples into the track loop through the flexure and corrupts 
the reference platform's desired LOS command to the CG loop. The CG loop follows 
this command thereby causing a LOS error to the system. 

The right-hand curves are plots of the RMS LOS error residuals verses fre¬ 
quency. They were generated by a frequency response computer program which 
applies the base motion PSD (Figure 6.3-3) to the azimuth loops in Figures 6.2-3, 


*Error rejection transfer functions show the sensitivity of the LOS to base motion 
disturbances in the frequency domain. 
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Figure 6.3-6: Servo Configuration Comparisons of Base Motion Error Rejection 

through 6.2-5. This generates an output PSD at the LOS which is integrated over 
frequency to obtain the mean-squared error. 

The three configurations show very small differences between their base 
motion LOS residual errors. At frequencies below 10 hertz, the mass-stabilized 
RMS LOS error is about .10 urad larger than the other two configurations, but at 
125 hertz, it differs by only .013 tc .025 urad. This is attributed to the 
improved mass-stabilized error rejection at the higher frequencies. Although 
this improvement is small, it acts over a much larger bandwidth (10-125 hertz) 
than the inferior low frequency error rejection ((J-l0 hertz). 

The desired IMC loop bandwidth was determined by applying the base motion 
disturbances to the gyro-stabilized servo for various bandwidth values and recording 
the RMS LOS errors. This was done for 100 to 400 hertz bandwidths with the 
results plotted in Fioure 6.3-7. 
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Figure 6.3-7 LOS Error vs. IMC Loop Bandwidth 
With a system error budget of 1 urad, the 400 hertz bandwidth was chosen. 
Its LOS error ot 0.5 urad allows enough margin for other error contributions. 
This bandwidth was also chosen to exceed the base motion bandwidth of 125 hertz 
to allow for unforeseen high frequency components. 

6.3.3 Friction Torque Disturbance 

Due to the limited resources and time for this study,,the friction dis- 

* 

turbance was modeled as a simple Coulomb friction torque acting at the gimbal 
bearings. Stiction and motor hysterysis were not addressed. 

Coulomb friction was modeled as a random process from which a probability 
density function (pdf) was derived (See Appendix 6.5.3). From the pdf, a PSD was 
generated to provide us with a frequency domain representation. The PSD from 
Appendix 6.5.3 is factored to yield, 

2AVG- 




/ cnvu * ( 2 a/it 

v j w + 2u' -jw + 2u 


) 


6.3.3.1) 


The factored form provides us with a friction acceleration transfer function that 
was applied to each of the servo configurations in the manner depicted in Figure 


6.3-8. 
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Figure 6.3-8 : Friction Torque Disturbance Input to Servo Loops 


A white noise process was input to the friction transfer function. Its 
output is a random acceleration disturbance that enters the servo loop and causes 
a residual LOS error. The RMS of the LOS error is generated by a computer program 


that performs the following linear mathematics. 


RMS 


« [ 1 f Sv (u>) | H'(jw)H'(-j«J) | du> ] 

/o r 


1/2 


( 6 . 3 . 3 . 2 ) 


H'(S) is the servo loop transfer function that relates the LOS angle.^Los* t0 

• • 

the input friction acceleration, for each of the candidate servo configura¬ 


tions. For the simple block diagram in Figure 6.3-8, 


H'(S) 


LOS 


( 6 . 3 . 3 . 3 ) 


S z + H(S) 


Substituting equations (6.°.3.3) and (6.3.3.1) into (6.3.3.2) yields, 


RMS 


ao 

C f ^^u 

'0 a)2 + 4 u ! 


+ H(jw) + H(-jw) 


d<u ] 


1/2 


( 6 . 3 . 3 . 4 ) 


Simplifying this equation yields a function of u, the average number of 
zero crossings per second. RMS LOS errors for A equal to unity were generated 
for each servo configuration and for varying values of u. That data is plotted 
in Figure 6.3-9 (note that two scales were needed to plot all three curves on 


the same graph). 












♦ 



(NO. OF ZERO CROSSINGS PER SECI 

Figure 6.3-9: Friction Acceleration Effects to Servo Confiouration LOS Errors 
(See Appendix 6.5.3for an interpretation of the Coulomb friction modeling and LOS 
results in Figure 6.3-9)- 

The data at u * 500 is highlighted because that is the estimated number of 


zero crossings per second that is computed from the base motion PSD (Figure 
6.3-3). The algorithm estimating u is derived in Reference 1 to equal 


ao 

/ u 2 S(u>)df 1/2 


* J 


S (at) df 


( 6 . 3 . 3 . 5 ) 


where S (w) equals the base motion angular rate PSD. 
v 

In comparing the three curves, it is apparent that the mass-stabilized 
configuration performs much more poorly to friction disturbances than the other 
two configurations. This is true for two reasons. 

1) In Section 6.3.2, it was pointed out that the mass-stabilized confi¬ 
guration does not reject low frequency disturbances as well as tne 
gyro-stabi11 zed configuration. This fact by itself is not sufficient 
to cause the extreme performance differences. 
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2) The shape of the disturbance input in the frequency domain is also an 
important factor. Figure 6.3-10 shows an asymptotic frequency plot 
of the friction transfer function. 


GAIN (06) 
20 LOG |w( JO))| 


V FREQUENCY 
7T 

Fioure 6.3-10 i Friction Disturbance Angular Acceleration Transfer Function 
The gain in Figure 6.3-10 represents a random angular'acceleration. To compare 
this disturbance to the base motion disturbance and apply it to the error rejec¬ 
tion transfer functions (Figure 6.3-6), it must be converted to an angular 
displacement through double integration. 

Therefore, ]_ W(S) = 1 2Aju_ (6.3.3.6) 

S 2 S 2 S 2 +2u 

and when plotted verses frequency, becomes 




Figure 6.3-11 : Friction Disturbance Angular Displacement Transfer Function 

It is quite obvious that this transfer function predominates at the lower 

frequencies and decreases very rapidly with increasing frequency (1_, 1 ). 

S 2 S 3 
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Therefore, this function shape, together with the inferior mass-stabilized error 
rejection performance, justifies the relative LOS errors between the mass- and 
gyro stabilized configurations. 

The complementary filter error rejection performance is slightly inferior 
to the gyro-stabilized results. This again is caused by two factors. 

1) It was pointed out in Section 6.3.2 that the complementary filter 
error rejection curve at low frequencies does not perform as well as 
the gyro-stabilized configuration but performs better at higher 
frequencies. 

2) The friction disturbance at u * 500 (corner frequency equals 160 hertz) 
has its predominant energy at low frequencies (Figure 6.5-9) . 

Since the friction error is essentially at low frequencies where the comple¬ 
mentary filter performs less favorably, its LOS error is larger. 

The curves in Figure 6.3-9 are normalized LOS errors for a frictior-to- 
inertia ratio of A equal to 1 rad/sec^. For our gimbals, the azimuth and elevation 
values of A were estimated to equal .12 rad/sec^ and .16 rad/sec^, respectively. 

At u equal to 500 zero crossings per second, the LOS errors for each configuration 
are tabulated below. 


CONK I GUMATION 

FRICTION LOS ERRORS 

( URAO) 

TOTAL (upad) 

AZ 

EL 

GYRO-STAB. 

.022 

• Uil 

cc 

*"> 

o 

• 

MASS-STAB. 

b • 79 

9-Hi 

11.64 

COMP. FILTER 

• UZJ 

• 0Z9 

.056 


Figure 6.3-12: Friction LOS Error Summary 
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6.3.4 Other Error Sources 

This section discusses the less significant error sources that complete 
the error analysis. The errors consist of 

1) point ahead 

2) dynamic following 

3) LVDT non-linearity 

4) gyro noise 

6.3.4.1 Point Ahead Error 

The point ahead error is the angular error in pointing the transmitted 
beam caused by the time delay it takes the beam to travel from one satellite to 


the other. The two way transmission time is computed as 
T = 2R 0 


( 6 . 3 . 4 . 1 . 1 ) 


where Rq = relative LOS range between satellites 
C = speed of light 

The point ahead angle equals the change in the LOS angle during the trans¬ 


mission time, T. It is computed from the LOS tangential velocity to be 


0 


which yields. 


V D 


( 6 . 3 . 4 . 1 . 2 ) 


2 V 


0 


D 


(6.3.4.1.3) 


Note that the point ahead angle is only a function of Vq. Therefore, the point 
ahead error is derived from the uncertainty in estimating Vq. 

From previous studies, the tangential velocity error, AVq, was assumed to 
be less than 1 m/sec for orbits where Vq equals 50 km/sec (for our LE0/6E0 orbits, 
maximum Vq equals 12 km/sec). Using a conservative aVq of 10 m/sec, 

9 

= £i!° = .066 urad I 6 - 3 - 4 -' - 4 ) 

0 c 3X10 3 
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for both axes which applies only to the transmitted beam's pointing error. 

6.3.4.2 Dynamic Following Erro r 

The dynamic following error is a DC error that is caused by tracking an 
accelerating target with a Type II track loop. Type II implies that there are 
two integrators in the track loop (either electronic as in the gyro-stabilized or 
kinematic as in the mass-stabilized). 

The following error is the amount of angular position error that must be 
output by the track sensor before the track loop responds to an acceleration 
disturbance. 

For example, in Figure 6.3-13 , a Type II track loop is driving a 
reference platform that is being perturbed by the friction acceleration disturbance, 




Figure 6.3-13 : Track Loop Block Diagram for Following Error 

• • 

The reference platform will not move until the motor's acceleration, @ m , 

• • 

equals or exceeds 9 p. From the figure, 

• • 

9 m = K T G T 6 e (6.3.4.2.1) 


Under steady state conditions, yielding 

8 s k t 0 
m N e 

• • 

To null the disturbance acceleration 9 p, 

_, the dynamic following error. 

<T" 


( 6 . 3 . 4 . 2 . 2 ) 

9 e must build up to or exceed 


• • 
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Formally, the transfer function relating the track error to 9p equals 

1 ( 6 . 3 . 4 . 2 . 3 ) 


8 e (s) 


7 


o (s) S fc +K G 


Under steady state conditions. 


1 im 
s*0 


e a (s) 


_ 1 


• • 


0 F (S) 


O 

where K T = (2"fc) as derived in Appendix 6A.2 

„l/2 

due to an accelerating target is 

.. J/2 


(6.3.4.2.4) 


Therefore, the steady state error 


6 


e f 


( 6 . 3 . 4 . 2 . 5 ) 


ss 


(2ir)‘f 


27 2 


From Appendix 6A.4, it was determined that the maximum total angular LOS 

acceleration between satellites equals U.2 urad/sec^. Setting the lead ratio R 

equal to 6, equation 6.3.4.2.5 is evaluated for track loop bardwidths of 1 and 

400 hertz to yield 

8 e (1) = .01 urad 
ss 


and ® e (400) = 0 urad 
ss 

The dynamic following error can mostly be eliminated by feedforwarding 
estimates of the LOS angular velocity and acceleration data* to the inner loops 
of the mass- and gyro-stabilized configurations. The feedforward signals 
effectively increase the track loop bandwidth by providing anticipation to the 
system. 

6.3.4.3 LVDT Non-Linearity 

An accurate approach to measuring small angular displacements is to use an 
LVDT. The LVOT measures linear displacements that are adjusted by the lever arm 
distance from the center of rotation to the LVDT, to output an angular displacement. 


*It is assumed this data does not contain excessive noise. Otherwise, it is 
best to not correct the error at all. 
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A typical LVDT manufactured by Schaevitz has a linear range of +.01 inches. 


The maximum IMC mirror linear travel equals +.0075 inches assuming an angular 
travel of +10 mrad. For these conditions, the LVDT is linear to +.25% of the IMC 
travel. 

The LVDT is being used for two applications. 

1) It measures the relative angles between the reference platform and the 
telescope (mass-stabilized configuration only). 

2) - It measures the relative angles between the IMC mirror and the telescope 

(for all configurations). 

• Since the reference platform operates about the null during steady state, any 
LVDT non-linearities will vanish as the angular deviation approaches zero. 

For the IMC mirror, the angular deflections are small because the IMC loop 
only responds to disturbances above 10 hertz. For the worst case, however, lets 
assume the mirror must remove the entire 12 urad of base motion RMS error. 
Therefore, the LVDT non-linearity is computed as 
®LVDT = (.0025)(12 urad) = .03 urad 
for each axis. 

6.3.4.4 Gyro Noise 

The GI-G6G gas bea. mg rate-integrating gyro built by Northrop is a strap- 
down space-qualified device with a life expectancy of years. It has a noise 
bandwidth of 5 kilohertz with an RMS error of 0.3 urad. The maximum useable band¬ 


width of 400 hertz results in a gyro noise error of 


© = [400(4^) l 1 1' 1 - .084 urad 

gyro 5000 


for each axis. 
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6.4 Servo Configuration Analysis Summary and Conclusions 


This section summarizes the error analysis and spacecraft reaction torque 
results for each servo configuration. It delineates the receiver and transmitter 
optical LOS errors and examines the overall system pointing error due to contribu¬ 
tions from both sources. Finally, a comparative summary table lists the pros and 
cons of each configuration to assist in recommending a servo approach. 

6.4.1 Error Analysis Results 

The error analysis results generated in Section 6.3 are summarized for 
each configuration in Figure 6.4-1. 


Figure 6.4-2 reduces the data from the above table to compare the servo con- 
figuration performances. It adds the la LOS tracking error to the point ahead 
and dynamic following errors to obtain a total LOS pointing error for each optical 
system. • 


Configuration 

Aperature 

Mean 

lo LOS 

Total 


Diameter 

LOS Error 

LOS Error 

LOS Error 


(in.) 

(urad) 

(urad) 

(urad) 

Gyro-Stabilized 

6 

.078 

.63 

i .71 1 


20 

.078 

.49 

.57 

Mass-Stabi1ized 

6 

.078 

11.65 

11.73 


20 

.078 

11.65 

11.73 : 

Complementary Filter 

6 

| .066 

.93 ! 

1.0 


20 

.066 

.52 

i .59 


Figure 6.4-2: Two-Axis Error Analysis Summary 


The tabulated data -oncludes that the mass-stabilized configuration does 
not satisfy the system LOS error budget requirement of 1 urad. The large LOS 















Error Source 


Sensor Noise 


Base Motion 
Disturbance 


Gimbal Friction 


Point Ahead 


LVDT Non- 
Linearity 


Gyro Noise 


Figure 6.4-1: Two-Axis Error Analysis Results 



Configuration 

Aperture 
Oiameter 


(In.) 

GS, MS 

6 


| 20 

CF 

6 


20 

GS 

6 & 20 

MS 

6 & 20 

CF 

6 & 20 

GS 

6 & 20 

1 MS 

6 & 20 

CF 

6 4 20 

GS,MS,CF 

6 & 20 

GS.MS 

6 & 20 

CF 

6 & 20 

GS,MS ,CF 

6 & 20 

GS 

6 & 20 


Azimuth 
LOS Error 
(urad) 


Elevation 
LOS Error 
(urad) 


Total lo 
LOS Error 
(urad) 





.56 

.17 


.475 



.802 

.240 


.475 


.488 


.004 


.488 


.462 



.462 


.022 


.031 


.038 


6.79 


9.45 


11.64 


.021 


.029 


.036 




.042 


.08 


.08 


.113 


Mean 


LOS E 

rror 

(urad 

) 


28 

.040 

09 

.013 





























































error was shown In Section 6.3.3 to be caused by its poor error rejection per¬ 
formance for low frequency disturbances. In addition, the friction disturbance 
is predominantly low frequency for u equal to 500 zero crossings per second. 

It was suggested that magnetic bearings be used to reduce the friction 
torque disturbance. In consulting HAC experts, it was pointed out that magnetic 
bearings are cumbersome to implement because backup roller bearings are required 
replacements to withstand excessive loading and power failure. 

The gyro-stabilized configuration shows the best LOS error performance for 
both 6 inch and 20 inch receivers. Its shortcomings are reliability, and the 
limited bandwidth of the gyro. For this study, the highest frequency error that 
the gyro must sense was base motion that occurred at 125 hertz which approaches 
its bandwidth limit. For higher frequencies, the gyro must be complemented by a 
higher bandwidth device such as an AOS, which adds further complexity. 

To simplify the design, NASA/GSFC suggested that the gyro be replaced by 
an ADS. The AOS is a high frequency device operating between 2 and 1000 hertz 
(Systron Conner). If only the ADS were used to measure iner'ial LOS position, 
the track loop bandwidth must be increased to 2 hertz to atteruate all errors 
below the ADS's 2 hertz limit. This approach improves reliability since the ADC 
is not a spinning device. 

The complementary filter^ LOS tracking accuracy is inferior to the gyro- 
stabilized configuration because it transmits a greater amount of sensor noise to 
the LOS. This additional LOS noise is due to the wide bandwidth track loop that 
eliminates the need for inner loops to attenuate base motion disturbances. 

However, by doubling laser power to increase the sensor SNR to 10, the 
complementary filter LOS errors are reduced to 

.82 urad 6 in. dia. 

.55 urad 20 in. dia. 

which is comparable to the gyro-stabilized results. 
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6.4.2 Reaction Torque Results 

The reaction torques to the spacecraft are predominantly caused by torque 
loads applied to the 2-axis gimbal. These reaction torques were generated for 
the friction and kinematic coupling disturbances via the computer simulation that 
generated the LOS errors. 

The simulation outputs azimuth and elevation gimbal RMS accelerations which 
are converted to RMS torques by multiplying them by their respective inertias. 

This data was generated for each configuration and tabulated in Figure 6.4-3 . 



•RESULTS FOR 500 ZERO CROSSINGS PER SEC 

f 

Figure 6.4-3: Reaction Torque Summary , 

The complementary filter RMS reaction torque is somewhat smaller than the 
other configurations. Upon observing the transfer functions relating gimbal 
angle to base motion disturbance, the complementary filter transfer function 

rolls of at 1 (at 10 hertz) verses the other configurations roll off at 1 

s 4 s 2 

(at 10 hertz). This means that the complementary filter coarse gimbal removes 

less LOS error than the other configurations resulting in less reaction torque. 

Figure 6.4-4 shows the spacecraft reaction acceleration PSD resulting from 

the base motion disturbance provided by NASA and applied to the gyro-stabilized 

configuration. The forward sum of this curve times the gimbal inertia results 

in the RMS value of .107 ft-lb appearing in the above table. 
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6.4.3 Comparative Summary and Recommendations 

* + 

The table In Figure 6.4-5 summarizes the advantages and disadvantages of 
each servo configuration. It seems appropriate to review them in light of the 
three requirements stated In Section 6.1.1 

1) pointing accuracy 

2) spacecraft reaction torque 

3) design simplicity and reliability 

For pointing accuracy, the gyro-stabilized approach outperforms each of 
the other configurations. It is able to satisfy the 1 urad 2-axis LOS pointing 
error budget with half the required power than the complementary filter approach. 

The mass-stabilized approach was unable to meet the error budget due to excessive 
LOS error residuals caused by friction. 

The spacecraft reaction torque results are very similar for all configura- 
ions. The complementary filter has more flexibility in splitting the LOS error 
by merely changing the corner frequency of its low-pass filter. This characteristic 
is advantageous only if the actual base motion disturbances are larger than the 
predicted values. Under this condition, the IMC mirror travel may be inadequate 
requiring the corner frequency to be increased. A problem caused by underestimat¬ 
ing the base motion is that the complementary filter track loop bandwidth must be 
increased to attenuate the additional disturbance. This in turn increases the 
sensor noise perturbing the LOS. For the gyro-stabilized approach, increasing 
the IMC loop bandwidth has no noise effects on the LOS. 

The complementary filter provides the most reliability of all the configur¬ 
ations because of its simple design. However, it was mentioned that the gyro- 
stabilized reliability could be improved if the gyro were replaced by a more 
space-qualified sensor such as the AOS. This change decreases complexity by 
eliminating excitation signal generators and demodulators associated witn gyro 
implementation. 
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The above results Indicate that the gyro-stabilized configuration provides 
the best tracking approach for the laser communication problem. Its reliability 
was Improved by replacing the gyro with an ADS as the Inertial LOS reference. 

Its design has been Implemented by HAC on numerous pointing and tracking systems 
and satisfies all system requirements. 










6.5 Appendices 

6.5.1 Justification for Soft-Mounted IMG Mirror 

Figure 6.1-2 in the text shows the IMC mirror to be soft-mounted to the 
telescope by means of a single post flexure. This mounting isolates the optical 
train from any telescope vibrations that may enter the loop. To illustrate its 
importance, a simple servo block diagram of the IMC loop is presented in Figure 



Figure 6.5.1-1:IMC Servo Block Diagram with Soft-Mounting Dynamics 
where, = DC gain 

Gjmc = IMC loop compelsation 

5 n = flexure damping 

u) n = flexure natural frequency 

The closed loop transfer function relating the LOS angle, ei_QS» t0 the base motion 


disturbance angle, Og, is 


LOS = 
9 B 


25 i.. S + u 
r n J n 


S'- + 2? u S + in ^ * K G 

n n n IMC IMC 


( 6 . 5 . 1 . 1 ) 


The IMC mirror becomes more isolated from the base (telescope) as 
becomes small. In the limit, 
lim e L0S = 0 

V° ®B 

and 9 LL )s decouples from 9g. 


( 6 . 5 . 1 . 2 ) 








The IMC mirror becomes more rigidly attached to the base as « n becomes 


large. In the limit, 

lim 9 LOS = t (6.5.1.3) 

W n ->® 0 0 

and 0 los couples directly with0g. 

Hence, we see that isolating the IMC mirror from the telescope isolates 
the optical train from any disturbances to the mirror base. 

6.5.2 Derivation of Servo Control Equations 
6.5.2.1 Type II Loop Gain and Compensation 

Given a Type II loop, we wish to determine the gain and compensation to 
provide ample phase margin at the open loop crossover frequency, f c . 

Without compensation, a Type II loop has zero phase margin ( $ = -130° at 
crossover). It is necesssary to add phase lead at f c in order to improve loop 
stability. This is done by incorporating a lead/lag network into the loop that 
is expressed as 


G c (S) = a_(6.5.2.1.1) 

±_ + 1 
R a 

with s = Laplacian variable 
a = compensating zero 
R = lead ratio 

We wish to derive the value of 'a' such that the lead angle, <*> c is a 
maximum at f c . We first express equation 6.5.2,1.1 in the frequency domain. 


_ , . * 1 

G c U“c) = _§_ 

j—C + 1 
Ra 


( 6 . 5 . 2 . 1 . 2 ) 


Multiplying numerator and denominator of 6.5.2.1.2 by the denominator's 


complex conjugate yields, 

“c ■ “C* - fti c 

G (ju.. ) = 1 + Ra? + 1 ( Ra ) 
c c —- 


(6.5.2.1.3) 
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(6.5.2.1.4) 


» 


and solving for the phase angle results in 

Z.G (jw ) = 4» = tan" 1 V*" 1 ) — 

c c c Ra + u> c z /a 


where # c = the desired phase margin at f c . 

Taking the tangent of both sides of 6.5.2.1.4yields 

tan $ = r.^" 1 — 

c Ra + ^/a 


(6.5.2.1.5) 


To solve for the maximum value of c requires taking the derivative of $ c 
with respect to 'a', setting it to zero, and finally solving for 'a'. 

Therefore, 


or 


d_ tan 0 = _ 1 (6.5.2.1.6) 

da c cos^ 

c 

p 

d ^ cos t{ ^ 

“da = (Ha 2 l w 2)2 [(Ra 2 + <-> c 2 ) « c (R-l)- “ c (R-l)(2Ra 2 )] = U (6.5.2.1.7) 

There are two ways for equation 6.5.2.1.7 to equal zero 


1) ^ 9U degrees 

2) the bracketed expression equals zero 
Since we do not wish to limit $ c to +_ 9U degrees, 

(Ra 2 + w c 2 ) = and 



( 6 . 5 . 2 . 1 . 8 ) 


Having solved for the compensating zero, we wish to now solve for the DC 
gain (Ka) which must satisfy the requirement that the open loop gain at the cross 
over frequency equals unity, i.e.. 


Ka 

jw f 

V + 1 

G(J <*> c ) 1 = 

4 . 1 

Ra + 1 


= 1 


(6.5.2.1.9) 


Upon substituting equati on 6.5.2.1.8 i nto 6.5.2.1.9 and simplifying, yields 


G(j w c ) 


. iy R u>' /2 = i 

w r 2 (R-l +1 ) i n 


(6.5.2.1.10) 
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Finally, solving for K a results in, 

l| 2 i} 

i K = «c 
a R1 / 2 

6.5.2.2 Complementary Filter Compensation 

This section derives the compensation networks for the complementary filter 
configuration. The servo philosophy is to split the track error signal into its 
high and low frequency components such that the high frequency error is attenuated 
by the IMC and tne low frequency error is attenuated by the coarse gimbal. A 
servo block diagram depicting this scheme appears in Figure 6.5.2-1 





where Gqq = coarse gimbal compensation 
Gimc s mirror compensation 
M = telescope magnification 

^ = IMC mirror flexure damping ratio 


U) 


m = effective IMC mirror spring rate 


9 


We wish to determine and Gj^ such that the open loop transfer function, 

LOS , equals an optimally compensated Type II track loop (Appendix 6.5.2.1), ie., 
9e 


0 X — + I 

e L0$ . K a a 


ee 


2 s 


(6.5.2.2.1) 


Ra + 1 


The open loop transfer function was computed using Figure 6.5.2-1 to be. 
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Since there are two unknowns and one equation, G^g was defined to be the gain K a 


times a fourth order low pass filter 

K j 

CG / ~: \ 2 



(6.5.2.2.3) 


where f =“n = corner frequency ( 1U hertz) 
n 

With Gqq defined, equation 6.5.2.2.2is set equal to 6.5.2.2.1 and soived for G^g 
to yield. 

To 0^+1 

a . 1 


6 IMC = K a 


s2 +W 4 M 2 
S 2 


(1 


L_ + 1 (siL + 2C n S + 1)2 M 
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Ra 
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The resulting servo block diagram appears in Figure 6.5.2-2 

6.5.3 Interpretation of Statistical Coulomb Friction Model (See reference 1 for 
derivation of model) 

Coulomb friction enters the CG loop as a result of the relative motion 
between the gimbal base and its bearings. Figure 6.5.3.1 shows this effect in a 
simple servo block diagram where 

“t = inertial gimbal angular rate 

“ v = inertial gimbal base angular rate 

• • 

F = friction acceleration 
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Figure 6A.3-1: Friction Acceleration Model 


In this model, it is assumed that the gimbal rate, “y, is zero, and the 
only input is the base motion rate, w v . Therefore, the relative motion between 


the gimbal base and its bearings is« v 
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w v is a random signal that can oe represented by the upper curve in Figure 
6.5.3.2, The friction acceleration corresponding to this rate is the constant 
amplitude A, that has a polarity equal to the polarity of w v . The lower curve in 
Figure6.5.3-2 represents the applied friction acceleration. This curve, known in 



Figure 6.5.3-2: Base Motion Rate and Friction Acceleration vs. Time 
probability theory as a “random telegraph signal", consists of a sequence of 
steps of constant amplitude (+A) that occur at random time intervals, AT. The 
probability that the signal changes sign exactly K times in T seconds is given by 
the Poisson distribution, 

P(K) = (uT) k e~ uT (6.5.3.1) 

k i 

where u = average number of zero crossings per second. 

The maximum value of P(k) occurs when k equals its expected value, uT (the 
average number of zero crossings). As uT increases, the probability of uT, P(uT), 
decreases and the variance about the mean increases. This is illustrated by 
plotting P(k) for two values of uT equal to 3 and 6 in Figure 6.5.3-3 

Increasing the average number of zero crossings implies that the base 
motion rate has frequency components spread over an increasingly larger portion 
of the spectrum. The larger frequency spread increases the sample size and 
decreases the probability of occurance for a given frequency value. 

















Figure 6.5.3-3: Poisson Distributions for uT * 3 and 5 


These probability distribution characteristics appear in the autocorrela 


tion function derived in reference 2 to be, 

R( T ) = A 2 e" 2u l T I 

and is plotted in Figure 6.5.3-4 fo,r two values of u. 

R(T) 


( 6 . 5 . 4 . 2 ) 
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increasing u 


Figure 6.5.3-4: Friction Acceleration Autocorrelation Function 
As u increases, the friction acceleration becomes less correlated with 
time (dotted curve). Increasing u subsequently reduces the PSD amplitudes as 
displayed in Figure 6.5.3-5. And, in the limit, 


1 i m S(w ) = 1 i m u 


= u 




(6.5.4.3) 
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Figure 6.5.3-5: Friction Acceleration PSD Variations with u 
During our final briefing to NASA/GSFC, there was concern that as u 
increases, more frequency components of friction would exceed the servo loop 
bandwidth and cause the LOS error to increase. This statement is correct but is 
not reflected in the data of Figure 6.3-8. Since this data was generated by a 
computer, integration limits are required. As u was increased, that portion of 
the friction disturbance PSD that exceeded the upper integration limit, was trun¬ 
cated. Since the friction angular displacement PSD (S p ( «)) rolls off at 1_ 

F s 6 

(1 times S.. (w) in Figure 5.5.3-5', increasing u truncates a greater portion of 


s 4 Op 

S fi (w) resulting in a smaller RMS LOS error. 

F 

6.5.4 Orbital Mechanics Equations to Compute Relative LOS Angular Accelerations 
This section develops the LOS angular acceleration equations between the 
two earth-orbiting satellites. One satellite (GEO) is in a geosynchronous orbit 
in the equatorial plane while the other (LEO) is in a low earth orbit at an 
inclination of 8.2 degrees from the spin axis. The orbital geometry is illustrated 
in Figure 6.5.4-1, 

An inertial orthogonal coordinate system is defined by the unit vectors 
(1, J, K) where T and J are in the equatorial plane and K passes through the 
Earth's spin axis. 
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LED ORBIT 



For relative LOS orbital computations, an in-plane polar coordinate system 
is most desirable where position, velocity and acceleration vectors are resolved 
into components along (u r ) and perpendicular to (u ) the LOS vector. 



Figure 6.5.4-2s In Plane Polar Coordinate System 
The following list defines all parameters used in the derivation. 
7|_, ?g s GEO, 6E0 position vector in inertial coordinates 

Vl, Vy = LEO, GEO velocity vector in inertial coordinates 

“q = GEO, GEO angular velocity vector in inertial coordinates 
The LOS vectors are defined as 

rQ - tg - r|_ LOS position vector 

V D s V G ‘ V L GOS velocity vector 

V n = LOS acceleration vector 






















The LOS position, velocity and acceleration vectors are represented in the 
polar coordinate system as follows: 



r D = r D u r 



• 

Vo * ro u r + “o r D u ® 


and 

J. . • • p — * • — 

Vq 8 ( r o - w |) r o) u r + («*>0 r D + 2 “D r 0^ u e 
, .( V D- f 0 2 ) 1/2 

(6.5.4.1) 

where 



D r 0 


and 

u 

(6.5.4.2A) 


u » —— 
r r D 

„ . V D ’ ^D u r (6.5.4.2B) 


u “D r 0 

We wish to solve for wq in equation 6A.4.I. By taking the dot product of 

6A.4.1 with the unit vector, u , results in 

9 

• • • 

V D • u = + 2(^0 (6.5.4.3) 

Substituting equations6.5.4.2A and B into 6 . 5 . 4.3 and solving for wq yields the 
desired expression 


& = rp Vn.VD --,r D V D .r D - 2^ 

D a, D r D r D 
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7.0 TRACKING AND BEACON LINK 


7.1 INTRODUCTION 

In the given study scenario of a low earth orbiting (LEO) satellite 

transmitting a high data rate link to a .satellite In geosynchronous (GEO) 

orbit, the most severe tracking burden Is placed on the LEO optical tracking 

system. The LEO has a 15 cm (5.9 inch) receiving telescope, and must generate 

2 

the pointing commands to direct the narrow communication beam (lOurad, 1/e ) 
to the GEO with about one microradian rms accuracy. By contrast, the GEO 
optical tracker has a receiving area of 50 cm (19.7 inches), and has more than 
a factor of two reduction In required pointing accuracy due to the larger 
beamwldth of the low data rate beacon (GEO to LEO). 

This study then focuses on optical tracking at the LEO, since If this can 
be accomplished, the GEO optical tracking can also be accomplished, due to the 
greater optical collecting area and reduced pointing requirements. 

7.2 NOISE EQUIVALENT ANGLE OF QUAD CELL TRACKING SENSOR 

The pointing and tracking servo analysis of section 6.0 requires 
knowledge of the tracking sensor noise equivalent angle (NEA) in order to 
determine the relative contribution of the optical sensor to the residual 
pointing error. The NEA can be considered as the Inverse of the rms signal to 
noise ratio referred to angle space. 

For a quad cell with the upper right quadrant designated as quadrant 1, 
and numbering clockwise, the azimuth and elevation tracking error signals are 
given by 


e z (t) * [Xj(t) ♦ X 2 (t)] - 

Cx 3 (t) ♦ x 4 (t)] 

Azimuth 

e L (t) x CX 1 (t) ♦ X 4 (t)] - 

[x 2 (t) ♦ x 3 (t)] 

Elevation 
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where X,(t) Is the electrical output of the 1 th quadrant. The assumptions 
made In this analysis are: 

1) Each quad element has Identical statistics, and Is statistically 
Independent of the other elements. 

2) The mean output of a quad element due to signal Is E[X^] a eG K<. 
where e Is the electronic charge, G Is the average detector gain, and 
K are the signal photoelectron counts. Since K $ are the signal 
photoelectron counts per element, must be computed on tne basis 
of one-fourth of the optical signal power Incident on the quad array. 
The noise mean Is assumed to be zero (non-zero mean noise Is 
discussed at the end of this section). 


3) The rms noise per element Is given by 



e 6 / F N K S + K N 


where F^ Is che detector excess noise factor, and are the 
total noise photoelectrons. For an RCA silicon reach-through 
avalanche photodiode, the excess noise can be approximated for high 
gain by 


f,'!’ 0-02 S 

where G is the mean detector gain. The excess noise is a measure of 
the randomness of the detector gain, and is unity for a detector with 
no gain. 

The noise photoelectrons are given by 
s F N * K B + K 0* + K DS * K T 
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where Kg are the photoelectrons due to optical background 
radiation, Kg are the photoelectrons due to gain-dependent dark 
current, K^ are the photoelectrons due to gain-independent dark 
current, and Ky are the thermal noise photoelectron counts. Since 
these photoelectron counts are for a single quad element, Kg must 
be computed on the basis of one-fourth of the total background power 
incident on the array. 

4) The Airy spot on the quad array is approximated by a uniform 
intensity disc with radius 

r Q = 1.22 (x/D r ) f 

where D is the effective diffraction limited receiving telescope 

diameter and f is the focal length of the focusing lens. 

We will analyze the elevation error signal, and the azimuth case is 
completely analogous. The expected va'iue of the error signal is 

E M* E { Xl>X 4)- ‘M} 

= eG ^ n /hJ *s T ^ A i + A 4 * *2 " A 3* 

where n is the detector quantum efficiency, hv is the photon energy, I is 
the signal intensity, t is the integration time, and A- is the area inter- 

t h * 

cepted by the spot in the i quadrant. 

Referring to the focal plane geometry of Figure 7.2-1, where tfi ^ 
and are the offset angles from null in azimuth and elevation, 

respectively, we obtain 

A 2 * A 3 * 4- r 0 6 - f h V r o * f2 *< 2 

A i , »c' r M r ! ,,f h'l'l - f2 *‘ 2 
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where o is defined in the figure. Taking the difference, we have 


A 1 + A4 - (A 2 + = rrr 0 


1 —— cos 

IT 


*(-e) 


77 J 1 - < f *|/r 0 > 


r 


This result agrees with that obtained by Gagliardi and Karp, Optical Com ¬ 
munications , page 391. Using the small argument approximation for the inverse 
cosine (to linearize the loop). 


COS"* X ~ ir/2 - X , X << 1 


the area difference becomes 

Aj * A, - (A 2 ♦ A 3 ) = (.r*) -t- 


(1.22 X/D r ) 


The expected value of the elevation error signal is then 
E (e L (t)] . eG(„/M I s (.4 -t- Jr *f TN - T 

16 eQP' 

* ir ebK S 11.22 x/D p ) 

The second step can be accomplished since I^r^ is the total power in 
the spot (intensity x area), which generates 4K $ signal photoelectrons (K $ 
in each element). 

Now examining the variance of the elevation error signal, due to the 
assumptions of independence and identical statistics per element we have 


Var e L (t) = Var Xj_ ♦ Var X 4 ♦ Var (-X 2 ) + Yar (-X 3 ) 

* 4 v * 4 * V 
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Figure 7.2-1. Noise Equivalent Angle Focal Plane Geometry 
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The mean signal to rms noise of the «n< 7 ':e measurement is then 


1 R (t) } 

Var Ce L (t)3 


8 _ h _ 

(1.22 x/D p ) 



where the rms signal to noise ratio is defined by 



Note that the azimuth offset angle does not appear in the expression for 
elevation error, so that the azimuth and elevation tracking loops are 
uncoupled. 

The noise equivalent an-.ie (NEA) is thus given by 



(1.22 x/D r ) 


SNR 


rms 




If the imaging system is not diffraction limited, then the quantity (1.22 
x/D r ) can be replaced by the radius of the angular optical- blur. Again, it 
should be emphasized that the SNR is computed on the basis of an element of 
the quad, which sees on average one fourth cf the total optical power incident 
on the array (signal and background). 

Recently, analyses of quad cell tracking have been performed which take 
into account the Airy intensity distribution of the focal plane spot rather 
than using a uniform intensity disc. [G. Tyler and D. Fried, J. Opt. Soc. 
Am., Yol. 72 No. 6, June 1982; R. Gagliardi, et. al., JPL Publication 81-40]. 
The difference with the present analysis is slight, replacing 1.22ir/8 in the 
previous NEA expression by 3n/16. 

It should be noted that if assumption (2) in the derivation for tracking 
NEA is relaxed to allow an identical noise mean in each quadrant, the NEA 
expression does not change. The DC currents due to thermal noise, dark 
current, and uniform background illumination are cancelled in the sum and 
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difference algorithm [GagHardl and Karp, Optical Communications ]. For this 
to occur, however, the four detector elements must be matched in terms of 
gain, responsivity, excess noise, preamp noise, and dark currents. 

The noise means therefore affect the tracking algorithm only to the 
extent they are nonuniform among the quadrants. Assuming the detector-preamp 
parameters can be matched resaonably well, the noise mean nonuniformity is 
primarily caused by optical background. An extended source could be offset in 
the quad FOV, or a star located In a single quadrant. For these cases the 
noise mean nonuniformity appears identical to the signal mean nonuniformity to 
be tracked. The only solutions to this are to overwhelm the worst case noise 
mean nonuniformity with signal, or to track a source with a spectrum shifted 
away from DC (such as a PPM data waveform) so that AC coupling may be used. 
In general, the latter approach is suggested, although if a gyro- stabilized 
tracking servo approach is used, there may be enough SNR available to overcome 
the noise mean nonuniformity. 

It should also be mentioned that shot noise nonuniformity between the 
quadrants (as in the case of a star in one quadrant) does not bias the 
tracking error. This is because the error variance of the angular measurement 
is the sum of the noise variances in the four quadrants. Thus a star of a 
given irradiance in one quadrant contributes one fourth the angular error 
contribution of an extended source with the same irradiance level present in 
all four quandrants. 

To compare alternative servo tracking system designs such as gyro- 
stabilized or complementary filter, it is desirable to know the dependence of 
the rms SNR on sample rate. Assuming the integrator output is directly 
sampled, the sample rate (f ) is given by f = 1/t , where t is the 

3 J 

integration time. If the noise is not dominated by thermal noise (valid for 
both PIN and APD detectors for t > 1.0 ysec) the noise photocounts can be 
written as K^ = n^ t, where n^ is the noise count rate (e“/sec) and r 
is the integration time. Similarly, the signal counts can be written 
K$ = n s t » where n~ is the signal count rate (e“/sec). The expression 
for SNR then becomes 


7-7 








1 



n 




F N ", * "» 


1 


Since is proportional to laser power, we conclude that for a given laser 
power, the rms SNR is inversely proportional to the square root of the sample 
rate. This is important for comparing the stabilized approach and the 
complementary filter approach, since the complementary filter approach 
requires a much higher sample rate, and therefore more laser power for a given 
SNR. 

7.3 QUAD CELL DEAD ZONE LOSS 

The typical quadrant tracking detector ’has a gap between the quadrants 
which generates no photoelectrons in response to optical intensity - in 
effect, a dead zone. ( If the detector is reverse-biased, some photoelectron 
events may occur in the quadrant gap, but for a worst cas.e analysis the gap 
will be considered as "dead".) When a quadrant detector has the optical spot 
centered, as in closed loop tracking, the spot has a vertical and horizontal 
stripe of dead zone obscuration. The fraction of total spot intensity lost 
due to the dead zone can be approximated as the ratio of area obscured by the 
dead zone to the total area of the spot. (This assumes a uniform intensity 
approximation to the Airy spot). 

Referring to the geometry of Figure 7.3-1, d is the width of the dead 

zone, and the diffraction limited optical blur is 2.44 -p- f , where D r is 
the receiving optics diameter and f is the focal length. r (If not diffraction 
limited, 2.44 x/D r can be replaced by the angular subtense of the optical 
blur). 

Let n denote the fraction of intensity lost in the dead zone gap. n is 
given by the ratio of obscured area to total spot area. Then n becomes 










1 

T 



GEOMETRY FOR DETERMINATION OF INTENSITY LOSS 
DUE TO QUADRANT GAP 


Figure 7.3-1. Geometry for Determination of Intensity Loss Due to 

Quadrant Dead Zone. 
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n « 


2 (2.44 x/D ) f d-d 2 

r 





where y Is defined by 

y = (d/f) / (x/D r ) 

r is proportional to the ratio of the angular subtense of the quadrant gap to 
the angular blur of the receiving optics. For small y, n ~ 1.04 y. Figure 
7.3-2 shows a plot of n versus y for n from 1 percent to 80 percent. 

Until realistic numbers are inserted, it is not apparent that the dead 
zone loss can pose a serious design constraint. Suppose we choose a quad cell 
with a 4.4 mil (112 pm) gap (for example the Hughes HP1N 200 Q quadrant 
detector), x = 0.85 pm, 0 R = 20 inches, and decide to allow only 1 percent 
loss. This would require a focal length of 6.6 kmJ 

Therefore very small quadrant dead zones are desirable, and to obtain 
practical focal lengths some intensity loss must be tolerated. The RCA 

quadrant avalanche photodiode (Section 7.5) has a smaller 75 pm dead zone. 
For D r = 6 inches (15 cm), x = 0.85 pm, d = 75 pm, and 50 percent dead zone 
loss a focal length of 26.5 meters is required. A program is underway at RCA 
(in Montreal) to reduce the dead zone to 30 pm. This improvement would reduce 
the required focal length to 10.6 meters, which is approaching a manageable 
length. Other programs are needed to produce small dead zone quadrant 

detectors. 

It might be suggested that dead zone loss could be decreased by defocus- 

ing the spot, making it larger in relation to the quadrant gap. But, as shown 

in section 7.1, the noise equivalent angle (NEA) of the tracking sensor is 
proportional to the optical blur size and inversely proportional to SNR. So 
if the spot were increased due to defocusing, the SNR would have to be 
increased to get back to the same tracking performance. Either way, the dead 
zone loss costs laser power. 
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7.4 QUAD CELL CLOSED LOOP FIELD OF VIEW 


It Is of Interest to determine the dynamic tracking range over which a 
quad cell tracker can remain closed loop, since If tracking Is lost some form 
of reacquisition must take place. If the optical spot is centered on the four 
quadrants, displacement of one spot radius will cause It to no longer overlap 
one or more quadrants. Unless all four quadrants have at least some 
Intensity, a quantitative error signal cannot be generated for both azimuth 
and elevation. Thus the maximum radial displacement for which closed loop can 
be maintained Is half of the angular blur. Let D $ denote the optical blur 
diameter In the focal plane, with angular subtense D s /f, where f is the lens 
focal length. Since the closed loop diplacement In the focal plane Is 
♦ Dj/2, we can consider the closed loop field of view to be 


e 


CL 


* D s /f 


Rewriting the expression for Noise Equivalent angle (NEA) from Section 
7.2 to call D^ = 2.44 {x/D r )f, the expression for the residual rms 
tracking error becomes 


a NEA 3 (ir/16) (D $ /f J/SNR 


rms 


= U/I6)(e n /SNR1 

ll rms 


Since a, 


is a radial error. If we call o. ICA = 2i 


as the residual 


NEA "NEA = tC NEA 

tracking error field of view, the closed loop tracking dynamic range is given 

by 


9 CL 

9 nea 


= (8/ir) SNR 

rms 
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Thus the closed loop tracking dynamic range Is proportional to the SNR. 
For SNR above about 50, other practical limitations to the residual tracking 
error begin to dominate, such as element nonuniformities, amplifier dynamic 
range, and focal spot asymmetry. The minimum tracking error practical to 
obtain Is about 0.01 (x/D), so that the maximum closed loop dynamic range for 
high SNR Is about 127. 

7.5 PHOTODETECTION ANALYSIS 

7.5.1 Means and Variances (Gaus*1 anj. This section will relate the signal 
and noise means and variances to the physical parameters of the optical 
receiver. This preliminary development is required in order to evaluate and 
compare the various detection strategies. 

The simple circuit model is shown In Figure 7.5-1. The receiver model 
described here Is a high impedence (integrating) front end, which is the 
lowest noise approach. Signal and background optical radiation is Incident on 
the optical detector, and their photocurrents (1 , 1. ) are amplified by 
the detector's internal gain (mean gain G). Also shown is a dark current 
which receives detector gain (1 Q ), as well as a gain-independent dark 
current (i^). The preamplifier has an equivalent thermal noise current 
source (i-r). The detector and preamp each have a capacitance, and the sum 
capacitance is C = ♦ C . R, is the Input resistance of the preamp, 

which is the physical source of the thermal noise. The Input resistance can 
be generalized to an equivalent temperature which includes the other noise 
sources of the amplifier. This equivalent temperature (T e ) Is related to 
the preamplifier noise figure (F) by 

T eq = (F-l) 290 


7-13 










ORT w *il.K 
OF POv'. Q •*-** ' l 




— 1 

G 's 

) 0 

G'u 

) 0 

G ‘o 

) a 


-ft 

k II 

‘.II 

_Ll 


* 

y » 7 « 

1 

>l l 

1 

i 




Slot 

r T 


OR 


■ 

WOR 0 



DECISION 


DETECTOR 


L 




INTEGRATOR 


Signal current 

CURRENT DUE TO iACKGROUNO RADIATION 
GAIN-OE’ENDENT Dark CURRENT 
GAIN-INDEPENDENT DARK CURRENT 
PREAMP EQUIVALENT THERMAL NO>SC CURRENT 
DETECTOR CAPACITANCE 

c»c d *c 

AMPLIPIER capacitance 
AMPLIFIER input resistance 
AVERAGE DETECTOR GAIN 





Direct detection receiver model 


Figure 7.5-1. Direct Detection Receiver Circuit Model 
















To find the equivalent thermal current source, recall that a resistor 

R l at a temperature T will generate a one-sided voltage spectral density 

4kT eq R L V 2 /Hz. or a two-sided density 2kT e q\ (where k Is 

Boltzmann's constant). The power delivered to R^ by this voltage source is 

2kT A „ W/Hz (two-sided). The power delivered __£o R. by the equivalent 

eq 2 2 L 

thermal current source Is 1 t R L (where 1 t is the mean squared 

current). Equating the expressions for power delivered to R^ yields 





or _ 

I 2 = 2kT /R. A 2 /Hz (two-sided) 
t eq L 

2 

The total mean-squared current noise level (ijoj > two-sided) in - time 
slot containing the optical pulse is 


i 


2 

TOT 


G 2 1 $ + 


e G 2 f B • e G 2 1 0 



♦ 




2 k T 


~TT, 


eq 



* 


2 

where G Is the mean-squared optical detector gain and e is the electronic 
charge. (The one-sided spectral density of shot noise due to current source I 
Is well known to be 2el.) 

Using Parseval's Theorem to take advantage cf the simple time domain 
response of an Integrator, 


hit) 



0 < t < T 
otherwise 
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wheretls the Integration time, the variance in a pulse slot Is given by 

«S 2 ■ ' T [ gZ ‘S 2 <8* B *V 'os* 2kT e, /e \] 

o 

The variance In a slot with no pulse (o„) Is obtained by setting 1 s 0 

0 \ n s 

in the expression for o*. The mean signal at the Integrator output 
(m $ ) Is given by 


m. 


f G 1 (t) dt = GIT 
o s 5 


If we define K = 1 T/e as the signal photoelectron counts per Integra- 
tlon period, then m * eGK . (For Gaussian statistics and fixed noise 
mean, the mean output for a slot with no pulse can be taken as m n r o, 
without loss of generality. For time-varying background sources the noise 
mean may have to be considered - Sections 7.7, 7.8). 


A quality factor, or SNR, can be defined by 


SNR 





♦ K 


N 


* 


where includes all the noise photoelectron counts, and F N Is the 
detector excess noise factor. The total noise counts are given by 



F N ^ K B + + K DS + K T 












where by this definition the noise counts depend on the excess noise factor. 
The noise photocounts due to Background radiation, gain-dependent dark 
current, gain-independent dark current, and thermal noise are given by, 
respectively, 

K B * ’b T/ ^ e ’ K 0 = *D T ^ e * K 0S = *D$ T ^ e g2, 
and 

K t = 2 kT^ v e h\. 

Rewritten in terms of photoelectron counts, the means and variances become 


m = e GK 



o 


(The noise mean will be addressed separately where required.) 



e 2 G 2 I F n ( K s * K b • K d ) * K 0$ ♦ K t 


[' 

[' 


= e G (F^Kj + K n ) 


e 2 G 2 I F n ( K b * K 0 ) ♦ K ds * K t 


= e 2 G 2 K 


N 


The photodetector excess noise factor, F,,, is a measure of the random¬ 
ness of the detector internal gain. It is defined by the ratio of the mean 
squared gain to the mean gain squared, i.e. 



Var (G) 


where Var(G) is the gain variance. For a reach-through silicon avalanche 
photodiode, the excess noise factor can be approximated as 
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F„ . k,„ G M2 - 1/GH1 - k eff ) 


where k is the electron Ionization rate ratio. For an experimentally 
measured k ^ * 0.02, and with reasonably high gain (G „ 100), the excess 
noise becomes 


F N ■ 2 + 0-02 G 

The photoelectron counts due to background radiation are given by 
K B = L R } e R A R N R b N F n £ 

where L„ is the background radiance, e R is the receiver planar field-of- 

view, A R is the receiving area, N R is the efficiency of all the receiver 
optics, b is the bandwidth of the optical bandpass filter, Np is the optical 
filter efficiency, n is the detector quantum efficiency, t is the slot 

integration time, and hv is the photon energy. If a localized object is 

p 

entirely within the receiver field-of-view, L„. [U/4)e R ] must be 

replaced by L,, the irradiance of the object. Care must be taken to 

properly scale objects which are only partially within the receiver field-of- 
view. 

7.5.2 Minimization of Thermal Noise . As shown in the previous section, the 

thermal noise photoelectron count (K T ) is given by 

2 2 

K x = 2 k T „ T /e G R. , where k is Boltzmann's constant, T is the 
T eq L eq 

equivalent temperature of the preamp, t is the integration time, e is the 
electronic charge, G is the mean gain of the optical detector, and R, is the 
input resistance of the preamp. As can be easily seen, it is desirable to 
make R. as large as possible to reduce the thermal noise photocount. 
However, the circuit must respond faster than the sample rate (1 /t ). Let C 
represent the total capacitance of the detector-preamp combination. The R^C 
product then imposes a bandwidth limitation on the receiver. If the circuit 


7-18 









Is modeled as a single pole RC filter, the 3 dB frequency Is fp = 1/2 itR.C. 
Setting the 3dB frequency equal to the sample rate (1/t ), we obtain fp = 
1/t , or T /R, = 2irC. Thus the thermal noise count becomes K T = 4» kT 

j « ■ ■ ■ 

C/e G 4 . Note that the thermal noise Is now independent of _r_. Thus the 
receiver has two fundamentally different noise sources, the background 

radiation and dark current noises which are t - dependent, and the thermal 

noise which is t - independent. It should also be noted that the minimum 

thermal noise is fixed by the capacitance of the detector - preamp combination. 

7.5.3 Optimum APD Gain . Recall from Section 7.4.1 the SNR was defined as 

SNR = K\ / [Fy (K s * K b * K 0 ) ♦ K os * K t ] 

and from Section 7.4.2 the minimized thermal noise term 

K t = 4* k T C / e 2 G 2 = «/G 2 , 

where s is now defined to absorb irrelevant terms. The excess noise factor 
was shown in Section 7.4.1 to be closely approximated by F., = 2 ♦ 0.02G. 
Thus it can be seen that as detector gain is increased, thermal noise is 
reduced, but excess noise is increased. Clearly an optimum value of gain 
exists to maximize the SNR. 

By differentiating the SNR expression with respect to gain, equating to 
zero, and solving, it is straightforward to show that the gain for maximum SNR 
is given by 

3 / -1 

G = ylOO fi / (K$ + Kg + Kq) . However, the gain which 

maximizes the SNR may not, in general, minimize a digital probability of bit 
error. This is because the expressions for error probability of an optical 
communication system are not necessarily dependent on SNR as an explicit 
parameter. 
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7.6 LINK EQUATION 


The relation for propagation of electromagnetic energy to a distant 
collecting area is applicable to RF and optical frequencies. Antenna gain 
values will be quite different for RF and optical system because optical 
apertures are much larger in wavelength units than RF antennas. Also, signal 
energy loss mechanisms will be quite different for RF and optical systems 
because the devices employed to Implement the systems are completely 
ddifferent. Therefore, although the range equation for RF and optical 
frequencies is conceptually identical, in practice the values of parameters in 
the equation are very different. 

The detection analysis results in determination of the optical power 
required at the detector input (P R ). The link equation translates the 
required P R into laser source transmitted power {Pj) From this viewpoint, 
the link equation can be written 

p T * Pr H - ) (A - ) ) mr 

where Dy and D R are the transmitting and receiving telescope diameters, x 
is the wavelength, R is the range, Ny and N R are the overall transmitting 

and receiving efficiencies, L p is mispoint loss, N p is the optical filter 
efficiency, T^ is the atmospheric transmission (if applicable), and T^ 

includes any power splitting. 

The reciprocal of the bracketed squared quantities containing Dy and 

D„ can be considered as the transmitting and receiving antenna gains. 

Because Dy and D R are very large in terms of optical wavelengths, the 

antenna gains will be much larger than for an RF system. However, the large 

antenna gain (narrow beam divergence) is the cause of additional pointing, 

tracking, and acquisition complexities for optical systems. Using 6y = 2.44 

x/D t allows the link equation to be expressed in terms of transmitting 

divergence, where e T is the full width of the main Airy lobe. (The 

2 

comparable 1/e width is about 2/3 of the Airy lobe width.) 
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The signal photoelectrons required from the detector output are given by 

Kj * (n/hv) P^ T 

where n is the detector quantum efficiency, hv Is the photon energy, and t Is 
the effective integration time. 

Table 7.6-1 summarizes the values of some of the link parameters which 
are assumed in this report, unless otherwise specified. 


TABLE 7.6-1. ASSUMED LINK PARAMETERS 


Parameters 

Range 

Receiver Diameter 
Transmitting Optics Efficiency 
Receiving Optics Efficiency 
Mispoint Factor 
Filter Transmission 
Filter Bandwidth 
Amplifier Noise Figure 


Value 

R = 44,500 km 
Dr = 0.15 m 
N T = 0.5 
Nr = 0.54 
Lp = 0.4 
Np = 0.7 
8 = 50 A 
F = 2 

(T eq = 290) 


Optical Wavelength 


x = 0.85 um 















7.7 PIN QUADRANT DETECTOR TRACKING 


Many companies manufacture silicon PIN quadrant detectors, Including 
Hamamatsu, Silicon Detector Corporation, Applied Solar Engery Corporation, 
United Detector Technology, and Hughes Industrial Products Division. A 
summary of commercially available quadrant detectors and their performance is 
given in Table 8.2-2. The Hughes data sheet is shown in Figure 7.7-1, since 
Its performance meets that of the other manufacturers, and has the most 
complete data sheet. The key performance parameters for the Hughes HPIN 200Q 
are: 


Dead Zone: 4.4 mils 

Element Size: 100 mils 

Quantum Efficiency: 83 percent 

Capacitance (per element): 4.5 pf 

.Dark Current (per element): 30 nA 


Since the HPIN 200Q is representative of commercially available quadrant 
detectors, its parameters will be used to calculate a tracking link budget. 

As required by the complementary filter tracking servo analysis (Section 

6.0), we will assume a sensor sample rate of 4 kHz. This is an effective 

integration time of t s 250 u$ec. With 4.5 pf capacitance per element of the 

quad, and assuming the quad element amplifier contributes another 5 pf, the 

total capacitance of the detector-preamp combination is about 10 pf. Using an 

amplifier noise figure of 3dB (Table 7.6-1), the amplifier equivalent 

temperature is T = 290, and the thermal noise counts are K T = 2.0 x 

1 I 

10 (See Section 7.4). For dark current Ig = 30 nA, the dark counts are 

K 0 = 4.7 x 10 7 . 

To determine the photocounts due to background radiation, we must 

determine the field of v' " of the entire quad array. Using dead zone loss n 
= 50 percent, with quadranv gap d = 4.4 mils, and for a 15 cm aperture \/D = 
5.6 urad, a focal length of 40 m is required. (A program to reduce the 
quadrant gap in quad cells is essential.) The quad array field of view is 
then e R = D/f, where D is the diameter of the quad array and f is the focal 
length. Using D = 200 mils and f = 40 m gives e D = 127 urad. 
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Referring to Table 7.7-1, which lists the astronomical sources of 
background radiation, the worst case Irradlance within 127 prad (neglecting 
the sun) would be due to the moon with Lj = U/4) L R = 5.7 x 

10"** WM"^ A'*. For parameters as assumed In Table 7-1, with quantum 
efficiency n = 0.83 and integration time = 250 psec, the background counts 
per quad element are 

K b = 4300 T $ 

where T~ Is the fraction of total Intensity routed to the tracking 
detector. (Note that in obtaining Kg for an element the total background 
power on the array was divided by four.) 

The total photoelecton counts due to shot noise are given by 

k n ■ k b * K o * k t ■ 6 - 7 x 10? 

Note that the counts due to background radiation are negligible, so the noise 
counts are not dependent on the tracking split T~. 

A block diagram of a tracking detection circuit is shown in Figure 7-7.2. 
This circuit can track cw or modulated sources. The output of each detector 
quadrant is amplified, integrated, thresholded, and then fed into the usual 
quadrant detector sum and difference circuit. The sum and difference circuit 
implements th° algorithms (Section 7.2) which generate the azimuth and 
elevation error signals. The detection probability densities are shown in 
Figure 7.7-3, which illustrates that by incorporating a threshold the noise 
mean is effectively removed. The signal mean just has to be enough larger 
than the noise mean to overcome the noise variance, so the detection quality 
depends only on the signal mean and noise variance (both signal present and 
signal absent). The threshold has eliminated dependence on the noise DC level 
(assuming the noise mean is fixed). 

An alternate* tracking detection scheme to that in Figure 7-6 for tracking 
a data waveform is shown in Figure 7.7-4. Each quadrant detector output is AC 
coupled and lowpass filtered, with the filter being sampled at the rate f . 
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These sample values are then summed and differenced (Section 7.2) to yield the 
angular error signals. The lowpass filter bandwidth must be greater than the 
sample frequency. The required sample rate is dependent on the selected 

tracking servo approach, as for example a 4 kHz rate for a complementary 
filter servo, and a 10 Hz rate for a gyro-stabllized servo. 

Since the detector output In AC coupled, this approach cannot track a cw 
source, but only a signal with a spectrum shifted away from DC. This Is 
characteristic of a PPM data waveform, for example. The lowest frequency 

component for PPM is about 1 /2 Mt, where M is the number of time slots and r 

is the time slot width. This is because the longest time period which the 

voltage baseline held is for a pulse in the first slot of a PPM word and a 
pulse in the last slot of the following word, a total time of approximately 
2Mt . The lowpass filter in Figure 7-8 effectively tracks the average power 
of the data waveform. For cases where the data spectrum is comparable to the 
sample frequency, a special scheme is required to avoid spectral interference 
(Section 7.10). 


TABLE 7.7-1. ASTRONOMICAL BACKGROUND SOURCES (x = 0.83 mn) 


Source 

Maximum 

Size (urad) 

Maximum 

Radiance 

(Wm-2sr-U-l) 

Maximum 

Irradiance 

(Wm-^A-l) 

Sun 

9300 

1800 

1.05 x 10- 1 

Moon 

8960 

4.5 x 10-3 

2.7 x 10- 7 

Venus 

85 

1.8 x 10-2 

1.0 x 10-1° 

Jupiter 

226 

4.5 x 10- 4 

1.8 x 10- 11 

Sirius 

2.8 x lO- 2 

7600 

4.8 x 10-12 

£. Orionus 

~1.1 x lO- 3 

-1.2 x 10 5 

1.1 x 10-13 
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Block Diagram of Tracking Detection Circuit. 
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which gives K $ * 4.1 x 10 4 . Note this Is three orders of magnitude less 
than the noise photocounts, Justifying the neglecting of signal noise for a 
PIN detector. 

For integration time of t « 250 usee (4 kHz sample rate), quantum 
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efficiency n « 0.83, and photon energy hv * 2.34 x 10 joules, the 
required power per element Is P R » (hv/nT ) « 46.2 pW. The total power 
required by the quad array Is then P R « 185 pW. 

Using the link parameters of Table 7.6-1, T^ = 0.1 (10 percent tracking 
split), 50 percent dead zone loss, and Airy lobe beamwldth 9y = 14 M rad, the 
required average transmitter power Is Py = 228 mW. Since the noise was not 
dependent on the tracking power split (background was negllble), by using a 
percent tracking split the average transmitter power can be reduced to 46 mW. 
Since commercially available laser diodes have only 20 mW cw power, even with 
a 50 percent tracking split and a 14 prad beam the link is not reliably closed 
with a PIN quadrant detector. To provide a more implementable and reliable 
link, a quadrant avalanche photodiode detector is evaluated in the following 
section. 

Note that for the PIN quadrant detector tracking analysis, a sample rate 
of 4 kHz was assumed, which Is representative of the complementary filter 
servo tracking approach. A stabilized servo system approach will require a 
much slower sampling rate. For example a gyro-stabilized tracking servo with 
a 1.0 Hz bandwidth, requires only f s 10 Hz sample rate. It was shown In 
Section 7.2 that SNR " Pywhere Py is the transmitted laser power 
(Py “ n^, the signal count rate). Thus to achieve a given SNR, a sample 
rate of f = 4 kHz requires a factor of 20 more laser power than f » 10 
Hz. Although the PIN quadrant detector Is not sensitive enough for a 4 kHz 
sample rate, the PIN quadrant detector is more than adequate for the gyro- 
stabilized servo sample rate. Thus to use a PIN quadrant detector for 
tracking requires adopting the gyro-stabilized servo appraoch, or obtaining 
6 - 9 dB more laser power. 
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Note that the optical background counts, even for a 50 percent tracking 
split, are only Kg ■ 2150 compared to Kg « 4.1. x 10 . As mentioned 
before, the background shot noise was negligible when compared with other 
noise sources. The background DC nonunlformlty Is also negligible, since the 
worst possible case Is 5 percent of the signal. Thus either tracking circuit 
In Figure 7.7-2 or 7.7-4 could be used with a PIN quadrant detector. 

7.8 AVALANCHE PH0T0DI00E QUADRANT DETECTOR TRACKING 

7.8.1 Quad APD Detector 

The only known commercial vendor of quadrant avalanche photodiode detec¬ 
tors at the present time is RCA Electro-Optics In Quebec, Canada. They have 
been involved In single element and quad array avalanche photodiodes (APD) for 
a number of years. Their single element APD is widely recognized as one of 
the best available. Hughes Industrial Products has recently fabricated some 
4 x 10 APD arrays, but these are not yet commercially available. 

Two versions of the quadrant APD have been fabricated by RCA. The old 
version has a common junction, and Isolation between elements Is obtained by 
segmenting. The advantage of this approach Is that there is virtually no dead 
space between the elements (Section 7.3). However this - device is very 
difficult to fabricate, and also difficult to radiation harden. The 
fabrication difficulties cause RCA to have little or no Interest In continued 
production of this device. 

The newer version of the quad APD has 4 separate junctions to form the 

quad element. This results In a dead space between quadrants of 75 pm, but 
there Is an ongoing program to reduce this to 25-30 pm. The entire quad array 
is the same size as their single element APD (1.5 nm), and If the four 

elements are tied together they achieve the same performance as their single 

element APD. A single bias must be used for all four elements. As the 

detector gain Is Increased by Increasing the bias voltage, the gain 
nonuniformity between elements increases. For a gain of 100, the gain 
nonuniformlty is +10 - 15 percent. For a gain of 50, the nonuniformity will 
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be only ^5 percent. In addition, DC levels on the amplifier input to each 
element can be used to compensate for nonuniformity. Thus nonuniformity does 
not seem to be a significant problem for gains of 50 or less. The new version 
of the quad APD will be assumed for performance analysis In this report. 

Since the four elements together give the performance of the single 
element APD, the assumed parameters (per quad element) are: 


Quantum efficiency: 

Capacitance (per element): 

Gain-dependent dark current (per element): 
Gain-independent dark current (per element): 


n = 0.77 
C = 0.5 pf 
I D = 2.5 x 10- 13 A 
I DS = 2.5 x 10" 9 A 


7.8.2 Quad APD Tracking Performance 


Assuming capacitance of 1.5 pf per amplifier, and 0.5 pf per quad 
element, the total capacitance per element is 2.0 pf. The thermal noise then 


6 2 

becomes (Section 7.5) Kj = 3.9 x 10 /G . For integration time f= 

Msec (sample rate 4 kHz), the gain-dependent dark photocount is K n = 391 per 

u & 2 

element. The gain-independent dark photocount is K DS = 3.9 x 10 /G per 
element. 

Assuming a dead space of 75 um, x/D R = 5.6 prad for D R = 15 cm, and 

50 percent dead zone loss, a focal length of 26.8 m is required. If the dead 
zone is reduced to 30 urn, the focal length becomes a more reasonable 10.7 m. 
The detector array field of view is given by e R = D/'f, where D is the quad 
diameter and f is the focal length. Using D = 1.5 mm, d = 30Min, and f = 10.7 
m, e R = 140 prad. Referring to Table 7-2 for background sources, the worst 

case irradiance (neglecting the sun) for 140 u rad is for Venus with L. = 1.0 

-10 -2 *-l 1 

x 10 WM A . The resulting background photocounts are then K„ = 

7600 Tj, where T^ is the fraction of received intensity sent to the 

tracking detector. 

For K„ = 


250 


6 = G C K 


= 391, 

3.9 x 10 6 


10 percent tracking split so that K„ = 760, 
guessing ~ 400, the optimum APD gain is 


G = 63 (Section 7.5). Since the optimum gain function is fairly broad, and to 
reduce gain nonuniformity, we will assume G = 50. If the tracking split is 
increased, so that background noise increases, this will be even closer to the 
true optimum. 
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For G = 50, the excess noise is * 3, and the noise photocounts 
become Kj = 1576, Kg = 391, = 1564, and Kg = 760 (10 percent 

tracking split). 

The total noise count Is 

K n = f N (Kg * K d ) * K os * K t = 6593 

For mean signal to rms noise equal to 5 (Section 6.0) we require 


SNR 




^ P N K S + K N 


5 


and thus require K^ = 445 signal photocounts per element. For n = 0.77 
and t = 250 psec, we require 0.55 pW per element, or total power for the quad 
array P R = 2.2 pW. Note that in the previous section, the PIN quad required 
total power P R = 185 pW, so the quad APD is 19.2 dB more sensitive . 
Assuming Airy lobe beamwidth = 35 prad (we used e- = 14 prad for the 
PIN quad), dead zone loss of 50 percent, 10 percent tracking split, and link 
parameters of Table 7-1, a transmitted power of 17.1 mW is required, which is 
less than the 20 mw cw currently available. Note that the 19.2 dB improvement 
in sensitivity (compared to the PIN quad) has been used to both increase the 
transmitted beamwidth and reduce the required transmitter power. 

For additional tracking link margin with a 20 mW source, we might want to 
increase the fraction of intensity sent to the tracking detector. 

4 

For 25 percent tracking split, K„ = 1900, K„ becomes 1.0 x 10 , and 
for rms SNR = 5, K- = 540 are required. This requires a total power on the 
array of P R = 2.6 pW. For e_ = 35 prad, dead zone loss of 50 percent, and 
25 percent tracking split a transmitted power of P- = 8.2 mW is required. 
This would give a 3.9 dB margin with a 20 mw cw source. 
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If a 50 percent tracking split Is used, Kg = 3800, becomes 1.57 x 
10 , and rms SNR = 5 requires K^ s 665. A received power P R = 3.3 pW Is 
then required on the array. For 6y = 35 urad (Airy lobe), 50 percent 
tracking split, and 50 percent dead zone loss, a transmitted power of Py = 
5.1 mW Is required. This provides 5.9 dB power margin for the tracking link. 

Note that, unlike the PIN quandrant detector case, the optical background 
counts are larger than the signal counts for the APD quadrant detector with 
SNR = 5. Although the K~ are adequate to overcome the shot noise due to the 
background DC current, the mean signal is not sufficient to overcome the DC 
background current Itself. To overcome the DC nonuniformity due to optical 
background, it is recommended that the AC - coupled tracking approach (Figure 
7.7-4) be used with a quadrant APD. Otherwise, much of sensitivity advantage 
over a PIN would be lost. 

In summary, the vast sensitivity improvement of the quad APD over the PIN 
quad makes the tracking link viable with a single 20 mW source. Since the 
beacon link at 10 Kbps will be seen to be compatible with a 20 mW source and 
50 percent tracking split (Section 7.9), up to 5.9 dB of margin can be 
allocated to the tracking link. However, just as in the PIN quad case, a 
program is needed to reduce the dead zone to 30 urn or less, to keep the 
required focal lengths reasonable. 

e 

7.9 BEACON LINK - TRACKING SPLIT 

Although a low data rate beacon link from GEO to LEO is assumed in the 
SOW for this report, no specific data rate is mentioned. As an example of the 
typical telemetry beacon link a data rate of 10 Kbps has been assumed, and 
discussion with NASA representatives has confirmed this is reasonable. 

The beacon link communication format will be taken as pulse position 
modulation (PPM), since it has become recognized as the most efficient 
modulation scheme for energy detecting optical communications. The PPM format 
is illustrated in Figure 7.9-1. 
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Since the peak to average power ratio of commercially available 20 mW 
laser diodes has been experimentally measured at Hughes to be In the range of 
8-16, we will choose the number of time slots to be M = 8. Hughes experiments 
have shown that commercial 20 mW diodes (RCA, Telefunken, Amperex) can 
maintain their 20 mw average power at least up to M = 8, for a peak power of 
160 mW. This may be true only for pulses shorter than 1.0 nsec, however. 

Since log^M bits are transmitted per pulse, and a pulse occurs every 
Mt seconds, where t is the time slot width, the PPM data rate is 

R q > log 2 M/Mx 

For R D = 10 Kbps, and M = 8, we have x = 37.5 usee. 

We will assume a single element APO such as the RCA C30902E, with a total 

detector-preamp capacitance of 2.0 pf. The thermal noise is then K T = 3.9 x 

6 2 ' 

10 /G , where G is the mean detector gain. For gain-dependent dark 

-12 

current Ig = 10 A, = I D T /e = 234. For gain-independent dark 

current I Q - = 10 nA, = 2.3 x 10®/G 2 . 

For a detector field of view e D = 40 yrad, the worst case background 

K 2 

(except the sun) is Venus (Table 7-2), with radiance L R = 1.8 x 10 
WM' 2 Sr' 1 A' 1 . Using the link parameters of Table 7-1, n = 0.77,T= 37.5 
Msec, and assuming a 50 percent tracking split, the background photocounts are 
Kg = 475. (A 50 percent tracking split is assumed to give the worst case 
required transmitter power.) 

For Kg = 475, K Q = 234, 6 = G 2 K T = 3.9 x 10 6 , and using 

K $ 450 yields an optimum gain around G = 70. Since the gain optimum is 
broad, and tabulated results exist for G = 100 (F^ = 4) we will analyze the 
system for G = 100. (It should be recognized that G = 70 should be used for 
actual system implementation.) 

For G = 100, the noise photocounts are K- = 393, K g = 475, K Q = 234, 
and Kg<j = 234, so that the total noise counts are K^ = F^fKg + K 0^ + 
Kg S + Kj = 3500. For error probability ' PE = 10"®, Kj = 470 are 

required, and P R = 3.8 pW is required at the detector. Still assuming a 50 
percent (worst case) tracking split, and e- = 3b prad (Airy lobe), the 
required peak transmitter power is P. = 3.0 mW. If the transmitter peak 
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power Is run at the cw level of 20 mW, the beacon link has 8.2 dB power 
margin. If the 20 mW average power source has 160 mW peak power at M = 8, as 
experimentally measured by Hughes, the peak power margin becomes 17.3 dB. For 
the 37.5 ysec pulses, the laser may be constrained to 20 mW peak with 8.2 dB 
margin. 

In summary, from a power standpoint the beacon link presents no design 
problems, even with 50 percent of the transmitted power allocated to the 
tracking system. 

7.10 TRACKING ON BEACON PPM 

As shown in the previous section, from a power standpoint there is no 
problem sharing the GEO beacon for both tracking and telemetry data at the 
LEO. However, there could be a separate problem due to interference of the 
beacon data modulation with the tracking signal, since for a tracking sample 
rate of 4 kHz and beacon data rate of 10 Kbps there is significant spectral 
overlap. In other words, if the tracking system samples at 4 kHz, sometimes 
it samples while a pulse is present, and sometimes when the pulse is absent. 
A 4 kHz low pass filter on the 10 kbps modulation will not remove all the 
signal fluctuation from the tracking system. Ironically, if the data rate 
were 100 Mbps or more, the tracking system could utilize the average power of 
the data signal, and would not respond to the high frequency modulation. One 
approach to avoiding this problem is to use a single quad APD for both 
tracking and beacon data modulation, rather than power splitting between a 
quad detector for tracking and another detector for beacon data. This can be 
accomplished by summing the output from the four quad elements to form the 
beacon data output, while the quad outputs are simultaneously summed and 
differenced in the usual way to provide the tracking error signals. A circuit 
for tracking with PPM would first examine the sum photocount of all four 
quadrants in each time slot, while storing the individual photocount from each 
quadrant in each time slot (4M storage registers). At the end of the PPM word 
the time slot with the maximum sum photocount is chosen as the slot containing 
the optical pulse, which demodulates the beacon data. Then the quadrant 
outputs from the time slot which has been determined to contain the pulse are 
summed and differenced to provide the tracking error signals. Note that the 
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circuit must first determine which time slot contained the pulse before the 
tracking signal can be generated, and the tracking error signals are updated 
at the PPM word rate. If the tracking sample rate is slower than the PPM word 
rate, the tracking error signals could be loaded Into a buffer register at the 
PPM word rate, and the buffer could then be read out at a slower rate. The 
PPM word rate should be a multiple of the tracking sample rate, so that 
readout does not occur while the buffer register Is being loaded. 

As an example, a system design for a 4 kHz tracking sample rate will be 
generated. The PPM word rate will be chosen to match the 4 kHz sample rate, 
so that 1/Mt = 4 kHz. 

If we then choose M & 8, the time slot width is t = 31.25 nsec. Since 
the pulses occur at the word rate of 4 kHz, and for M = 8 three bits are 
transmitted per pulse, the beacon data rate is 12 Kbps. 

We will assume a quad APD, such as analyzed in section 7.8. The noise 

photocounts from section 7.8 are Ky = 1576, Kg = 391, K DS s 1564, and 

where T^ was the tracking power split. Assuming 
T$ = 1.0, since there is no power split in this approach, and scaling the 

noise counts for = 31.25 p$ec instead of 250 psec (except K, - section 

7.5.2), the noise counts per element become Ky = 1576, Kg = 49, Kg<. = 

196, and Kg = 2625. The total noise counts are 

To obtain an rms SNR = 5 as needed- for tracking, 
K- = 534 are required. Thus 5.2 pw are required per element, for a total 

received power of P. = 20.8 pw. Assuming 50 percent dead zone loss, 

Oy = 35 prad (Airy lobe), and link parameters as in Table 7-1, the required 
peak transmitted power is Py = 4.0 mw. Note this ^s even less required 

transmitter power than for the 50 percent tracking split case in section 7.8, 
which required 5.1 mW. However, we still must verify that this is also enough 
power to demodulate the beacon data. 

To determine the communication performance with the summed quad elements 
as the communication detector, we use the fact that the sum of independent 
Gaussian random variables is Gaussian, with mean given by the sum of the means 
and variance given by the sum of the variances. Therefore 


K B * 2.1 x 10 T $ , 


k n = wv 


k ds +k t s 9794 ‘ 
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m~ * 4 eGKj 

•s * 4 A2 <Vs * k n' 

Note that previously derived equations for error probability can be used by 

I I 

considering in the equation to be 4K$, and in the equation 

as 4K n< The noise counts per element per integration time are K N = 9794, 
as computed for tracking. The error probability (PE) is then based on 

= 4 K n = 3.92 x 10 4 . For PE = 10" 6 , K* = 1500 are 

required. Since the required = 1/4 «<. = 375 per element are less 

than the = 534 required for tracking, Py = 4.0 mw is sufficient peak 

transmitter power to perform both the tracking and communication functions. 

Note that since this sytem tracks a data waveform, each quadrant 

dectector output can be AC-coupled. This avoids the problems of DC 

non-uniformity due to optical background. 

Thus using a single quad APD detector to perform both tracking and beacon 
data detection is an appealing alternative to power splitting with separate 

detectors. This technique assures that the beacon data does not spectrally 

interfere with the tracking signal. 

7.11 TRACKING WITH AN ACQUISITION ARRAY 

If an array is used for acquisition, which will be seen in section 8 to 
have significant advantages quad cell acquisition, this array could also 
perform the tracking function. The advantage of tracking with an array is to 
combine two separate acquisition and tracking subsystems into a single 
subsystem. A separate tracking detector, it optics and electronics are 
eliminated, as well as a bifurcated mirror and electronics to perform the 
handover function. 
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A video tracker examines the video output of the array, and generates the 
angular error signals to the tracking servo. The array and video tracker must 
provide angular error samples at a rate compatible with the tracking servo (4 
KHz for complementary filter, 10 Hz for gyro-stabilized). 

Array tracking Is made possible by using a video tracker, which Is an 
existing Hughes product line. The video tracker digitizes the video ouptut of 
each pixel within a tracking gate, and tracks the desired Image In either 
centroid or correlation mode. In correlation tracking, the digitized Image 
from the previous frame Is correlated with the Image within the current frame, 
to deduce the apparent image motion. Centroid tracking establishes a video 
threshold, and tracks on the centroid of the Image intensity that exceeds 

threshold. For tracking a distant optical communication terminal, centroid 
tracking would be used, with the threshold set to eliminate unwanted back¬ 
ground, but allowing the laser to appear to a bright spot. The tracking 

system would then align the line of sight with the center of the laser spot. 

Existing CCD array and video trackers operate at frame rates up to 60 Hz 
(interlaced), which is compatible with the requirements of the gyro-stabilized 
servo. If a CID array were used instead of a CCD array (Section 8.2*1), the 
x-y addressability of the CID array could be used to define a tracking gate. 
This tracking gate would consist of a much smaller subgroup of pixels which 
could then be read out at a much more rapid frame rate. The video tracker 
already incorporates a similar tracking gate in Its processing, so that its 
tracking algorithm does not have to be applied to the entire video field. A 

5 

400 x 500 array reads out 2 x 10 pixels In 1/30 of a second. At this rate, 

1500 pixels could be read out in 250 u$ec (4 KHz rate). Thus a 38 x 38 pixel 

track gate could be implemented, which is adequate since the distant receiver 
is unresolved within a pixel. The CID sacrifices about 3 - 5 dB in low light 
sensitivity compared to the CCD, but to eliminate star background the CCD 
performance will be background-limited anyway. 

For a 400 x 500 array, staring into a 0.5* (8.7 mrad) uncertainty cone 
(see Section 8.0), the pixel field-of-view is 8.7 x 10'^/400 = 22 urad. 






Setting the optical blur circle over 2x2 pixels gives the angular blur 
subtense of 44 urad. For reasonably high SNR, the tracking jitter Is closely 
approximated by 

o - 0.685 a/(SNR-3) 

where the a Is the angular suDcense of the optical blur, and SNR Is the ratio 
of mean signal to rms noise. For a = 44 urad, an SNR = 33 Is required to 
achieve o = 1.0 urad rms jitter. If a 1000 x 1000 array were used, with a = 
18 urad, an SNR = 15 is required. (If a complementary filter tracking servo 
is used, the sensor jitter contribution must be less than 1.0 urad.) 

Suppose a 400 x 500 CCD array is used with a gyro-stabilized servo at a 
10 Hz frame rate. The DC background due to thermal noise and dark current is 
fixed, and the other levels ride on top of these, so setting a threshold 
eliminates them. The DC optical background may or may not be present in a’ 
pixel, so that if the optical background Is eliminated with a threshold, the 
signal must be able to exceed this threshold and the shot noise variance even 
with the DC optical background absent. With these considerations, the SNR can 
be defined as 


(mean signal - DC optical background) 

SNR = - 

rms noise 

(This assumes that the array can only output the accumulated charge per pixel, 
rather than a modulated signal.) 

For integration times of r = 100 msec (10 Hz frame rate) the dark current 
nonuniformity dominates the noise variance. (The assumed parameters here are 
discussed In Section 8.5.3). For dark nonuniformity of 1.0 nA/cnr (15 C), 
12 um x 18 um pixel area, and r = 100 msec Integration time, the rms dark non- 
uniformity count Is K qn = 1350. 

For pixel F0V = 22 urad and Venus background, the optical background DC 
level Is Kg = 2.5 x 10^. SNR = 33 Is required for 1.0 urad gyro-stabil¬ 
ized jitter performance. Thus 
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or Kj = 3.0 x 10 . For 25 percent quantum efficiency, this Is a required 
power of 27.0 pW. Using the link parameters of Table 7.6-1, with Oj = 35 
Mrad (Airy lobe) and T- *• 0.5 (50 percent tracking split), an average power 
of 21.0 mW is required. (A 12 Kbps beacon link was shown to be compatible 
with a 50 percent tracking split with an APD detector). Therefore tracking 
using the acquisition array should be considered as a viable system option. 

7.12 SUMMARY - TRACKING AND BEACON LINK 

Section 7.0 has examined the problem of tracking a data beacon from the 
GEO satellite to less than one microradian rms accuracy in order to derive 
pointing commands for a high data rate LEO transmitter. The noise equivalent 
angle of a quadrant tracking detector was derived explicitly stating the 
assumed signal and noise parameters per element. The tracking detector 
quadrant gap has been identified as an area needing further technology 
development, but which does not preclude operational systems with currently 
available devices. The angular dynamic range over which closed loop tracking 
can be maintained was shown to be proportional to the rms SNR in an element of 
the quadrant detector. 

Photodetection analysis and link equations were presented to allow 
calculation of tracking (and acquisition) link budgets. It was shown that a 
PIN quadrant detector was not sensitive enough to close the tracking link 
budget to 1.0 urad rms accuracy unless a gyro-stabilized tracking servo is 
used or laser of 200 mW or greater are available. The quadrant APD detector 
is about 19 dB more sensitive than the PIN, which allows 1.0 yrad tracking 
accuracy with the complementary filter servo approach, 35 yrad GEO beamwidth, 
a commercial 20 mW laser diode, and 6 dB link margin. 
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For tracking a GEO 10 Kbps beacon data link three different configura¬ 
tions were shown to be viable. A 50-50 optical beamspllt between an APD 

quadrant tracking detector and single element APD data detector satisfied both 
functions (with adequate margin) with a conmerclal 20 mW laser. Or the 
quadrant APD tracking detector may also provide the data detection by summing 
the four outputs. This may be necessary If the tracking sample rate has 

significant spectral overlap with the beacon data, as would be the case with a 
4.0 KHz complementary filter sample rate and 10 Kbps beacon data. Finally, It 
was shown that If an array were selected as the acquisition detector. It could 
also perform the tracking function. A CCD (or CID) array with standard video 
frame rates is compatible with the gyro-stabilized tracking servo. A CID 
array with a circuit to create a 40 x 40 pixel tracking gate could provide 
updates at the 4.0 KHz complementary filter servo sample rate. 
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8.0 ACQUISITION 


8.1 INTRODUCTION 

Acquisition for optical communication systems Is the process by which two 
terminals Initially locate each other and progress to mutual fine tracking and 
two way communication. Each terminal Initially knows the location of the 
other within some angle of uncertainty, which must be reduced to the order of 
microradians In order to fine track with narrow communication beams over 
Intersatellite distances. This may involve handover between sensors within 
the acquisition subsystem before eventual handover to the fine tracking 
system. In one extreme the acquisition sources and sensors may be wholly 
separate to the acquisition subsystem, or In the other extreme the 
communication sensors and sources could be utilized for acquisition. 

We have performed acquisition analysis with several driving goals: 

a) minimum number of components In addition to those for communication 

b) reliable acquisition times less than 5-10 minutes 

0 

c) single stage for reduction of Initial uncertainty to fine tracking 

d) avoid mechanical placement or removal of optical elements in the 
optical train. 

Goal (a) is motivated by keeping the system as simple as possible. Goal (b) 
results since the low earth orbiting (LEO) satellite in our study orbits the 
earth every two hours. The acquisition time should be small compared to the 
roughly one hours conmunication time available per orbit. Since the earth 
blocks the LEO as seen from geosynchronous (GEO) orbit, acquisition must be 
performed every LEO orbit, or about twelve times a day. Thus the acquisition 
sequence should be very reliable, since It must be performed frequently. Goal 
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(c) again reflects the desire for system simplicity, to minimize sensor 
handover sequences. Goal (d) Is motivated by the difficulty of maintaining 
microradian optical alignment when mechanically adding or removing lenses or 
beamsplitters. Thus we have resisted using a lens to spoil the communication 
beam to serve as an acquisition beacon, but Instead recommend that a separate 
beacon be used. 

The acquisition subsystem has received the least attention In past 
analyses of optical communication systems. Since no conmunlcatlon can occur 
until acquisition takes place, the overall system availability Is ultimately 
dependent on acquisition performance. 

The determination of the acquisition uncertainty angle Is a critical 

first step In designing an acquisition subsystem. Since the acquisition time 

2 

will be shown to be proportional to «y N( . ^ ®UNC the acc l u ^ s ^ t ^ on 

uncertainty angle), and since background radiation photocounts are 


proportional to e 


UNC 


It is desirable to keep e 


UNC 


as small as 


possible. In general, however, this Is not an Independent design variable, 
but depends on the specific spacecraft and Its control systems. 

In general, the acquisition uncertainty Is dependent on the spacecraft's 
own attitude uncertainty, the orbital ephemerls uncertainties of the position 
and velocity of the other terminal, relative timing errors, etc. However, 
using the specific scenario of Landsat for the LEO and TORS for the GEO, we 
found that the attitude uncertainties were significant fractions of a degree, 
while the other errors were of the order of tenths of milliradians, so that 
attitude uncertainty clearly dominates. In fact, the TORS yaw axis attitude 
uncertainty at the ground station is the single dominant influence to the 
acquisition uncertainty angle. 

Based on conversations with Hup Oamare (NASA GSFC), simulations reveal 
maximum attitude uncertainties of 0.08 degrees In the roll and pitch axes, and 
0.2 degrees in the yaw axis. These are the maximum magnitudes of error before 
the next sensor update. The yaw axis sensor update rate is much slower than 
the other axes, resulting In Its larger attitude uncertainty. Glmbal servo 
analysis reveals that the yaw axis uncertainty may be fully coupled Into the 
line of sight uncertainty. 
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Assuming the uncertainties of the three axes are Independent, the root 
square sum has magnitude 0.23 degrees, i.e. ^ 0.23 degrees. To be slightly 
conservative, we have assumed that the uncertainty region is a 0.5 degree 
diameter cone, and our acquisition analysis is based on this value. 










8.2 ACQUISITION SENSORS 


8.2.1 Arrays 

Sensors for the acquisition subsystem fall into two broad categories - 
quad cells or arrays. Quad cells will be discussed In the next section. 
Arrays can be charge - coupled devices (CCD), charge - Injection devices 
(CIO), or MOS photodiode arrays which use charge - coupled readout 
technology. Although a linear array could be scanned across the uncertainty 
solid angle, we have chosen to avoid mechanical scanning and focus on area 
arrays. 

The area arrays function as any energy detector, since the photocurrent 
resulting from incident light is integrated to yield a collected charge within 
each pixel. The various arrays differ in how this charge is stored and 
ultimately read out. The CCD readout is similar to an analog shift register, 
in which the charge is transferred from pixel site to pixel site. Voltages 
which control potential wells are manipulated at the proper time sequences to 
shift the charge along the row. The endmost vertical column of charge falls 
into a register which is then serially read out before the next column of 
charge is clocked in. Figure 8.2-1 schematically depicts the CCD readout 
process. Early CCD devices were limited by charge loss during charge 
transfer, such that a point of light could become a horizontal streak. 
Current devices, however, have transfer efficiencies of 0.9999 - 0.99999, and 
this is no longer the low light limitation. The CCD readout is serial, and 
timing information is retained to restore the pixel to its proper location in 
a video display. Although the CCD readout technique requires readout of the 
entire array to address any specific pixel, there are two distinct advantages 
to optical sensitivity which results. 

First, since the image sites collect and transfer the charge, no separate 
charge storage areas are needed. The CCD focal plane filling efficiency is 
virtually 100 percent, so that no optical Intensity is wasted by illuminating 
non-responsive charge storage sites. Second, since the charge is read out 








IALLEL IMAGING REGISTERS 


I ON-CHIP 

AMPLIFIER 


Schematic diagram of CCD readout technique. 


ichematlc Diagram of CCD Readout Technique. 




.. 





* A 



















along one-dimension, electrical connections and associated capacitance are 
minimized. The sensitivity is directly reduced as the capacitance is 
increased. So the combination of low capacitance and 100 percent focal plane 
light responsivity cause the CCD's to be the most sensitive area array 
detectors, at the expense of x-y pixel addressability. 

The charge-injection device (CID) array uses an x-y coincident voltage 
readout technique, which allows a particular pixel location to be addressed. 
The two-dimensional nature of the readout results in more electrical 
connections and higher capacitance than for a CCD, so the CID sensitivity is 
poorer than the CCD. In parallel injection CID's It is possible to perform 
non-destructive readout, i.e. leaving the charge in the pixels by deferring 
the injection operation. Otherwise, injection is used to empty the charge 
storage capacitors after line readout. In sequential injection CID's the 
photon-generated charge is injected into the substrate and the resultant 
displacement current is sensed. Figure 8.2-2 shows schematically how CID 
readout is accomplished. Each sensing site consists of two MOS capacitors 
which are coupled to facilitate charge transfer. A row is selected for 
readout by setting its vertical scan register. Signal charge at all sites in 
that row is transferred to the column capacitors. The charge is then injected 
into the substrate by driving each column voltage to zero, in sequence, using 
the horizontal scan register and signal line. Charge in' unselected rows 
remains stored in the row capacitors during the column voltage pulses. Thus 
any row can be selected and read out, but having both vertical and horizontal 
voltage signal lines connected to each sensing site is the reason CID capacit¬ 
ance is greater than for a CCD. The x-y coincident voltage configuration of 
the CID is also adaptable to random pixel addressing, in which only the 
selected pixel receives coincident x-y voltages and injects its charge into 
the substrate. This is not seen to be a significant advantage over CCD's for 
acquisition, since frame rates are already fast (10 msec or less), and a 
relatively small percentage of aqcquisition time is expended after array 
illumination. However, if the acquisition array is used for tracking also, 
the CID pixel addressability allows a pixel tracking gate to be formed with a 
small group of pixels. These may be then read out at a much higher frame rate 
as dictated by tracking requirements (Section 7.11). 
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The MOS photodiode Image sensor Is an array of photodiodes with an MOS 
transistor scanning circuit embedded into the same monolithic structure. The 
photodiode array loses some responslvlty due to added capacitance from the x-y 
coincident voltage readout as In the CID. The photodiode array sacrifices 
additional responslvlty due to sharing the focal plane between the MOS 
scanning transistors, a loss of up to 40 percent. The MOS photodiode array 
was not Initially thought to allow large video arrays, due to difficult 
fabrication of Its complicated focal plane structure, and because the 

Individual cell size can require up to four times the area of a CCD cell. 

Recently, however, photodiode arrays of about 400 x 500 have become 

commercially available (Figure 8.2-3), which Is comparable to current CCD 
arrays. 

A summary of commercially available imaging arrays is shown in Table 
8.2-1. Array sizes of up to about 400 x 500 are available in both CCD and 
photodiode technology, and up to about 300 x 400 for CID. 

Unfortunately, most manufacturers specify the array sensitivity in 

luminous units, which takes into account the spectral responsivity of the 
human eye. However, a good Idea of the relative sensitivity of the CCD, CID, 
and photodiode array formats at standard video rates of 30 frames per second 
can be obtained from the first four entries in Table 8.2-1. The minimum 
illumination for the Fairchild and Sony CCD's are 5 lux and 3 lux, respec¬ 
tively. The GE CID requires 11 lux, a factor of 2-3 times more illumination 
than the CCD. The Hitachi photodiode array requires 30 lux, about a factor of 
three more than the CID, and a factor of ten more than the CCD. The relative 
radiometric sensitivity follows the same ratios, since less than 1.0 pw 
sensitivity for CCD and 10.0 pW sensitivity for photodiode arrays has been 
measured at Hughes. 

Fairchild and RCA both have intensified SIT-CCD arrays available. The 
Fairchild array has a microchannel plate illuminating a phosphor tube, which 
is then coupled by fiber optics into the CCD. The RCA version is more 
sensitive, since the intensifier tube launches directly into the CCD, 
resulting in secondary electron "avalanche gain" of about 1500. The RCA 
SIT-CCD has the same sensitivity as the well known SIT vidicon, which is 
200 - 300 times more sensitive than the unintensified CCD. However, the 
10 - 15 kV voltage requirement of the intensifier tube may make it undesirable 
from some space applications. 
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TABLE 8.2-1. COMMERCIALLY AVAILABLE IMAGING ARRAYS 


Manufacturer 

_ 

Type 

No. of Pixels 

Minimum 

Illumination 

FAIRCHILD 

CCD 

(SIT-CCO) 

380 x 488 

5 LUX @ F 1.4 

SONY 

CCD 

384 x 491 

3 LUX @ F 1.4 

GE 

CID 

290 x 416 

11 LUX 

HITACHI 

MOS P.D. 

384 x 485 

30 LUX @ F 1.4 

RCA 

SIT CCD 
(CCD) 

403 x 512 

1.1 x 10- 2 LUX 

EG &G RETICON 

CCPD 

256 x 256 

IQ -3 uJ/cm2 

TEXAS INSTR. 

CCD 

292 x 290 
(1024 x 1024 for 
space telescope) 

• 

- 2.0 x 10- 3 M W/cm 2 

HUGHES IPO 

CCD 

(radiation- 

hardened) 

1 

100 x 100 
(324 x 324) 

1 

< 1.0 pW 



















Texas Instruments Central Research Laboratory Is building two large 
space-qualified CCD arrays for space programs at the Jet Propulsion 
Laboratory. An 800 x 800 CCD array Is being fabricated for use on the Galileo 
program, and a 1024 x 1024 CCD array Is being built for Space Telescope. The 
space qualification Included the packaging, documentation, and vibration 
testing aspects, as well as radiation shielding for the Galileo array. Both 
the arrays read out at frame rates of 1-2 frames per second. 

In summary, solid state arrays are commercially available In sizes up to 
400 x 500, with sub-plcowatt sensitivity for the CCD's. Since CCD arrays with 
standard video rates will be shown to be compatible with acquisition sub¬ 
systems for satellite optical communication (Section 8.5), experimental 
testing by NASA personnel Is highly recommended. 

8.2.2 Quadrant Detectors 

A summary of commercially available quadrant detectors is shown in Table 
8.2-2, and the Hughes HPIN quadrant detector data sheet is shown in more 
detail in Figure 7.7-1. The Hughes quad detectors seem to offer the smallest 
dark current and capacitance for a given ratio of array diameter to dead zone 
width, which is important for field of view reduction during acquisition. The 
RCA quadrant APD is listed with a dead zone of 2.9 mils (75 urn), but a program 
is ongoing to reduce that value to 1.‘2 mils (25 pm). Unless otherwise 
specified, the Hughes quadrant detectors will be assumed for the PIN quad 
case, and the RCA for the APD quad case. 
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8.3 ACQUISITION SOURCES 




Two basic types of semiconductor laser sources for acquisition are 
presently commercially available: fundamental transverse mode single laser 
sources, and multi-transverse mode multi-laser Incoherent arrays. Commerci¬ 
ally available single transverse model lasers are summarized in Table 8.3-1. 
Average power In the 20 - 30 mW range Is presently available, and device power 
Improvement has been rapid. Several manufacturers have now reported 100 mW In 
the laboratory. Probably devices with 50 mW cw will be commercially available 
In the 1986 - 88 time frame. Perhaps by that time frame some version of the 
Xerox multiple stripe coherent array will also be available. Although 1.0 
watt cw has been reported in the laboratory for the Xerox array, lifetimes 
have been only several hundred hours. It Is difficult to extrapolate the 
power a reliable coherent array might obtain. In addition. Incoherent power 
combining schemes with diffraction gratings or dlchrolc mirrors might be used, 
but practical considerations may prohibit using a power combiner for both the 
communication and acquisition sources. 

The second type of acquisition source is the multiple-laser incoherent 
array, in which several multi-transverse mode lasers are combined in the image 
plane of the optical system. Laser Diode Labs (LDL) Is the industry leader in 
this technology, and high peak powers are available. Some devices physically 
position the lasers as close together as possible, and others use fiber optics 
to couple the emission from all the lasers into the Image plane. Figure 8.3-1 
summarizes the characteristics of the LDL LS-410 array. Although 1.0 kW peak 
peak power Initially seems impressive, a practical problem arises in collimat¬ 
ing the array. 

The basic problem with the LDL arrays is one of collimating an extended 
multiple-laser source with the broad angular emission characteristic of a 
single laser. The efficiency with which the optical system captures the 
emitted light is given by 
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• UP TO 1.2KW PEAK OPTICAL SOURCE POWER. 

• UP TO 0.6 AVERAGE SOURCE POWER OUTPUT. 

• 1 TO 5 KHZ REPETITION RATE. 

• 115 VAC. 60/50 HZ OPERATION. 

• 650 OR 904 NM STANDARD PEAK EMISSION 
WAVELENGTH. 

• CUSTOM UNITS AVAILABLE INCLUDING 800 
TO 880 NM WAVELENGTH SELECTION, 28 VDC 
OR 220 VOLTS/50 HZ INPUT POWER. 


DESCRIPTION 

The LS-410 and LSA-410 are intended as high pulsed power laser 
illuminators. The unit is capable of emitting over 750 warts peak 
optical power into the air with an 8 milliradian beam divergence at a 
5kHz repetition rate Up to 0.5 watts average optical power is 
available. The unit is complete, providing the laser optics, pulse 
circuitry control and f>ower supply in rugged and functional housings. 
The pulses may be tnggered from an internal or external source. 
Custom designs are available on special order. 



Laser Diode Laboratones high power illuminators are a line 
operated completely se'f contained system ideal for applications 
in Infrared viewing systems. Infrared target designators, etc. 

Figure 8.3-la). Laser Diode Labs High Pulsed Power Array. 
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KRFOfcMAKCt CKA.F’ACURU.; 


Peak power output . 

Pea* wavelength emission 
Spectral width 

(50% points). 

Average pcwer output at 
5KHZ (approx.). 


1000 watts ... 
904nm. 

3 5nm. 

0 6 watts . 

170x190 mils 


1000 watts 
830 nm 

3 5nm 

0 4 watts 
200x190 mils 


• • ; 

• • • 

Peak current. 40 amps . Up to 40 amps 

Optical pulse duration . 130ns . 90ns max 

Repetition rate (internal 

or external trigger) . 1 — 5KH2 . 1 — 5KHZ 


Peak power output 
(w»th Fl 8 lens) . 
Beam divergence . 
Size . 


* 


750 wans . 750 watts 

8mr. 8mr 

Laser transmitter . * 9 5' OD x 20 2 

long 

Power supply. * 17" x 16 x 8 

Total weight approx . * 30 lbs 


Figure 8.3-lb). Laser Diode Labs High Pulsed Power Array Data Sheet. 
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T NA « If exp (-x 2 / 2 o x 2 ) exp (-y 2 / 2 o y 2 ) dx dy 
” A ~ B 


exp (-x 2 / 2 o x 2 ) exp (-y 2 / 2 o y 2 ) dx dy 

•00 -00 

= erf (A//To x ) erf (B /-/To ) 



where erf (x) Is the error function, A and B are the half angles of acceptance 
of the numerical aperture of the optical system, and a and o are the 

p * j 

l/e half angles of the laser Gaussian beam divergence. The optical system 
must be fast (low f number) to capture a significant fraction of the widely 
divergent light. To capture 60 percent of the emitted light, F 2.5 op'.ics are 
required (experimentally measured). The far field beam divergence is then 
given by 


e FF = V f 

where D- Is the source diameter in the image plane and f is the lens focal 
length. To collimate the beam to the 0.5* (8.7 mrad) acquisition uncertainty, 
wi h F 2.5 optics for 60 percent collection efficiency and 200 mils array 
source diameter, a nine inch diameter lens is required. Even if a more 
complicated collimating optical system is used, the output aperture diameter 
must be nine inches or greater. Assuming only a six inch transmitting 
diameter is available on the GEO, and keeping F 2.5 optics, divergence angle 
is increased a factor of 1.5. The 600 watts peak that was col’ected delivers 
5.0 pW peak to a 15 cm collecting aperture on the LEO. However, this is an 

_3 

average power of only 2.3 x 10 pW, since the maximum pulse width is 90 ns 
and the maximum pulse rate Is 5 kHz. The average power Is too low to reason¬ 
ably detect, so a viable detection scheme would have to time gate around the 
optical pulse, thereby utilizing the laser peak power. This will be discussed 
in Section 8.4.6. 
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8.4 QUADRANT DETECTOR ACQUISITION 


8.4.1 Field of View Reduction 


There are three fields of view associated with a quadrant tracking 
detector. The field of view (FOV) over which a source can cause a spot to 
fall on the active area of at least one quadrant is normally considered to be 
the quadrant detector FOV. This Is given by « D/f, where D is the 

diameter of the quad array, and f is the focal length of the optical system. 
As shown In Section 7.4, the closed loop field of view is the FOV over which a 
source can cause a spot to overlap all four quadrants, such that both azimuth 
and elevation error signals can be obtained. The closed loop field of view is 
given by = D^/f, w ^ ere s P ot diameter in the quadrant 

detector focal plane. The third field of view for a quadrant detector is the 
effective FOV determined by the residual tracking jitter. For a quadrant 
detector used for acquisition, the residual tracking error is associated with 
the quadrant detector noise equivalent angle (NEA), which was discussed in 
Section 7.2. 

Analogous to defining o CL /© NEA = ^ °NEA^ as the c ^ osed loop 

tracking dynamic range in Section 7.4, we can define Opoy/e N r A as the 
dynamic range of FOV reduction for a quadrant detector used, for acquisition. 
Since from Section 7.4 we have 

»NEA * <•«> <V f > ' 

the acquisition FOV reduction ratio (R) is given by 

R s e F0V / ®NEA = (8/ir) (D/I V SNR rms 

Thus the obtainable FOV reduction is proportional to the rms SNR in a "wad 
element, and also proportional to the ratio of quadrant array diameter' (D) to 
optical spot diameter (D<.). 
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The obtainable FOV reduction compared with the ratio of uncertainty angle 
to tracking detector field of view determines tne number of acquisition stages 
required to handover to fine tracking. If the quadrant acquisition array Is 
configured so that e FOV = e UNC , where e UNC Is the angular uncertainty In 
location of the opposite terminal, then the following stage must have a field 
of view greater than or equal to o ( ,r» of the preceding stage. The FOV 
reduction ratio (R) then determines the field of view ratio between successive 
stages 

*®F0V-l /9 F0Y-2*» Slnce 0W EA-1 = ®F0V-2’ 

We have attempted to use single stage acquisition to keep the subsystem 


as simple as possible. This suggests that the single stage acquisition sensor 
must have FOV reduction capability such that R >_ (e^ ^ •mu T ^’ w * iere 


FOV-T 


e 


UNC is the acquisition angular uncertainty and Opg V t is the field of 


view of the tracking detector. Although it would be desirable to handover to 
the closed loop field of view of the tracking detector (e CL ^), this is not 
necessary since a multi-step angular reduction can move the spot from the 


quadrant array field of view 
(Ocl_t^ (Section 8.4.2). 


(e 


FOV 


,) to the closed loop tracking 


FOV 


8.4.2 Servo Angular Step Reduction 


When a quadrant detector is used for acquisition, one approach would be 
to set the closed loop field of view equal to the acquisition angular 
uncertainty. The FOV reduction would then be given by the expression for 
tracking d/namic range (Section 7.4). 



I e NEA = (8/ir) SNR pms 


From Section 8.4.1, if the total quadrant array field of view (e F gy) is set 
equal to the acquisition angular uncertainty, the FOV reduction is given by 


R = ®F0V f °NEA = (D/C V SNR rms * 
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The FOV reduction is increased by the ratio of the quad array diameter (D) to 
the optical' spot diameter (D^). For single stage acquisition the extra 

factor of D/D s is required, and we set e UNf * ®F0V 

Then it becomes necessary to adopt an algorithm to transfer the optical 
spot from being within a single quadrant in Op„y to being within e^, 

where the spot overlaps all four quadrants. This is the same problem that 
occurs when the acquisition sensor hands over to a single quadrant in the 

tracking quad array. The servo system knows only that the spot is within one 
quadrant, but not how far in angle the spot is from Leing centered. A linear 
angular step could be used, in which the servo steps x microradians toward the 
illuminated quadrant, checks to see if an error signal exists, then steps x 
microradians again, until eventually the spot falls with and an error 

signal occurs. 

An angular reduction scheme which is faster than the linear step 

reduction is to use a log step reduction. The first step will be half of the 
quadrant FOV,'the second step one-fourth of the quadrant FOV, etc., until the 
spot falls within e- . If the spot changes quadrants, then the next step is 
in the opposite direction, but still half of the previous angular step. 

We need to provide an estimate of the number of steps required for the 
log angular reduction to transfer the spot from Opgy to o^.. Define the 
ratio of the quad array diameter (0) to the dead zone (d) as 6 = D/d. For 
<mall dead zone loss (Section 7.3), the loss is given by n „ y = (d/f)/(x /D). 
Rewriting the expression with the spot diameter given by D$ = 2.44 xf/D gives 
n = 2.44 d/D $ . 

The number of servo increments (N) is given by 

°F0V / ®CL = 
or 

D/D $ = 2 N 

Rewriting in terms of d and n» and using 6 = D/d, we obtain 

N = log 2 U«/2.44) 
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For 50 percent dead zone loss (n « 0.5), and using a * 83.2 for the Hughes 
HPIN 444Q quadrant detector, Ua/2.44) = 17, so that N = 5 is required. 
Although sometimes fewer steps will be needed. In general N * 5 will be 
required for the HPIN 444Q and n = 0.5. 

8.4.3 Detection and False Alarm Probabilities 


For a quadrant detector acquisition sensor, the transmitter dwell time 
must be sufficient for the quadrant detector to perform N angular reduction 
steps which transfer the spot into the closed loop tracking field of view 
(Section 8.4.2). Since during that time the spot is located in a single 
quadrant, the probability that a single servo step is accomplished correctly 

3 

u p i (i - p f a ) , where P. is the detection probability in the 
qjadrant with the signal present, and P^- is the false alarm probability in 
a quadrant with no signal. The probability that N servo steps are accom¬ 
plished correctly is p|j (1 - P^ a )^ N . This is the probability that 
the acquisition quadrant detector achieves closed loop tracking when it is 
illuminated by the transmitter. 

The overall acquisition sequence also requires that no false alarms occur 
from the quadrant detector during the time it is not illuminated by the trans¬ 
mitter. Since we have required N servo steps to occur within- a dwell time, we 
have set t d = Nt , where t^ is the dwell time and r is the effective 
integration time cf the quad detector. During a dwell in which the quadrant 
detector is not illuminated, there are four quad elements, each with N chances 

4N 

for false alarm, so the probability that no false alarms occur is (1 - P f ) . 
If the transmitter addresses a total of M acquisition cells, then we require 
no false alarms to occur in the M-l dwell times during which the quadrant 
detector is not illuminated. The probability that no false alarms occur 
during the M-l non-illuminated dwell times is then (1 - Pf ) 4N ^ M "^. 

The probability of a completely successful acquisition sequence ) 

is the probability of N successful servo steps during the dwell in which the 
quadrant detector was illuminated, and no false alarms during the other 
dwells, and is given by 
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P ACQ ' P d“ 11 ' P f« )3N (1 - p fa )4N( "' U 


p K (1 P \4NM-N 
= K d u ■ K fa 


If we then, require equal contributions from missed detection and all the 
chances for false alarm, we set the threshold such that 

P - (1 - P )4M-1 

K d " 11 K fa J 

Using the binomial expansion this becomes 

P d 1 - (4M-U P fa 

For a given detection probability with equal contributions from missed detec¬ 
tion and all chances for false alarm, the resulting false alarm probability is 
then given by 


Pfa = (1 - P d )/(4M-1). 

8.4.4 CW Integrating Approach 

There are three basic modulation/detection alternatives for quadrant 
detector acquisition. The first is the cw integrating approach, in which the 
acquisition transmitter is a cw laser, the detector Integrates the signal 

9 

photoelectrons, and then thresholds out most of the noise. The second 
alternative is to intensity modulate the laser transmitter with an audio 
sinusoid (tone). The optical quadrant detector is followed by a narrow 
bandpass filter around the audio tone, envelope detected, and thresholded. 
This approach is discussed In Section 8.4.5. The final alternative is to 
time-gate the receiver to the laser pulse width. This approach Is discussed 
In Section 8.4.6 in the context of the high peak power, low duty cycle, Laser 
Diode Labs laser array, as well as a single mode laser source. 
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In the cw Integrating approach, the quadrant acquisition detector 
receives cw Illumination from the transmitter during the dwell time. The 

dwell time (tj allows N angular steps to occur, so t d = N t , where r Is 

the Integration time per servo step. The servo sample rate (1/r) then 

determines the required detector bandwidth. Setting (1/r) equal to the 
bandwidth determined by the RC product of amplifier Input resistance and total 
detector-preamp capacitance gives the minimum thermal noise (Section 7.5). 
For the Hughes HPIN 444Q quadrant detector the capacitance Is 15 pf per 
element. For a preamp noise figure of 3 dB (T = 290), the thermal noise 
photoelectron counts are K-* 3.0 x 10^ (independent of r). For the HPIN 

444Q dark current of Ip = 100 nA per element, the dark photoelectron count 
rate is given by 

n Q = I D /e = 6.3 x 10* 1 e"/sec 

For the 15 cm LEO collecting aperture, worst case irradlance in the 8.7 mrad 
(0.5*) field of view due to the moon (L. = 2.7 x 10’^ Wm" 2 A - * - Table 
7.7-1), quantum efficiency n = 0.83, and other parameters as in 'Table 7.6-1, 
the background count rate per element is n g = 8.15 x 10 10 e'/stc- (Not® 
the background flux splits into four detectors.) For optical signal power at 
the detector (P r ), the signal count rate is n $ = (n/hvjP*.. Defining the 
rms SNR as the square root of the power SNR In Section 7.5, we have for a PIN 
quadrant detector 




\/ K S + K B + K D + k t 





n S + n B + n D + K T /t 


i 
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The count rates and counts are related by K, o n^T » except that the 
minimum K, Is Independent of as shown In Section 7.5. From Section 8.4.3, 
the probability of successful acquisition Is 




( 1 - 


P fa } 


4NM-N 


9 


where the threshold Is set such that P 


fa 


(1 - P ( j)/(4M-1). For overall 


acquisition probability of P 


ACQ 


0.99, and with N * 5 servo steps to 


achieve closed loop, we require Pj * 0.999. The false alarm probability 

then depends on M, the number of acquisition cells addressed by the trans¬ 
mitter. However, for a given P^, the required SNR Is not extremely sensi¬ 
tive to P^ . Thus we will assume a nominal P~ = 10" 8 , and adjust it 

later for the M for each specific case. 

o 

For P^ = 0.999 and P* * 10 , an rms SNR = 9.3 is required. 

[Kingston, Detection of Optical and Infrared Radiation ], The required signal 
count rate versus integration time is 


n s 


T (msec) 

9.3 x 

10 8 

0.1 

2.5 x 

10 8 

1.0 

7.9 x 

10 7 

10.0 

2.5 x 

10 7 

100.0 

7.8 x 

10 6 

1000 


g 

Note that the signal count rate is less than n = 10 e /sec, which is a 

s io 

negligible shot noise source compared to n D = 8.2 x 10 e*/sec and 

ll - ° 

n D = 6.3 x 10 e /sec. 

For T 1.0 msec, the dark current dominates the noise contributions, so 
that the signal count rate is n $ ~ 1 //t" for t >_ 1.0 msec. Rewriting the 
signal count rate using P^ = n (hv/n), the required power versus inte¬ 
gration time is given by 
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P r (pW) 

X (msec) 

258 

0.1 

69.3 

1.0 

21.9 

10.0 

6.9 

100 

2.2 

1000 


Again, P R ~ 1 for r l.o msec, due to dark current dominated noise. 

Note that the large dark current Is due to requiring a physically large 
detector to get a high diameter to dead zone ratio (a). If the dead zone were 
physically smaller, for a fixed a the array diameter could be smaller, which 

i 

decreases both the capacitance and dark current In proportion to the reduction 
In detector area. Since we were required to assume worst case background 
(except the sun) the best detector performance will be background noise 
limited. A quadrant APD would not help In this respect since the background 
noise is multiplied by the APD excess noise (F., = 3 for G = 50). However, a 
PIN quad with a small dead zone and dark current less than 10 nA per element 
could be background noise dominated. 

Ftr the Ideal background - dominated performance case, the required power 
versus Integration time is given below, with the power savings relative to the 
dark-limited HPIN 444Q. 

Power savings 

P R (pW) T (msecJ_ (dB) _ 

158 0.1 2.1 


27.2 

1.0 

4.1 

7.5 

10 

4.7 

2.3 

100 

4.7 

0.75 

1000 

4.6 


8-25 










The required power goes as P« „ 1/^r”for t ^ 10 msec. Thus a small dead 
zone PIN quadrant detector which could achieve a given a with smaller array 
diameter (and correspondingly less dark current) could save 4.7 dB In laser 
power. 

The parametric evaluation of the HPIN 444Q quadrant detector performance 

_q 

Is summarized In Table 8.4-1, where the Initial assumption of P^ * 10 
has been corrected to the true value for the number of acquisition cells (M) 
for that case. Note that t. > Nt , and N * 5 angular reduction steps are 
assumed. 

Several features of Table 8.4-1 are noteworthy. First, even though Inte¬ 
grating longer reduces the required power at the detector, since the dwell 
time Is also Increased, the number of acquisition cells must be reduced for a 
given acquisition time. Reducing the number of acquisition cells requires 
Increasing the transmitted beamwidth, which In turn Increases the required 
transmitter power. Thus integrating longer actually Increases the required 
transmitter power for a given acquisition time. 

Second, the cases for r < 10 msec would be difficult to Implement, since 
the angular step reduction would require very large accelerations and torque, 
as well as greater than 100 Hz bandwldths. Third, to use an acquisition 
beamwidth Oj < 30 prad Is probably Impractical. This Is because the 
spacecraft rms two axis attitude jitter Is about 17 prad, and'some fraction of 
that jitter Is coupled by glmbal friction Into line of sight pointing error, 
reducing control of beam overlap for a narrow beam. Then If cases of < 
30 prad and r < 10 msec are Ignored, the most viable case from the standpoint 
of laser power Is T^q * 2500 seconds, * 58 prad (Airy lobel and P^ = 
22 mW. Since present single transverse mode sources have 20 - 30 mW power, 
there would be essentially no margin. If a new PIN quadrant detector were 
developed, such that It could be background limited instead of dark current 
limited, a power savings of 4.7 dB would result for r * 10 msec. This would 
In effect provide the power margin of 4.7 dB with an acquisition time of 
Tacq * 2500 seconds (42 minutes) and P^ = 22 mW transmitter. Note that by 
restricting laser power to current commercially available levels, an un¬ 
realistic acquisition time of T^— * 42 minutes Is required. In fact, to 
obtain any power margin a developmental quadrant detector had to be postulated. 
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TABLE 8.4-1. HPIN 444Q QUADRANT 


DETECTOR ACQUISITION PERFORMANCE 


a ) ^ACQ * 2500 sec 


T (msec) 

0.1 
1.0 
. 10 
100 
1000 


b) T fi£q * iooo sec 
t( msec) 


0.1 

1.0 

10 

100 

1000 


c ) t ACQ = 500 sec 
T (msec) 


0.1 

1.0 

10 

100 

1000 


d) t ACQ * 250 sec 
T (msec) 


0.1 

1.0 

10 

100 

1000 


M 


5 x 

10 6 

5 x 

105 

5 x 

10 4 

5000 


500 



M 

2 x 

10® 

2 x 

10 5 

2 x 

10 4 

2000 


200 


M 


106 


105 


10 4 


103 


100 


M 



5 x 10$ 
5 x 10 4 
5 x 10 3 
500 
50 


«T(urad) 



9.2 

29 

92 

291 

919 


13 

41 

130 

411 

1300 


»l(nrad) 


18.4 

58 

184 

581. 

1840 


Pj(mW) 


3.1 

7.3 

22 

70 

210 


Pr(mW) 

7.2 

18 

56 

169 

513 


PftmW) 

14 ' 
37 
111 
334 
1015 


Py(mW) 


28 

72 

220 

659 

2000 











If the burden were shifted to the transmitter laser, assuming that 100 mW 
cw were available, an acquisition time of 1000 seconds (17 minutes) could be 
obtained, with 2.5 dB margin (Py > 56 mW required). If a 100 mW transmitter 
and a developmental quadrant detector (background limited) were assumed, 
Ta(Q " 5°° seconds (8.3 minutes) could be obtained with about 4.3 dB power 
margin (4.7 dB less than 111 mW required). However, this Is probably the 
upper bound on cw Integrating quadrant detector performance for acquisition, 
and the minimum acquisition time Is 8.3 minutes (100 mW transmitter, develop¬ 
mental background - limited quadrant detector). 

The results In Table 8.4-1 can be described (to first order) in compact 
form. The acquisition time Is given by T^.q « M t^, where M Is the number 
of acquisition cells and t^ Is the dwell time. Since t^ * Nr , where N Is 
the required number of angular reduction steps and r Is the Integration time, 
we can write T^q « HNt . The number of acquisition, cel Is (M) Is given by M 
■ (®uN(/ 9 T^’ wfiere ®unc aC( l u ^ s ^ t ^ on angular uncertainty, and 


Is the acquisition transmitter beamwldtn 
be written 


Thus the acquisition time can 



At this point It can be seen that the acquisition time Is proportional to the 

square of the uncertainty angle, l.e. the uncertainty solid angle. 

Using the link equation (Section 7.5), It can be seen that the required 

transmitter power Is proportional to the required detector power, and propor- 

2 

tlonal to the transmitted beamwldth squared, Py „ Oy. Since we showed 
P R . 1/yT", for r >_ 10 msec whether background or dark limited, we can 


wrl 


te Py ~ 9j/fr~, or Oy ~ Py , Replacing #y In the expression 


for T^.q we obtain 


ACQ 


-No 


UNC 


f> f T 
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Thi’s the acquisition time varies Inversely with laser power (since broader 
beamwldths can be used). Is directly proportional to the uncertainty solid 
angle, and Increases as the square root of the integration time. This 
relationship Is modified slightly when the dependence of on M Is taken 
Into account. 

The data from Table 8.4-1 1» plotted In Figures 8.4-1 through 8.4-7. 
Figures 8.4-1 through 8.4-5 plot required power (P-.) versus the Inverse of 
acquisition time (1/T^q) for fixed Integration times (t) from 0.1 through 
1000 msec. As can be easily seen, the deviation from Inverse proportionalIty 
due to the dependence of P fa on M Is slight. Figures 8.4-6 and 8.4-7 plot 
Pj versus t for fixed acquisition times of T^q «< 250 and 2500 seconds, 
respectively. The linearity of the log-log plots In Figures 8.4-6 and 8.4-7 
Indicates that the deviation from P- *yr In equation (1) Is small. 

The analysis of the cw integrating approach thus far has considered only 
the shot noise generated by the thermal, dark current, and optical background 
noise sources, and not the DC currents themselves. This Implicitly assumes 
that the DC currents do not vary In time, and can therefore be thresholded out 
(see Figure 7.7-2). While a threshold can Indeed eliminate the DC currents 
due to thermal noise and dark current (since this can je calibrated Into the 
receiver design), the DC current due to optical background changes as 
background sources move In and out of the receiver field of view. Star 
sources will be present In one quadrant, but not the others. Therefore the DC 
current due to optical background must be considered tor the cw Integrating 
approach. 

On the basis of shot-noise limited performance, the cw Integrating 
approach will be seen to be about 4.5 dB superior to the sinusoidal intensity 
modulation approach, and could be considered an upper bound to performance. 
However, the sinusoidal Intensity modulation approach Is not affected by 
optical background DC, and for practical systems will be superior to cw 
Integration for this reason. 
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Figure 8.4-1. Required Power (Pj) versus the Inverse of Acquisition 

Time (1/Tacq) f° r Integration Time r= 0.1 msec. 
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Figure 8.4 


1/T 


ACQ 
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2. Required Power (Pj) versus the Inverse of Acquisition 
Time (1 /T/\cq) for Integration Time r= 1.0 msec. 
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Figure 8.4-6. Required Power (Pr) versus Integration Time (t) for 

Acquisition Time T^q = 2500 Seconds. 
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Figure 8.4-7 


Required Power (Py) versus Integration Time (t) for 
Acquisition Time T/^q = 250 Seconds. 













It Is useful to consider an example of the impact of optical background 
DC on the cw Integrating approach. As In the tracking case, a uniform DC 
background cancels out In the sum and difference circuit, so the DC non¬ 
uniformity between quadrants Is the real Issue. However, the nonuniformity 
can be as large as the worst case background, since even an extended source 
could Illuminate only one of the four quadrants. 

When background DC Is considered, the SNR must be redefined as mean 
signal minus DC background divided by rms noise, l.e. 








Again defining count rates (n^) by 
time, the SNR becomes 


n, T, where Z Is the integration 


(n<- - n 0 ) 

SNR_= — —- — 

n S + n B + n D + K T /T ' 

Taking the same count rates for the HPIN 444Q as before (n n * 6.3 x 10** e‘/sec, 

in 7 _ 

n R * 8.15 x 10 e /sec, and K T * 3.0 x 10 e ), ,requiring SNR = 

D 10 
9.3 and T= 0.1 msec, the required signal is n^ * 8.25 x 10 , a 19.5 dB 

Increase compared to the case for shot noise only. (For longer T , the dB 

Increase Is even larger). If the optical filter bandwidth were reduced from 

50 to 10 angstroms (note that e R = e UNC * 0.5* and cannot be reduced), 

reducing the background count rate to n D = 1.63 x 10*^ e"/sec, n c = 

in “ ** 

1.72 x 10 e /sec are required (SNR = 9.3, 0.1 msec), a 12.7 dB 

increase over shot noise alone. 

For the LEO to GEO uplink In this study, with a large 0.5* diameter 
uncertainty angle, and acquisition required about twelve times dally, It Is 
probably not practical to wait for Intense background sources to exit the 
field of view. 
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Therefore a modulation scheme which eliminates the optical background DC 
Is necessary for quadrant detector acquisition. Although about 4.5 dB worse 
than cw Integration In shot noise performance, the sinusoidal Intensity 
modulation approach will be seen to eliminate the optical background DC 
problem. The performance In Table 8.4-1 (Figures 8.4-1 through 8.4-7 ) can be 
considered as the upper bound on shot-noise limited PIN quadrant detector 
acquisition, for cases where severe background sources can be avoided. 

8.4.5 Sinusoidal Intensity Modulation Approach 

This method of quadrant detector acquisition requires biasing the 
transmitter laser at Its half power level, and intensity modulating the laser 
with an audio sinusoid. For detection, the optical detector is followed by a 
narrow bandpass filter around the audio tone, and an envelope detector then 
determine the amplitude of the sinusoid (Figure 8.4-8). A threshold then 
removes most of the noise. We need to determine the performance of the 
modulation/detection scheme. 

The transmitter laser Is biased at half power (P~) and sinusoidally 
Intensity modulated, so that the power waveform is described by P(t) = 

P^d+cos u $ t). Note that the laser must have a cw power capability 
P CH * 2Pj. Although It might be argued that since sinusoid modulation 
Is 50 percent duty cycle, the diode average power could be set equal to the 

cw power, this Is valid only for frequencies above about 1.0 MHz. Audio 

frequencies appear to be cw to the diode. 

The photon Induced current at the optical detector Is 1 *t) = (neG 

P r /hv) (1 ♦ cos u s t) where P r Is the received power (related to P$ 

by the link equation - Section 7.6) and the other parameters are defined In 
Section 7.5. The average value Is 1 = neG P /hv, and this Is also the 

amplitude of the sinusoid term. 

The two-sided mean-squared current noise level (per element) Is given 
by (from Section 7.5) 


* 
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QUADRANT DETECTOR DEMODULATOR FOR ACQUISITION WITH 
SINUSOIDAL INTENSITY MODULATION 


ERROR 


DIFFERENCE 

CIRCUIT 


T-THRESHOLD 


Quadrant Detector Demodulator for Acquisition 
with Sinusoidal Intensity Modulation. 

























• ’tot ■ « G 2 ( s * e 6 2 1 B * . G 2 1 D 

♦ e1 K * 2 k T eq / R l (A 2 /Hz) 

The noise spectral density following the bandpass filter of transfer function 
Hgppff) Is 

V f) * ^TOT H BPF (f) 


Using a quadrature expansion of the noise we can write n(t) = n.(t) cos 

c 


(u $ t ♦ ) - n $ (t) sin (u $ t ♦ ^ ). 


The spectra of n„ and n, 


are 


obtained by shifting the one-sided bandpass spectrum to baseband, 


S „c ,fl * V fl 


2 ’tot I l W ,f 1 


where Hg p p(f) Is the bandpass filter function shifted to baseband. For a 
bandpass bandwidth from f - B/2 to f ♦ B/2, the noise power after the 

V* C 

bandpass filter Is given by 


/ 

B/2 


B/2 


2 1 TOT ^BPF (f) 


2 ^0T B 
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The envelope detector Is a nonlinear element, and alters the Gaussian 
probability distribution of the noise. When no signal Is present, the 
envelope Is given by 

n—p 

a(t) * y n /» 4 n e » an< ^ ^e enve 1°P e obeys Rayleigh statistics. 

The Rayleigh density for the envelope a Is given by 
P a (x) « (x/c 2 ) exp (-x 2 / 2 a 2 ) 

When the signal Is present, the statistics become Rician. The envelope of a 
sinusoid of amplitude 1 $ added to Gaussian noise Is 

2 2 l 1/2 

♦ n c (t) ] Z ♦ n 2 (t) J 
and the Rician density for the envelope Is 

* 

P a (a) « (x/o 2 ) exp [ - (x 2 + 1 $ 2 ) / 2 o n 2 ] I Q (Ax / o n 2 ) 

which Is valid for x > 0, and 1 Is the modified Bessel function. 

— 0 

False alarm and detection probabilities for Rayleigh and Rician statis¬ 
tics are tabulated In I ntroduction to Radar Systems by M. Skolnlk (McGraw-Hill 
1962) as well as many radar texts. 

1 ~9 

Since the power in the signal sinusoid In P $ - ( /2) 1‘, Skolnlk 

defines the signal to noise ratio S/N as 



S/N . (1/2) 1 2 / o 2 

s n 
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Inserting the expressions for l Jt o n 2 and Ijgj , and rearranging, gives 


1/2 (n P R / hv) 


S/N 


♦ < B * V 2F N B/e * 


2 V 4 k T eq B 


e G‘ 


e 2 G 2 R, 


Using the acquisition probability expression from Section 8.4.3, for acquisi¬ 
tion probability of * 0.99, equal contributions from missed detections 
and all chances for false alarms, and N « 5 servo steps, we require Pj « 
0.999 and P^ « (1 - P rf )/(4M - 1), where M Is the number of acquisition 


cells. 

Since the required S/N Is not extremely sensitive to Pr,, we will use a 

-12 ra 
nominal P. b 10 , and later revise the S/N for the exact M and P^ a 

for that case. For P d m 0.999 and P^ b 10’* 2 , S/N b 17.5 dB Is 

required (Skolnlk). 

Again we assume parameters for the Hughes HPIN 444Q quadrant detector, 
and worst case background Illumination due to the moon with L T * 2.7 x 10" 7 

-2 t -1 1 

Wm A . We will also assume that the bandpass filter has a bandwidth of 
100 Hz. Note, however, that the detector bandwidth has to respond to the 
sinusoid modulation frequency f , and minimum thermal noise Is computed on 
the basis of detector bandwidth f (section 7.5.2). Therefore f should 
be as low In frequency as possible, yet must be high enough to get outside the 
detector */f noise. We will assume f ■ !0 kHz, which Is a low enough 
frequency that thermal noise Is not the dominant noise source for the HPIN 
444Q. (Tht relative contributions of background, dark current, and thermal 
noise are In the ratios 1.6, 12.5 and 6.3, respectively). 

Note that If a modulating frequency such as f b 1.0 MHz were chosen, 
so that the diode average power could be set equal to Its cw power (Instead of 
P cw b 2P $ ), the thermal noise would be Increased a factor of 100 (compared 
to f $ b 10 kHz.) For fixed capacitance, R^ must be reduced to Increase 
the bandwidth, and the thermal noise spectral level Is 2 K T e /R^. 
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Signal Induced noise Is negligible due to the high contributions of the 
other noise sources. For S/N > 17.5 dB (P d ■ 0.999, P^ a « 10"^), 

P R ■ 42 pW Is required. 

Table 8.4-2 summarizes the acquisition performance of the Hughes HPIN 
444Q with sinusoidal Intensity modulation. P^ Is the cw power requirement 
of the laser transmitter, given by P^ ■ 2 P , where P $ Is the amplitude 
of the sinusoid In the power waveform. Table 8.4-2 was calculated with a 
fixed dwell time t^ > N/B, where N Is the number of servo steps for angular 
reduction, and the servo response time Is 1/B, where B Is the bandwidth of the 
bandpass filter. For N « 5 and B - 100 Hz, t^ ■ 0.05 seconds. 

TABLE 8.4-2. ACQUISITION WITH SINSUSOIDAL INTENSITY MODULATION 

(HUGHES HPIN 444Q) 


T ACQ ,snI 

Oy(nrad) 

M 

P rw 

CW 

1 

2.1 x 10 4 

• 

20 

• 

4.2 x 10 5 

1 

« 

9.2 1 

1.0 x 10 4 

29 

2.0 x 10 5 

. 17.9 j 

5027 

41 

1.0 x 10 5 

35.3 

2000 

65 

4.0 x 10 4 

87.4 

956 

94 

1.9 x 10 4 

171 ! 

471 

; 134 

9412 

338 

188 

212 

3760 

835 

93 

• 

302 

1853 

1650 


Note that the sinusoidal Intensity modulation (SIM) approach (Table 
8.4-2) has poort*' performance than the shot noise limited cw Integrating (CWI) 


approach (Table 8.4-1, for t . 10 msec). T 


ACQ 


1000 seconds requires about 
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164 mW (Interpolating by T^q _ Vp cw )» wh ^ e f° r CWI T^q ■ 1000 
seconds Is achieved with 56 mW, a savings of 4.7 dB. Similarly, T^q = 500 
seconds requires about 318 mW for SIM, and 111 mW for CWI, a 4.6 dB savings. 
CWI would be a superior approach tr SIM If optical background DC were fixed 
and could merely be thresholded out. However, as described In the previous 
section, the CWI approach pays a 12.7 dB or greater penalty for optical 
background moving in and out of the receiver field of view, as would occur for 
practical satellite receivers. Therefore, by eliminating the background DC 
current, SIM Is seen to be the superior approach for quadrant detector 
acquisition. 

Figure 8.4-9 plots required power versus the inverse of acquisition 
time (1/T acq ) from Table 8-4 for the HPIN 444Q. Since the log-log plot 
deviates only slightly from a slope of unity, the inverse proportionality is 
only slightly modified by the dependence cf P f on M. 

A savings of 2.0 dB in laser power from Table 8.4-2 is obtained if a 
developmental quadrant detector is assumed, which has a smaller area (for a 
given a) so that dark current is negligible. (This is less than the 4.7 dB 
savings which results by removing dark current from the shot noise limited CWI 
case, since SIM has additional thermal noise from the requirement that the 
detector bandwidth accomodate the 10 kHz tone.) 

It is interesting to qualitatively compare the shot noise limited cw 

integrating approach (CWI) with the sinusoidal intensity modulation (SIM) 

approach. SIM requires biasing the laser at half power, so that the sinusoid 

amplitude is given by half the laser's cw power capability - effectively a 3 

dB loss. The receiver pays an additional 3 dB penalty, since the power in a 

2 

sinusoid of amplitude A is 1/2 A . However, the Rayliegh/Rician statistics 
give the same detection and false alarm probability performance at about 1.5 
dS less in rms SNR than for Gaussian statistics. The net effect is that shot 
noise limited CWI could achieve the same acquisition time with about 4.5 dB 
less laser power. However as previously mentioned, the optical background DC 
current CWI problem eliminates this advantage. 
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Figure 8.4-9. Required Power (Pew) versus the Inverse of the Acqulsiti 

Time (1/T/^eg) for the HPIN 444Q Quadrant Detector with 
Sinusoidal Intensity Modulation. 


















Since for the SIM approach If dark current Is removed thermal noise 
becomes Important (unlike shot noise limited CWI, which was background limited 
with dark current removed) it Is of Interest to see If a quad APO could 


improve the SIM approach. 

Assuming the developmental RCA quadrant APD with dead zone reduced from 
d = 75 um to d = 25 urn, « = 60. The number of servo angular steps required 


are N * log- (n «/2.44), which gives N = 4 for 50 percent dead zone loss. 


For P ACQ * 0,99 ’ 
Section 8.4.3, we 


using the detection and false alarm expressions from 

_g 

require P d = 0.999, and assuming a nominal P^ = 10 , 


S/N s 16.5 d8 is required (Skolnik). 


Again assuming moon background, and gain G = 50 (F^ = 3), performance 
is background limited as might be expected for an APD. Note It Is desirable 


to become background limited with a PIN detector If possible, because the 
background is multiplied by F N for an APD, which is a factor of /FT 
additional laser power required. In a 100 Hz bandwidth, for G = 50, 19.7 pW 
are required for S/N = 16.5 dB, but reducing the gain to G = 25 (F^ = 2.5) 
reduces the required power to P- = 18 pW. 

Table 8.4-3 summarizes the acquisition performance of SIM with a quadrant 
APD. The dwell time is tj = N/B = 0.04 seconds. Comparing Table 8.4-3 to 
Table 8.4-4 (HPIN 444Q), at = 160 seconds the APD quad requires 374 mW, 
while the PIN quad requires about 872 mW (using P cw ~ 'l/ T A r Q t0 * nter " 
polate), a difference of 3.7 dB. T^q = 500 seconds requires about 136 mW 
with the quad APD and 318 mW for the PIN quad, also a 3.7 dB difference. 


Figure 8.4-10 plots required power (P-^) versus the inverse of 
acquisition time (1 / t acq^‘ Again the log-log plot has nearly unity slope, 


for a near linear relationship between P_ w 


and (1/T AC q)- 
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Figure 8.4-10. Required Power (Pcy) versus the Inverse of Acquisition 

Time (1 /Tacq) for an RCA Quadrant APD with Sinusoidal 
Intensity Modulation. 





























TABLE 8.4-3. ACQUISITION WITH SINSUSOIDAL INTENSITY MODULATION 

(RCA QUADRANT APD) 


TACQ(sec) 

ejturad) 

M 

Pew 

7660 

30 

1.9 x 10 5 

10 

3830 

42 

9.6 x 10 4 

: 20 

1878 

60 

4.7 x 10 4 

39 

734 

96 

1.8 x 10 4 

96 

360 

J 137 

9000 

189 

180 

194 

4500 

374 

II 70 

311 

1750 

914 

34 

445 

850 

1785 


Let us also derive a compact (first order) expression for acquisition 

time with sinusoidal intensity modulation. As for the shot noise limited cw 

integrating case, T^p = M t d , where M is the number of acquisition cells 

and t d Is the dwell time. The dwell time Is given by t d = N/B, where B is 

the servo bandwidth, but we will consider t d fixed at Its shortest value 

(0.05 seconds for the HPIN 444Q) to avoid the dependence of required power on 

noise bandwidth. Again M = (®Mur/®T * where Is the. angular 

acquisition uncertainty and e T is the transmitted beamwidth. Since the 

required transmitter power is proportional to ® T (link equation - section 

2 

7.6) we can write P^ = y ® T , with the result 


T ACQ = y 9 UNC V P T’ 


Again the acquisition time is proportional to the uncertainty solid angle, and 
Inversely proportional to laser power. This relationship is slightly modified 
when the dependence of P f . on M Is taken into account. 
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8.4.6 Time-Gated Acquisition 

For this detection scheme the acquisition laser source Is operated 
pulsed, and the acquisition receiver establishes time slots equal to the pulse 
width. The receiver Integrates over each time slot, and compares to a 
threshold to determine If a pulse was present In that time slot. Since the 
transmitter and receiver do not share time synchronization during acquisition, 
the desired detection performance must be achieved in the worst case for which 
the pulse energy Is shared equally between adjacent time slots {1.e. can only 
utilize half the pulse energy). This scheme thus utilizes a higher bandwidth, 
more complex acquisition receiver than the cw integration (CWI) or sinusoidal 
Intensity modulation (S/M) approaches. S1i.ce the lowest frequency component 
of the pulsed waveform Is approximately (l/T ), where T is the pulse 
repetition period, the detector output can be AC-coupled to avoid the problem 
of the optical background DC component. 

The detection and false alarm probabilities for the time-gated acquisi¬ 
tion (TGA) scheme are somewhat different than for CWI or SIM (Section 8.4.3). 
As before, there are N servo steps which must be correctly performed to center 
up the Illuminated spot on the quadrant detector. The dwell time (tJ is 
then given by t d = N x , where t Is the servo integration time. Now assume 
there are R laser pulses per servo integration time, where R s T /T rf p. The 
probability that at least one of the R pulses is correctly detected is given 
by [l-(l-P d ) R ], where P d is the detection probability within a time 
slot. (This is one minus the probability that none of the R pulses were 
detected.) 

The total number of time slots per servo integration time is given 
by t / Tp, where tp is the laser pulse width (same as time slot width.) 
Since there are 4 detectors in the array, the number of chances for false 
alarm is 4 t / r . Since N successful servo steps are required, the 
probability of a successful event during the illuminated dwell time is 

ci - (i-p d ) R ] N d-p fa ) 4N t/t p 
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The probability that there are no false alarms during the (M-l) non- 
llluminated dwell times (M Is the number of acquisition cells) Is 


Thus the probability of successful acquisition Is 


P ACQ * tl-ll-p/j" •l-’V 


4NM V T p 


Since the acquisition transmitter and receiver do not share time synchroniza¬ 
tion, the desired acquisition probability must be obtained with only half the 
pulse energy within the time slot. 

For equal contributions from missed detection and false alarm, we set 
l-(l-P d ) R = (l-P fa ) 4M T/ T p 
Using the binomial expansion, this becomes 


l-(l-P d ) R = 1-C4M T / T p )P fa 
U-P d )R 

P fa = T4 M't/t"J 

P 

First we examine the case using the Laser Diode Labs (LDL) LS-410 
Incoherent laser array (Section 8.3). The maximum repetition rate Is 5 kHz 
(T s 200 psec) and the maximum pulse width Is 130 nsec. 

For a 9S percent acquisition probability, and equal contributions from 
missed detection and false alarm, we require 

[l-(l-P d ) R ] N = \j 0 . 99 " 
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For N s 5 servo steps (Section 8.4.2), servo Integration time t * 0.01 seconds 

(B . 100 Hz), and T ' ■ 200 usee (R « t /T rep - 50 pulses), the required 

detection probability Is Pj = 0.13. This Illustrates that the detection 

probability requirement Is low since we only need detect one out of 50 

pulses. Assuming a maximum M > 2.0 x 10^ acquisition cells (M = (13.7 x 

10" 3 /30 x 10“ 6 ) 2 ), t * 0.01 sec, and t = 130 nsec, the required 

-14 P 

false alarm probability Is Py « 10 . For P^ * 0.2 and P^ a * 

10” , an rms SNR « 7 Is required. 

For an RCA quadrant APD, with values as In Section 7.8.1, G * 50 
(Fjg « 3), moon background (Table 7.7-1), o R a 0.5*, x p = 130 nsec, 2.0 
pf per detector-preamp element, and parameters as In Table 7.7-1, the noise 
counts are K~ = 0.2, K^ = 0.8, Kg s 38,600 and Ky 8 1576. For rms 
SNR s 7 we require 



which then requires K$ = 2500. Since up to half the pulse energy Is wasted 

due to lack of time synchronization, the pulse : tt nerate K^ = 5000. The 

required power Is P R s K s (hv/n *C ) a 11.5 nW. Since in' Section 8.3 we 

found that the LDL array can be collimated to deliver only about 5.0 pW peak 

power to the detector, the array has about 33.6 dB less power than required. 

Now let us consider a single mode laser with 20 mW average power, such as 

the RCA C86030E. We will assume the same pulse width ( T p = 130 nsec), but 

for duty cycles of 12 percent or less the RCA can deliver up to 160 mW peak 

power. The false alarm requirement Is the same, but since now T « 1.0 

4 rep 

nsec, there are 10 pulses per servo Integration time. The required 

detection probability Is only P d 8 6.9 x 10”*, with P f 8 10” 14 as 

before. Using P d 8 0.01 and P^ = 10” , an rms SNR s 5.5 Is required. 

With the same noise counts as just evaluated, and SNR = 5.5, K^ 8 2000 are 

required. Thus the total photoelectrons per pules are K^ s 4000, due to 

lack of time synchronization. This corresponds to a required power of P R s 

9.2 nW. 
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For P R » 9.2 nW, ® T = 30 prad (airy lobe), and other parameters as In 
Table 7.7-1, a transmitted power of Pj ■ 2.6 W Is required. This Is about 
12 dB more than the 160 mW peak power capability of the RCA diode. However, 
it can be seen that the LDL array, In spite of Its 1.0 kW peak power, Is about 
20 dB Inferior to a single mode laser. This Is because the large source size 
of the LDL array makes the output difficult to collect and collimate. 


8.4.7 Quadrant Detector Acquisition Sequence 


To avoid the earth as a background source, the geosynchronous (GEO) 

satellite Is designated the acquisition transmitter, and the low earth 

orbiting (LEO) satellite Is the acquisition receiver. Initially, the GEO 
scans Its transmitted beam over the 0.5* uncertainty cone, while the LEO 
acquisition quadrant receiver stares at the entire 0.5* uncertainty cone of 
the GEO. Although the GEO scan Is continuous, the time to complete the scan 
can be estimated by T^q = M t^, where M is the number of acquisition 
cells and t^ Is the dwell time. M is given by M = (®wN(/®r ’ where 
®UNC the ac< l u ^ s1t1on uncertainty angle and e T is the ’transmitted 

beamwidth. This scan time can be loosely called the acquisition time, since 
It Is by far the dominant contribution. 

When the LEO quadrant detector senses that one quadrant has been 

illuminated, the servo angular step reduction sequence (Section 8.4.2) is 

Initiated. Once closed loop tracking Is obtained the field of view is then 

reduced to (Section 8.4.1). The (worst case) time for the acquisition 

sequence to this point Is M tj. 

It should be noted that the SNR which determines is different than 

the SNR on which the false alarm and detection probabilities are calculated. 
This Is because prior to angular step reduction, the light spot is in a single 
quadrant (typically), but once closed loop tracking begins, the spot becomes 
divided between four quadrants and additionally sees dead zone loss. For 50 
percent dead zone loss, this represents a factor of eight reduction in power 

per quadrant for closed loop tracking. This corresponds to a factor of 8 

(9 dB) reduction in rms SNR, and a factor of 64 (18 dB) reduction in power 
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SNR. The rms SNR per quadrant for acquisition (cw Integration approach) was 
about 9.3 for P, ■ 0.999 and P fa - 10 . After a factor of eight loss, 
the rms SNR Is only 1.2, which Is used to compute # NEA . For SNR » 1.2 the 
field of view reduction ratio. Is R » 3.1 D/D^ « 1.27 n a. For 50 percent 
dead zone loss, R ■ 0.64 a. Thus the quadrant detector diameter to dead zone 
ratio Is the critical parameter for FOV reduction. For the Hughes HP IN 444Q 
a « 83.2, so that R « 53.2. Reducing the 0.5* (8.7 mrad) by R * 53.2 gives a 
GEO uncertainty FOV at the LEO of « 164 prad. 

From the point of having achieved a closed loop ® NEA estimate at the 
LEO, the acquisition sequence proceeds as follows. By the time the servo 
angular step reduction Is completed at the LEO, the GEO transmitter has 
essentially moved to address the next acquisition cell. The GEO acquisition 
quadrant detector must have a field of view at least three or four acquisition 
cells wide, so that the transmitter can move on to the next acquisition cell, 
without having to Include more time for the return beam within each dwell time. 

The LEO now knows that the GEO Is located within « NEA = 164 prad, and 
directs its 10 trad (1/e^) communication beam to scan e N r A . With a dwell 
time of t^ = 0.05 seconds, the GEO can perform angular step reduction with 
Its acquisition quadrant detector. This scan time will take approximately 
T$can s 300 (0.05) = 15 seconds. 

As soon as the GEO quadrant detector is satisfied that It has been 

4 

Illuminated (e.g. two successive quadrant hits) the GEO acquisition trans¬ 
mitter Is directed to stop scanning, and is gradually directed toward © NEA 
during the dwell time, as determined by the GEO quadrant detector. (g nea Is 
about four acquisition cells reduced by R » 53.2). At this time the 35 prad 
beacon transmitter is turned on, which shares the aperture with the acquisi¬ 
tion transmitter by orthogonal polarizations. If © NEA (GEO) Is large, a 
small scan may take place, otherwise e^ EA is Illuminated continuously. 
Using the beacon laser allows the LEO acquisition quad to get a boosted rms 
SNR by the ratio of the acquisition transmitter beamwidth to beacon beamwidth 
ratio squared. The LEO acquisition quad obtains additional FOV reduction due 
to higher SNR, so that the spot falls through a bifurcated mirror onto the LEO 
tracking quad. The tracking quad performs an angular step reduction, which 
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then results In closed loop tracking. The 10 prad LEO communication beam then 
continuously Illuminates the LEO. Both the acquisition quadrant FOV reduction 
and tracking quadrant angular step reductions on the LEO require a total of 
about 0.1 seconds. The GEO then performs FOV and angular step reduction on 
Its quadrant detectors, requiring another 0.1 seconds, resulting In mutual 
tracking and communication. 

The amount of time which should be added to the times In Tables 8.4*1, 
8.4-2, and 8.4-3 Is 15 seconds plus 0.1 seconds plus 0.1 seconds plus 0.6 
seconds, where the last term represents four d/c time of flight delays. Thus 
about 16 seconds should be added to the acquisition times In Tables 8.4-1, 
8.4-2, and 8.4-3. 













8.5 ARRAY ACQUISITION 


8.5.1 r'jld of View Reduction 


For array acquisition, the linear field of view (FOV) reduction Is given 
by /nT where N Is the number of pixels In the array. With a 400 x 500 
array, the FOV can be reduced by a factor of 400, from an uncertainty angle of 
0.5* (8.7 mrad), to 22 urad. This Is nearly eight times greater than the FOV 
reduction of 53 which Is practical for a quad cell with low SNR. 

If the array Is being used with a video tracker, additional FOV reduction 
can take place. If the optical blur Is set to subtend two by two pixels (for 
good tracking performance) the angular tracking jltte 1 * (for high rms SNR) Is 
given by 



0.685 a 
SNR - 3 


where a Is the angular subtense of the optical blur. We define © NEA * 2 
oj similar to the quad cell case. Since the blur angle Is twice the single 
pixel angle, the uncertainty Is = /TT a/2, where i/n" is the one 

dimensional number of pixels In the array. Defining the FOV reduction ratio 
as R * ®UNC^®NEA’ ^ reduction Is R « 0.365 /n* (SNR -3). 

8.5.2 Dete ction and False Alarm Probabilities 

All CCD-type Imaging arrays effectively Integrate the photocurrent In 
each pixel - they vary only In how the collected charge Is stored and read 
out. The choice of detection scheme Is to threshold the collected charge In 
each pixel, removing much of the noise. A simple circuit on the video output 
of the array electronics can accomplish the thresholding (a greatly simplified 
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version of current Hughes video trackers, for example.) This Is similar to 

the Integration and threshold scheme used for tracking (Sections 7.4, 7.8) and 

to the cw Integration scheme for quad cell acquisition (Section 8.4.4). The 

detection performance can be characterized by the probabilities of detection 

and false alarm, as for the quad cell acquisition schemes. 

Let the array consist of N pixels ( /¥ x /n~ ), and the acquisition 

uncertainty consist of M acquisition cells (\Tm x^M ), where M Is given by 

2 

M « (®uN(/*r as ^ or < l ua< * rant detector cases. The array stares 

Into the entire uncertainty field of view, where In the worst case, It must 

wait (M-l) seconds before It Is Illuminated ( t Is the Integration time). 

Note that the Integration time Is equal to the dwell time, since no step 

angular reduction Is required as for the quad cell. Since there are M-l 

chances for all N pixels to Incur a false alarm, and one chance for N-l false 

alarms, there are a total of MN-1 chances for pixel false alarm, but a single 

• * 

chance for detection. 

MN 1 

The probability of a successful acquisition Is P^q * P^ (1 - P^) 
where Pq Is the detection probability and P fa Is the false alarm proba¬ 
bility. The expression for P^q results since we require a single correct 
detection, and no false alarms. For equal contributions to error by missed 
detection and all chances for false alarm, the threshold Is set such that 

p M p \MN-1 
P D * 11 • P fa 

Using the binomial expansion, this Is equivalent to setting 


fa 


(1 - P D ) / (MN-1) 


An upper bound can be obtained on MN as follows. For a 400 x 500 array 

5 

N b 2.0 x 10 . Since e-j. • 30 urad Is about the minimum beamwldth 
desirable to use for acquisition, and the uncertainty angle Is 0.5* (8.7 


mrad), an upper bound on M is obtained by M 
An upper bound on MN Is then about 4 x 10 . 


(e UNC /# T l 


1.9 x 10' 
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-13 

For* P^gQ « 0.99, we require P Q ■ 0.995, and then P p# ■ 10 

based on the upper bound of MN. Since the required SNR Is not extremely 

-13 

sensitive to P pa , we will use the nominal P fa ■ 10 and later correct 
the SNR for the exact case. For P^ > 0.999 (better than required) and P p 
■ 10* , an rms SNR > 10.6 Is needed. 


8.5.3 Array Acquisition Performance 


As mentioned In Section 8.5.2, nominal detection and false alarm 

-13 

probability of Pp - 0.999 and P. ■ 10 will be assumed, which results 

In better than « 0.99 overall acquisition probability. The false alarm 
probability will later be adjusted for the number of acquisition cells In the 
specific case, according to the P FA expression In Section 8.5.2. 

We will assume an array size typical of currently available commercial 
CCO arrays (Section 8.2.1), and use 400 x 500 pixels, so that N ■ 2.0 x 10^. 
Since the overall acquisition uncertainty Is e UNC « 8.7 mrad (0.5*), the 

single pixel FOV Is (8.7 mrad / 400) » 22 urad. Note that the planet 
background sources (Tab.e 7-2) subtend at least 4x4 pixels, and a simple 

Image processing circuit could discriminate against planets since they are 
extended sources. Note that stars are unresolved sources, and could not be 
discriminated. However, we will assume that there Is no extended source 


discrimination, so that the array must detect with the worst case background 
Irradlance (except the sun) in the 22 urad FOY. 


Using e D = 22 urad, the extended sources have the following Irradlances 

_2* _1 K .7 .12 

In Wm A : sun, L, ■ 6.8 x 10 ; moon, L. & 1.7 x 10 ; Venus, 

-12 1 -13 1 

Lj ■ 6.8 x 10 ; Jupiter, L. ■ 1.7 x 10 . These can be compared to 

Lj « 4.8 x 10"^ for Sirius and Lj « 2.1 x 10"^ for Mira. Note that 

the extended sources and the stars have very comparable Irradlances for this 

22 urad FOV. If the acquisition uncertainty were larger, or there were fewer 

pixel elements In the array, the pixel FOV would be larger, and the extended 

sources would clearly dominate the Irradlance. Even at the pixel FOV of only 

22 urad, Venus Is the strongest Irradlance. For the parameters (Table 7.6-1) 

0 R » 15 cm, N r . 0.54, b - 5oi, N p . 0.7, n ■ 0.25, hv « 2.3 x 10* 19 
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-12 -2 -1 

joules, lj = 6.8 x 10 Km A and t = 10 msec, the background photocounts 

are K 0 * 2.5 x 10 . Note that the assumed quantum efficiency of 25 percent 

Is somewhat conservative, since some CCO arrays approach 50 to 60 percent for 
the 0.8 - 0.9 urn region. 

The maximum number of thermal electrons are Ky = 300, which Is both a 
commercial manufacturer's specification and value obtained from the Hughes CCD 
group. The thermal electrons are Independent of Integration time (t). 

The typical CCD array has a dark current (at 15*C) given by I n = 20 

2 u 
nA/cm . (This could be lowered by cooling.) For typical pixel array 12 pm 

x 18 urn, and = 10 msec, the photocounts due to dark current shot noise are 

3 

K D — 2.7 x 10 . However, an array has an additional dark current problem 

that single photodetectors do not have - dark current nonuniformity. That is, 

the variation in dark current from pixel to pixel will cause a "fixed pattern 

noise" for long integration, given uniform illumination. The typical CCD 

o 

array dark nonuniformity is 1.0 nA/cm. (15 C), which for 12 urn x 18 pm pixel 
area and t * 10 msec integration time, results in = 135 rms photocount 

variation. 

The rms photocount variation due to shot noise is 
(K s ♦ Kg ♦ K d ♦ M 1/2 = (3 x 10 4 ) 1/2 = 174. 

Root-square-summing the shot noise and dark current photocount variation gives 
o = C(174) 2 ) ♦ (135) 2 ] 1/2 = 220.. 

The CCD array has an additional noise problem with the DC current due to 
optical background. The DC components due to dark current and thermal noise 
can be removed with a fixed threshold. Optical background sources move in and 
out of the array field of view, however. The standard approach to removing 
background is to establish a threshold above the worst case background level. 
Since background would not always be present in every pixel, the signal is 
required to have enough signal to exceed the threshold without the noise mean 
added to it. From this viewpoint, the rms SNR is given by 


SNR 


rms 
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where o Is the combined rms variation due to shot noise and dark current 

nonuniformlty. 

4 

For x * 10 msec, o * 220 and K R = 2.5 x 10 . To achieve rms 

4 D 

SNR = 10.6, Kj = 2.7 x 10 are required, corresponding to a received power 

P R » K s (hv/n t) =2.4 pW. 

Using the link parameters of Table 7.7-1, with 6y = 130 prad (Airy 

lobe), and T s = 1.0 (no split), a transmitted power of Py = 12.9 mW is 

required. A laser source with Py = 51 mW would satisfy the acquisition 

requirement with 6.G dB power margin. 

2 4 

For 9y = 130 prad there are M = (13 mrad/130 prad) = 10 

acquisition cells. (Note a transmitter needs Oy = 13 mrad Airy lobe width 
to adequately cover the 8.7 mrad uncertainty - see Section 8.6.1). The dwell 
time could be set equal to the integration time if the GEO acquisition 

quadrant detector had a field of view several acquisition cells wide, so that 
it could see a response after moving on to address the next cell. Alterna¬ 
tively the servo time to address the illuminated pixel could be included in 

the dwell time. Since the former approach provides faster acqusition, we set 
the dwell time equal to the array integration time. * 

4 

For Oy = 130 prad, there are 10 acquisition cells, and t d = t = 10 
msec dwell time, so that the acquisition time (neglecting handover) is T^q 
= M tj = 100 seconds, with P A( .q = 0.99 and 6.0 dB power margin with a 50 
mW source. For comparison, acquisition with an AP0 quadrant detector was 
T ACQ = ^34 seconds with no power margin (SIM approach), and assuming a 96 mW 
cw source (see Table 8.4-3). The array has a factor of 7 shorter acquisition 
time with a simultaneous 9.0 dB power advantage. Using the tremendous 

advantage of the array over the quadrant detector is the only way to achieve 
acquisition times below two minutes, for currently available laser power. 
Table 8.5-1 shows some alternate parameters for array acquisition performance. 

If the dwell time is set equal to the integration time, and integration 
time considered fixed, the acquisition time with an array is given by 

T ACQ = *®UNC 1 9 T )2 


8-59 










TABLE 8.5-1. CCD ARRAY - ACQUISITION PERFORMANCE 


10 msec 


Pj (mW) 

07 (urad) 

T acq (sec) 

50 

260 

25 

100 

[ » _ 

! 260 

25 

50 

130 

100 


Margin (dB) 


100 


184 


50 


6 . 























Defining P- = y e^, where y absorbs the other link equation parameters, 
the acquisition time becomes 

t acq * y9 unc t/p t 

The acquisition time Is proportional to the uncertainty solid angle, and 
Inversely proportional to laser power., just as In the quadrant detector case. 
Thus the factor of 7 shorter time and 9.0 dB power advantage of the array over 
APD quadrant detector acquisition (SIM) can be considered a 17.5 power 
advantage. 

It Is also of interest to knew the effect of varying the Integration 

time. Since thermal noise was >ot a dominant noise source, the shot noise 

variance (counts squared) Is proportonal to t . The dark current non- 

2 

uniformity variance is propor:cnal to t . The rms SNR can be rewritten as 



2 ? 
where og N is the variance due to dark nonuniformity and is the 

variance due to shot noise. Since the shot noise variance is equal to the 

2 2 2 
number of counts, o N _ K^. The dark current variance is oq N = k^. 

Defining = n^T, where n Is the count rate, the SNR becomes 



(n$ - n B )x 

2 t 2 

(n N X + n DN ) 




♦ 


V T > 



Thus the dark nonuniformity Is the dominant noise source for long Integration 
times, and negligible background, as has been experimentally observed. 
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(Note: T"< 1.0 second to avoid charge saturation.) Note that the signal count 
rate must exceed the DC background count rate without dependence on Integra- 
tlon time. For 1 nA/cm dark current nonuniformity, and 12 urn x 18 urn pixel 
area, the dark nonuniformity count rate Is n^ = 13500 e‘/sec. The 
thermal noise (Kj = 300) Is negligible compared to the background and dark 
current shot noise. The noise count, rate can then be set as n N = n B + n D a 
2.8 x 10® e”/sec (n g = 2.5 x 10®, n Q = 2.7 x 10®. For SNR = 10.6, the required 
power versus Integration time Is (n = 0.25) 


(pw) 

T (msec) 

2.8 

1.0 

2.5 

10.0 

2.4 

100.0 

2.4 

1000.0 


The effect of Integration time is slight because the required signal count 
rate Is dominated by the fixed DC background rate. 

In summary, the array acquisition performance was seen to be superior to 
the APD quadrant detector with SIM by a 17.5 dB laser power' advantage, which 
could also be considered a factor of 56 shorter acquisition time. Since the 
CCD array Is required to allow intersatellite laser comnunication acquisition 
times less than two minutes, It is highly advisable to begin space qualifying 
currently available CCD arrays. 

8.5.4 Array Acquisition Sequence 


The acquisition sequence with an array is somewhat simpler than with a 
quadrant detector. Initially, the GEO acquisition transmitter scans a 130 - 
260 urad (Airy lobe) beam over the 0.5* uncertainty cone, which requires 
M = 2500 - 5000 acquisition cells to be addressed. (Again, although the scan 
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Is continuous, "addressing acquisition cells" Is a convenient description.) 

The 400 x 500 CCD array acquisition receiver on the LEO stares Into the 0.5* 

uncertainty cone. When a pixel Is Illuminated on the LEO array, the glmbal Is 

2 

directed to point the 10 prad (1/e ) communication beam to the 22 prad FOV 
location of the illumln- ated pixel. Since the APD quadrant tracking detector 
on the LEO has a 140 prad FOV, subsequent Illumination from the GEO will fall 
through the central hole of a bifurcated mirror onto the tracking detector. 

Although the GEO transmitter has by now moved on to address the next 
acquisition cell, the GEO acquisition quadrant detector has an P 0V which 
Includes several acquisition cells. Thus It Is able to see the communication 
beam from the LEO even though the scan has continued. The 10 prad beam from 
the LEO does a small spiral scan in the 22 prad pixel FOV to Illuminate the 
GEO. 

When the GEO quadrant detector recognizes that It has been illuminated, 
the GEO acquisition transmitter Is directed to stop the scan, and the 35 prad 
(Airy lobe) beacon transmitter is turned on. Due to the high power and low 
divergence of the communication beam, even though the signal is intermittent 
during the spiral scan, the GEO quadrant detector should be able to perform 
angular step reduction and locate the LEO to less than the 35 prad beacon 
beamwidth. 


Once illuminated by the beacon, the LEO tracking quadrant detector 
performs an angular step reduction to obtain closed loop fine tracking. The 
communication beam then illuminates the GEO acquisition quadrant detector 
continuously, so that closed loop tracking Is obtained, and the spot falls 
through the central hole of a bifurcated mirror onto the tracking quadrant 
detector. Once the angular step reduction is performed at the GEO tracking 
quadrant detector, both terminals have established fine tracking, and • 


communication can take place. 

The total acquisition time is given by the expression T^g = M t^ 
from Section 8.5.3, plus the handover time. The time to achieve f’ve steps of 
angular reduction at the GEO acquisition quadrant detector is determined by a 
small scan time. Since the 10 prad LEO communication beam is illuminating a 
22 prad uncertainty, there are about 3x3 cells to address. For a dwell time 


of 0.05 seconds the total time to scan 9 cells is about 0.5 seconds. 








When the 35 prad GEO beacon is returned to the LEO, the LEO tracking 
quadrant detector performs an angular step reduction reaulrlng about 0.05 
seconds. LEO fine tracking allows continuous Illumination of the GEO, where 
another 0.05 seconds Is required for GEO tracking detector angular step 
reduction. 


The total handover time to be added to T^g = M 


t^ from Section 8.5.3 


Is then 0.5 ♦ 0.05 ♦ 0.05 + 0.6 e 1.2 seconds, where the last 0.6 second term 
accounts for four time of flight delays (4 R/c). 










8.6 ISSUES FOR ACQUISITION BEAM SCANNING 


8.6.1 Beamwldth 


The beamwldth for laser communication systems Is defined as the width of 
the commonly observed far field Airy profile. The Airy profile has an angular 
Intensity dependence given by 





2J 1 (2.44* Y /a) 



2.44 tty/a 


where I Is the on axis Intensity, a is the full angle beamwldth to the 
first null of the Airy pattern and y Is the radial angular variable. 


The full beam divergence to the first null of the Airy disk is given as 


a = 2.44 x/D (2) 

e 

\ 

where x Is the wavelength and D Is the diameter of the antenna (primary 
mlrror). 


Optimization of the size of the Airy beam profile superimposed upon the 
basic pointing uncertainty area is needed to provide sufficient power at its 
boundary while also minimizing the necessary transmitter power. We will use 

o 

T ®UNC as the uncerta1nt y area t0 example a staring acquisition 
routine. It Is clear, for example, that the first dark ring of the Airy 
profile cannot be even nearly coincident with the edges of e UN( . since little 
power would be available there. We have chosen to maximize the power received 
by a collector at the edges of the pointing uncertainty region relative to the 
transmitted beamwidth, since the least Intensity Is available at the edge. 
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The optimization proceeds by using the range equation and the Airy profile. 
The range equation yields the power collected by the receiver and can be used 
to deduce an appropriate relation between a and e UN £.' The range equation Is 
effectively given by 



where P Is the power collected by the receiver, and C 0 represents the 
appropriate multiplicative factors, such as, receiving antenna diameter, 
system optical losses and reciprocal range squared. 

The optimization may be determined by setting dP(x)/da = 0 with x = 
2.44»y /a, then setting y coincident with the edge of the pointing uncertainty 
region. Thus, the optimized beamwidth a. is. 


a A * 2,08 ®UNC 



which places e^p close to the 1/e Intensity point on the Airy profile as 
can be determined from Equation 1 and Equation 2. 


8.6.2 Beam Overlap 

0 

We assume for simplicity and ease of Implementation a line by line raster 
scan of the pointing uncertainty area. Other scan geometries are certainly 
possible; such as, a spiral scan or a circular scan. However, gimbal and 
electronic requirements are simpler with a raster scan. A line width of the 
raster scan is defined by with the beamwidth of the Airy profile at its 1/e 
Intensity points (a.). The beam overlap from one line to the next is the 
fraction of the beam relative to the 1/e Intensity point which must be caused 
to traverse a portion of the previously scanned area. This is necessary to 
provide sufficient coverage of all points within the area of pointing un¬ 
certainty for the pulsed operation. The general linear beam overlap is 
denoted as c with 0 <_ c < 1 and Is defined here for the absence of residual, 
uncompensated, or unknown drift rates. 









The required beam overlap Is dependent upon whether the acquisition laser 
Is operated CW or pulsed. A laser operated CW would provide continuous 
emission along Its scan path so that c (cw) = c (cw) « 0. The situation 

a Jr 

Is different when the laser Is operated pulsed as continuous coverage Is not 
provided as the scan proceeds. The beam positions must be placed In such a 
way to preclude missed coverage. The required overlap is determined by the 
geometry of the situation, and In general, with pulsed operation, e x (p) « 
e (p) > 1/2 (1 - /2/2) s 0.15. 

A single pulse into an element of the area of uncertainty may not be 
sufficient for acquisition detection. An evaluation of the signal detection 
probability must be made to establish whether more than one laser pulse per 
area element Is necessary. The operation and the signal to noise ratio (S/N) 
required for the servo system must also be considered. The angular scan rate 
is associated with the laser pulse rate (r), and assuming m laser pulses per 
area element, 

a = (1 - E x ) .r/m (5) 

The number of dwell points (n) within the area of intltial pointing 
uncertainty Is greater than the geometrical factor (e UN « /a A ) 2 by the 
Inverse of the beam overlap, and 

n a {{ W a A )2 (1/(1 - c x > (1 - (6) 

8.6.3 Beam Scan 

The time required for the acquisition beam to fully scan the pointing 
uncertainty Is given by the number of scan elements within the area of 
pointing uncertainty, and the dwell time at each scan element. The number of 
scan elements was given by Equation 6, and the dwell time can be considered 
as t(J = Nj/B, where N, are the angular servo step reductions (quadrant 
detector) and B Is the servo bandwidth. The scan time Is then given by 








T o ■ 'W-A 1 2 < ‘ 1 ( ‘ I ",/» < 7 > 

"*x ”'y 

The values of the overlap factor when the laser Is operated CW or pulsed 
were noted earlier. 

The angular rate along a raster scan line Is associated with the beam 
divergence and the dwell time, so that 

e « a A /t d (8) 

The laser pulse rate required during the scan Is dependent upon the S/N 
required for acquisition. One the S/N Is defined for the acqulstlon sensor 
Integration period (usually t * 1/B) the number of pulses required per scan 
element may be determined. Equation 8 may then be set equal to Equation 5 to 
define the required laser pulse rate, so that, 


r s 


m 


<*-„> ‘d 


(9) 


8.6.4 Spacecraft Induced Angular Errors 

The pointing angular uncertainty area and the scan time may be modified 
by the relative angular motion between spacecraft terminals after correction 
by the spacecraft attitude and pointing control systems. This residual 
angular motion Is the uncompensated drift rate, e . This error is due to 

fluctuations and uncertainties In spacecraft attitude, ephemeris and 
electronic controls. Since the vector direction of the angular changes caused 
by e Is unknown, the impact of Its presence during the scan time may be 
Interpreted as an increase in by an additional term © . The angular 

change is given by the product of the uncompensated drift and the scan time to 

that element. The maximum value over a complete scan Is given by 


e 


m 


V e u 
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The uncompensated drift rate « u may also be viewed as causing an 
Increase In the time necessary to scan the original uncertainty area with a 
beamwidth reduced by a factor (1-c), where c Is related to © u . The 

additional beam angular overlap Is thus a A (l-c x )(l-c ). The change 

In the angular uncertainty area may be oriented along the y- dimension 
without toss of generality and used to define the additional overlap needed 
per scan line. The additional bean angular overlap Is als <;qual to the ratio 
of tne angular change e (along y) and the number of scan lines along y, 

Ny * ©ufjc/ * x )i ^ that 


and 0 < c < 1. 


1 


c ■ 


1 +(9 UNC /T o 9 u ) 


( 10 ) 


. The scan time T determined earlier will be Increased to a larger value 
(Tj) due to the overlap factor c needed to account for the uncompensated 
drift rate e u . We may find the modified scan time by equating T$ to T 0 
/(1-e), and 



T (1*T 
0 0 


V^NC* 


(id 


T- Is equal to T as defined In Equation 7 when 6=0. We recognize 
^ o u 

that Tj as defined would be different If the probability distribution of 
finding the distant satellite within were non-uniform. 
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8.7 RECOMMENDED ACQUISITION/TRACKING CONFIGURATION 


As described In Section 7.0, the quadrant avalanche photodiode detector 
provided about 19 dB more sensitivity than the PIN quadrant, and the APD quad 
Is the recommended choice for tracking detector for that reason. The APD bias 
for each quadrant must be controlled more accurately than for the PIN quad¬ 
rant, but this Is not a serious electronic burden. 

As described In Section 8.0, the APD quadrant detector cannot provide 
reasonable acquisition times, even for 100 mW laser power and no power 
margin. The CCD array has a 17.5 dB power advantage over the APD quad. Some 
form of array Is needed as the acquisition detector on the LEO, and since the 
CCD Is most sensitive it Is the recommended choice. Although adequate CCD 
technology for video cameras Is currently available, and the acquisition 
requirements are compatible with video frame rates, a development program Is 
necessary to space qualify a 400 x 500 CCD array. Since CCD's are being spac. 
quallfled for Galileo and Space Telescope, the 400 x 500 video CCD arra» 
should pose no special problems. 

A suiimary of the recommended acquisition/tracking configuration Is given 
In Table 8.7-1. It Is also suggested that NASA personnel seriously Investi¬ 
gate using the acquisition array for tracking. This could be a CCD array with 
a gyro-stabilized servo sample rate, or a CID with a pixel tracking gate to 
meer. the complementary filter servo sample rate requirements. 
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TABLE 8.7-1. RECOMMENDED ACQUISITION/TRACKING CONFIGURATION 


GEO - Acquisition transmitter 
LEO - Acquisition receiver 


GEO 

20 Inch Receive Telescope 

1) Acquisition Quad Cell 

_I_ 

Bifurcated Mirror 

2) Communication Detector 


5-10% power split 


3) Tracking quad APD < 

2.5 Inch Transmit Telescope 

1) 130 iirad (Airy lobe) acquisition beam 

2) 35 iirad (Airy lobe) beacon beam 
(share aperture by polarization) 


LEO 


6 Inch Receive Telescope 
1) Acquisition 400 x 500 CCD Array 


Bifurcated Mirror 
"1 

2) Quadrant APD for tracking and 
beacon data demodulation 



6 Inch Transmit Telescope 

1) 15 iirad (Airy lobe) 
communication beam 


» 

i 


8-71 


















8.8 SUMMARY - ACQUISITION 


The acquisition uncertainty angle was found to be dominated by attitude 
uncertainties on the TDRS spacecraft, especially due to the yaw axis. The 
resulting uncertainty was approximated by a 0.5* diameter cone. 

Acquisition sensors were investigated, including CCD, CID and photodiode 
arrays, as well as PIN and APD quadrant detectors. Commercially available 
arrays seem to be compatible with acquisition requirements, while quadrant 
detectors may require development In the areas of smaller quadrant gap, and 
smaller device size to reduce capacitance and dark current. 

Acquisition sources were evaluated. Including the single mode cw lasers 
and the high peak power Laser Diode Labs incoherent array. Although the 1.0 
kW peak power of the LDL array was appealing, It was found that Its large 
source size and divergence made It difficult to collect and collimate its 
emission. Single moda cw lasers (or coherent multi-stripe arrays. If they 
become available) are recommended as the choice for acquisition source. 

issues for quadrant detector acquisition were examined. Including field 
of view reduction and servo angular step reduction. The probability of 
successful acquisition was derived based on detection probability when 
illuminated by the transmitter, and false alarm probabilities while not 
Illuminated. The cw Integrating (CWI) approach was found to yield the best 
shot noise performance for quadrant detector acquisition, but is seriously 
limited by the DC component of optical background in a practical system. 
Sinusoidal Intensity modulation (SIM) avoids the DC background problem, but is 
about 4.5 dB worse in shot noise performance than CWI. Acquisition times for 
SIM exceed ten minutes with a quadrant APD detector and a 100 mW source, with 
no power margin. A time-gated acquisition scheme was found to have no 
advantage over SIM. The quadrant detector acquisition sequence was identified. 

The field of view reduction for array acquisition was derived, as well as 
the probability of a successful acquisition. Commercially available CCD 
arrays were found to allow acquisition times under two minutes, with a 50 mW 
source and 6.0 dB power margin. The acquisition sequence using an array was 
described. 














Acquisition times were found to be directly proportional to the 
uncertainty solid angle, and Inversely proportional to laser power, for both 
array and quadrant detector acquisition. 

Issues for scanning the acquisition transmitter beam were Identified, 
Including optimum beamwidth, beam overlap, scanning with a pulsed laser, and 
spacecraft uncompensated drift. 

A recommended acquisition/tracking configuration was proposed, which used 
an acquisition array on the LEO, and a quadrant APD for both tracking and 
beacon data demodulation on the LEO. Further Investigation Is recommended to 
evaluate performing the tracking function with the acquisition array and a 
video tracker. 
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9.0 SYSTEM WEIGHT AND POWER 


9.1 ESTIMATE RATIONALE 

% 

The weight and power estimate was made following the baseline concept of 
a two axis gimballed optical assembly with support electronics off-glmbal. 
The configuration approach was described in Section 4.0. Most component 
weight estimates are based upon earlier engineering work in this area or upon 
existing hardware developed for similar space electro-optical applications. 
The baseline concept is noted in Figure 9.1-1. The basic optical layout is 
reviewed in Figure 9.1-2. The terminal concept design parameters are 
summarized in Table 9.1-1. 

The estimate rationale is based on the use of beryllium for the material 

for the main optical elements. Both light-weight cervit and beryllium optics 

were considered. The weight of the optical elements scales with the diameter 

(D) of the elements as measured in inches. The weight of Light-weighted 

2 7 

cervit optics scales as 0.0279 D * and the weight of beryllium optics 

2 7 

scales as 0.006 D * . The weight scaling of the light weighted cervit and 
beryllium optics is shown in Figure 9.1-3. It is clear that beryllium optics 
offer substantional weight savings, especially for diameters greater than 
about 10 inches. 
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TABLE 9.1-1. TERMINAL CONCEPT DESIGN PARAMETERS 



Property 

GEO 

1 LEO 

• 

Concept 

t Separate Transmit 

and Receive Optics 

• 

• Separate Transmit 
and Receive Optics 

• 

Transmit Optics 

• 2.3 Inch Diameter 

a 

• 6 Inch Diameter 

• 

Receive Optics 

t 20 Inch Diameter 

• 6 Inch Diameter 

t 

Primary Material 

• Beryl 1iurn 

• Beryl 1iurn 

• 

Acq Sensor 

• Si-Quad 

• Si-Array I 

• 

Track Sensor 

• Si-Quad (APD) 

1 

• Si-Quad (APD) 

• 

Conn Sensor 

• Si-APD 

• SI-APD 

• 

Laser Source 

• GaAlAs 

• GaAlAS 

• 

Transmit Wavelength 

• Beacon (LDR) = 0.85 pm 

• Comm (HDR) * 0.83 pm 

t 

Receive Wavelength 

• Conn (HDR) - 0.83 pm 

• Beacon (LDR) = 0.85 pm 

• 

Gimbal 

• Two Axis-Optics on 
Gimbal 

• Two Axis - Optics on 
Gimbal 

• 

Az Range 

• ♦ 20* 

• ♦ 175* 

• 

El Range 

• ♦ 20* 

t ♦ 140° 

• 

Weight (with 
redundancy) 

211.5 lb 

175.2 lb 

• 

Power 

t 143.1 watt 

• 137.8 watt 

i; | 
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DIAMETER, INCHES 


Figure 9.1-3. Telescope Weight Scaling with Size for Be and Cervlt Material 

















9.2 WEIGHT AND POWER 


The weight estimates for the LEO and GEO terminals are noted In Table 
9.2-1 and 9.2-2. The estimate Is divided Into two columns; one without 
redundancy, one with redundancy. The LEO weight without redundancy Is 118.1 
lb while that with redundancy Is 175.2 lb. The GEO weight without redundancy 
Is 155.5 lb and that with redundancy Is 211.5 lb. The primary difference 
between a system without redundancy and one with redundancy Is reliability, 
life on orbit and risk; Since a space system Is normally desired to be long 
lived, reliable and of low risk, the higher weight estimates are more 
reasonable. Some additional weight may be expected when the system details 
regarding the spacecraft Interface are more well defined. We consider these 
weight estimates to be conservative. 

The associate power estimates for the LEO and GEO terminals are noted In 
Table 9.2-3 and 9.2-4. The estimated power for the LEO spacecraft Is 137.8 
watt and the estimated power for the GEO spacecraft is 143.1 watts. We 
caution the reader that our weight estimate for data electronics may by 
somewhat low as a result of the uncertainty regarding the modulation and data 
formatting scheme to be implemented on the (as yet) unspecified satellite. 
Additionally, redundant electronics weights. Including power supplies, are 
often optimistic early In a program development as spacecraft interfaces and 
support are not known. 

Weight and power estimates were based upon existing Hughes spacecraft 
communications systems. These included the Ku-band terminal, the shuttle 
microprocessor, the DMSP SSM/I pwave radiometer and other systems currently 
under development. The data electronics was sized by assuming that the 500 
Mbps data originated as the LEO satellite or single channel baseband digital 
unformatted form. This data was transmitted from the LEO to GEO satellite. 
It was assumed that the 500 Mbps digital data stream was fully recovered in 
baseband form on the GEO satellite. This Is somewhat more complex and more 
massive than Is required for an rf system where the baseband data is not 
recovered. 
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TABLE 9.2-1. LEO WEIGHT ESTIMATE 
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Subtotal 107.4 159.3 

Contingency (10 percent) 10.7 15.9 




































TABLE 9.2-2. GEOSYNC HEIGHT ESTIMATE 
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Subtotal 141.4 

Contingency (10 percent) 14.1 

Total 155.5 





























TABLE 9.2-3. LEO POWER ESTIMATE 


Component 


No Redundancy 
(Pounds) 


Redundancy 


OPTICS 


POINT 

AND TRACKING 


STRUCTURE 


SUPPORT 

ELECTRONICS 


THERMAL 


6 Inch dla Be telescope 
Coupling optics 
Tracking sensor 
6 Inch dla Be telescope 
Coupling optics 
Laser transmitter 

Acqulsltlon/trackliig array 
Tracking IMC (2 axis) 
iMC drive 

Point ahead IMC (2 axis) 

IMC drive 

Acquisition/tracking electronic 
Servo control electronics 
Two axis Be glmbal 
Torquers/resolvers (A z, El) 

Terminal packaging 
Spacecraft interface 
Radiation shielding 

Communication electronics 

Beacon data electronics 

Microprocessor 

Spacecraft T and C 

Power conditioning (83 percent) 

Cabling 

Torquer power supply 

Passive radiator 
Ducting 

Active control 
Thermal blankets 

Subtotal 

Contingency (10 percent) 


0.3 

6.0 

6.0 

15.0 

20.0 

12.0 


24.0 

4.0 

10.0 

3.0 

22.1 

2.5 


3.0 


130.1 

13.0 


Total 


143.1 






















•TABLE 9.2-4. LEO POWER ESTIMATE 




No Redundancy 


Area 

Component 

(Pounds) 

Redundancy 

OPTICS 

t 20 Inch dla Be telescope 

• Coupling optics 

• Communication sensor 

0.2 

a 


• 2.3 Inch dla Be telescope 

• Coupling optics 

• Acquisition/beacon sources 

a 

2.0 

• 

POINTING 

• Acquisition/tracking array 

0.3 


AND TRACKING 

• Tracking IMC (2 axis) 
t IMC drive 

• Point ahead IMC (2 axis) 

6.0 



• IMC drive 

6.0 



• Acquisition/tracking electronics 15.0 

t 


a Servo control electronics 
• Two axis Be glmbal 

20.0 | 

• 


• Torquers/resolvers (Az, El) 

12.0 


STRUCTURE 

• Terminal packaging 

• Spacecraft Interface 

• Radiation shielding 



SUPPORT 

• Communication electronics 

20.0 

• 

ELECTRONICS 

• Beacon data electronics 

4.0 

• 


• Microprocessor 

10.0 

a 


• Spacecraft CMC and TLM 

| 3.0 



a Power conditioning (83 percent) 
a Cabling 

21.3 



a Torquer power supply (83 percent) 2.5 

I 

THERMAL 

a Passive radiator 
a Ducting 
a Active control 
a Thermal blankets 

3.0 

1 

1 


Subtotal 

125.3 



Contingency (10 percent) 

12.5 



Total 

1 

137.8 

1,1 



EDSG - 75/4i - JOB and GSM 
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